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Submitted  to  the  Department  of  Physics  on  December  12,  1980 
in  partial  fulfillment  of  the  requirements  for  the  Degree  of 
t  Master  of  Science  in  Physics. 

ABSTRACT 


The  resolution  of  a  two-element  interferometer  increases 
in  proportion  to  the  element  separation  distance,  so  that 
incorporation  of  interferometer  elements  on  satellites  will  enable 
very  long  baseline  interferometry  (VLSI)  arrays  to  achieve  greater 
resolution  of  celestial  sources  than  has  been  possible  using 
elements  on  the  earth.  The  interferometric  response  is  the  visibility 
function,  which  is  the  Fourier  transform  of  the  source's  brightness 
distribution  function.  Since  a  satellite  typically  moves  faster 
in  its  orbit  than  does  a  point  on  the  rotating  earth,  and  since  the 
integration  time  period  for  cross -correlating  the  received  signals 
will  generally  be  larger  for  a  system  which  includes  satellites, 
the  interferometric  response  will  actually  be  an  averaged  visibility 
function  as  opposed  to  a  point-wise  sampled  visibility  function. 

Since  the  baseline  magnitude  is  finite,  the  visibility  function  is 
not  known  throughout  its  domain  and  is  therefore  effectively 
truncated.  Truncating  and  averaging  the  visibility  function  adversely 
affects  the  brightness  function  resolution,  and  understanding  these 
effects  is  therefore  necessary  to  more  fully  realize  the  source 
resolution  capabilities  of  an  interferometer  incorporating  elements 
on  satellites. 

This  thesis  investigates  the  effects  on  brightness  function 
resolution  when  the  visibility  function  is  truncated  and  averaged. 

The  causes  of  truncation  and  averaging  are  presented,  and  the  basic 
Fourier  transform  relation  between  the  visibility  and  brightness 
functions  is  reviewed.  Basic  properties  of  the  two-dimensional 
Fourier  transform  are  reviewed,  and  the  standard  Fourier  inversion 
method  is  presented.  This  method  is  based  on  a  generalization  of 
the  Sampling  Theorem  and  does  not  account  for  truncation  or 
averaging  effects.  The  effects  of  truncation  and  averaging  are 
illustrated  for  a  double  Gaussian  model  source  using  the  standard 
Fourier  inversion  method.  Theorems  concerning  the  Fourier  transforma¬ 
tion  of  a  truncated  and/or  averaged  function  are  then  developed. 
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and  an  algorithm  based  on  this  study  is  presented.  This  algorithm 
attempts  to  enhance  resolution  by  using  a  derived  relationship 
between  an  average  value  of  the  visibility  function  and  the  true 
source  brightness  function.  The  algorithm  is  applied  to  several 
one-dimensional  test  cases  to  illustrate  its  potential. 
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.CHAPTER  I 


INTRODUCTION 


1 .1  Problem  Statement  and  Purpose  of  This  Study 
1.1.1:  The  Role  of  VLBI 

Interferometric  techniques  in  the  radio  region  of  the 
electromagnetic  spectrum  have  been  developed  to  the  present  state  at 
which  this  technology  provides  the  most  accurate  means  of  determin¬ 
ing  positional  and  structural  information  on  extragalactic  objects 
that  are  emitters  in  this  part  of  the  spectrum.  An  interferometer 
is  an  instrument  consisting  of  two  or  more  receiving  terminals  (or 
"elements")  which  compares  signals  received  at  each  element  from  a 
common  source  in  order  to  obtain  information  on  the  nature  and 
position  of  the  source.  As  will  be  shown  later,  the  resolution  of 
an  interferometer  exceeds  that  Dossible  with  a  single  element,  and  the 
maximum  resolution  is  directly  proportional  to  the  greatest  separation 
distance  between  elements.  Present  technology  has  enabled  systems 
with  element  separations  on  the  order  of  intercontinental  distances 
to  be  used,  resulting  in  angular  resolutions  on  the  order  of  mi  1 1 i - 
arcseconds.  These  systems,  in  which  the  widely-separated  elements  are 
not  in  communication  with  one  another  during  the  actual  conduct  of  an 
experiment,  are  known  as  Very  Long  Baseline  Interferometry  (VLBI) 
sys  terns . 


The  importance  and  utility  of  VLBI  as  a  research  tool  is 
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immediately  realized  by  a  consideration  of  the  diverse  applications 
which  it  serves.  Highly  accurate  source  position  measurements  and 
capabilities  for  high  resolution  of  source  structure  are  of  great 
value  in  astrometry  and  astrophysics.  Enhanced  knowledge  of  the 
kinematics  of  celestial  bodies  and  of  diffuse  objects,  such  as  HI 
and  HI I  regions  and  molecular  clouds,  can  be  used  for  refined 
testing  of  dynamical  theories  and  to  achieve  greater  understanding  of 
astrophysical  processes  and  source  parameters  (mass,  angular  momentum, 
etc.).  Examples  of  some  applications  of  VLBI  along  these  lines 
include  resolution  of  close  but  discrete  water-vapor  masers  in  our 
own  galaxy,  accurate  positional  determinations  for  ALSEP  transmitters 
on  the  moon  yielding  information  on  lunar  kinematics,  and  tests  of 
general  relativity  by  measuring  changes  in  relative  quasar  positions 
by  deflection  of  radio  waves  in  the  sun's  gravitational  field  [Shapiro 
(1976)].  Additionally,  VLBI  may  become  an  Important  tool  in  accurate 
spacecraft  tracking  [Addleman  (1978),  Treinish  (1978)]. 

VLBI  can  also  be  applied  to  geodesy  and  to  studies  of 
dynamics  of  the  earth's  crust  [Counselman  (1976),  Whitney  (1974)]. 
Shapiro  and  Knight  (1970)  enumerate  and  discuss  the  geophysical 
applications  of  VLBI  and  indicate  the  attainable  levels  of  accuracy 
for  determinations  of  various  geophysical  parameters. 

If  a  point  source  is  observed  with  a  two-element  interferom¬ 
eter,  accurate  determination  of  the  source  position  is  predicated  on 


accurate  knowledge  of  the  baseline  vector,  which  is  the  vector  between 
the  two  elements.  We  can  therefore  see  that  accurate  knowledge  of  a 
point  source's  location  can  in  turn  be  used  to  provide  information 
on  the  baseline  vector.  In  actual  practice,  by  using  one  or  more 
different  baselines  and  observing  a  number  of  different  sources,  both 
source  positions  and  baseline  lengths  may  be  treated  as  unknowns 
and  solved  for  by  using  the  interferometry  data.  (See  Whitney  74  for 
more  details  on  this  method.)  In  this  way,  distances  between  points 
on  the  earth's  surface  will  be  ascertained  to  the  centimeter  level 
of  accuracy,  and  hence  provide  geodetic  information  and  insight  into 
associated  dynamical  processes. 

We  have  illustrated  the  importance  of  VLBI  as  a  research 
instrument,  and  seen  the  order-of-magnitude  resolution  capable  with 
present  baseline  lengths.  Given  that  Interferometer  resolution 
improves  with  increasing  baseline  length,  a  natural  Improvement  on 
existing  systems  would  be  to  incorporate  elements  on  satellites 
thereby  greatly  increasing  baseline  length  and  enhancing  resolving 
power.  However,  there  are  some  disadvantages  incurred  when  a 
satellite  is  used  as  an  interferometer  element.  In  this  thesis,  we 
shall  be  concerned  with  one  aspect  of  these  inherent  problems. 

1 .1 .2  Purpose  and  Scope 

The  purpose  and  scope  of  this  thesis  will  be  an  examination 
of  the  adverse  effects  of  truncating  and  averaging  a  function  on  the 
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resolution  of  its  conjugate  Fourier  transform.  We  will  examine  these 
effects,  and  then  develop  theorems  which  relate  the  truncated,  and/or 
averaged  function  to  its  exact  Fourier  transform.  We  will  then 
develop  an  algorithm  based  on  these  theorems  which  attempts  to  achieve 
at  least  partial  compensation  for  the  effects  of  truncating  and 
averaging  V(u,v)  and  possibly  giving  a  better  approximation  to  the 
true  brightness  distribution  function  than  does  the  standard 
Fourier  inversion  method. 

1.1.3  Adverse  Effects  on  Resolution 

Antenna  gain  is  proportional  to  the  physical  aperture  of 
an  antenna.  [Kraus  (1966)]  Since  the  radio  signals  from  celestial 
sources  are  generally  quite  weak,  large  physical  apertures  are 
normally  employed  in  the  radio  telescope  antennas  which  serve  as 
elements  in  an  interferometer  array.  To  yield  Intelligible  informa¬ 
tion  from  a  celestial  signal,  the  signal  must  not  be  less  than  the 
fluctuations  of  the  noise  which  is  present.  The  noise  results  not 
only  from  instrumental  sources  which  may  be  reduced,  but  also  from 
fundamental  causes  such  as  background  noise  in  the  sky  and  quantum- 
statistical  limitations  which  cannot  be  circumvented.  A  large 
physical  aperture  providing  high  antenna  gain  is  therefore  an 
Important  factor  in  achieving  an  acceptable  ratio  of  signal  to  rms 
noise  (signal  to  noise  ratio,  or  SNR)  for  the  system.  The  SNR  Is 
enhanced  when  the  received  signals  are  Integrated  during  the 
process  of  cross  correlation.  Increasing  the  integration  time  period 


by  a  factor  of  N  results  in  enhancing  the  SNR  by  the  factor 
[Kraus  (1966)]. 


We  have  seen  that  antenna  size  and  integration  time  are  two 
primary  parameters  affecting  the  SNR.  For  an  earth-based  interferometer, 
the  physical  aperture  may  usually  be  constructed  adequately  large  to 
keep  the  required  integration  time  period  relatively  small  while 
maintaining  adequate  SNR's.  The  integration  time  period  here  is 
considered  small  relative  to  the  time  required  for  the  projection 
of  the  baseline  vector  S  normal  to  the  source  vector  sQ  (a  unit 
vector  towards  the  source,  which  will  here  be  considered  infinitely 
distant)  to  change  by  an  appreciable  amount.  (We  shall  presently 
clarify  the  term  "appreciable".)  The  maximum  interferometer 
resolution  occurs  when  the  source  vector  is  perpendicular  to  the 
baseline  vector  and  is  directly  proportional  to  baseline  length  |S| ; 
however,  the  resolution  for  a  typical  case  in  which  the  source 
vector  is  not  perpendicular  to  the  baseline  vector  depends  upon  the 
projection  of  S  normal  to  sQ.  What  we  have  heretofore  ignored  is 
the  fact  that  $  and  sQ  are  in  relative  motion  due  to  the  earth's 
rotation  (other  motions  such  as  the  earth's  orbital  motion,  the 
motion  of  the  solar  system  barycenter,  etc.,  can  be  ignored  for  the 
infinitely  distant  source  which  we  are  idealistically  considering  here). 
Therefore,  the  resolution  of  the  interferometer  is  really  a  time- 
dependent  function. 

The  brightness  distribution  function,  B(x,y),  characterizing 
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a  source  can  be  conceived  as  the  function  which  specifies  the 
radiation  intensity  distribution  as  a  function  of  positional  coordin¬ 
ates  x  and  y  on  the  plane  of  the  sky.  (In  general,  B(x,y)  is  a 
frequency  dependent  function;  however,  we  will  consider  it  to  be 
independent  of  frequency  in  this  thesis.)  Since  B(x,y)  is  an 
intensity,  it  is  constant  along  a  ray  path  in  free  space  and  hence 
the  brightness  function  measured  at  the  earth  is  identical  to  the 
brightness  function  at  the  surface  of  the  source,  and  for  this 
reason  the  brightness  function  provides  important  information  on  the 
physical  parameters  characterizing  the  source.  As  we  will  show  later, 
the  interferometer  response  is  a  function  V,  known  as  the  visibility 
function,  which  is  the  fourier  transform  of  the  brightness  distribu¬ 
tion  function.  The  visibility  Is  a  function  of  the  interferometric 
resolution,  and  we  shall  see  that  the  arguments  of  V  are  in  fact 
two  orthogonal  resolution  components,  designated  u  and  v.  As  mentioned 
above,  the  interferometric  resolution  is  time  dependent;  hence,  u  and 
v  are  time  dependent  parameters.  The  domain  of  the  visibility  function 
Is  the  u-v  plane,  in  which  u  and  v  serve  as  orthogonal  coordinates. 

As  the  baseline  vector  S  moves  relative  to  the  source  vector  s„,  we 
know  that  the  interferometric  resolution  varies,  and  hence  we  must 
be  moving  along  a  path  in  the  u-v  plane.  We  will  see  that  these 
paths  are  generally  ellipses  in  the  u-v  plane.  As  £  moves  relative 
to  sQ  ,  the  interferometer  is  measuring  the  visibility  function 
V[u(t),  V(t)]  (implicitly  a  function  of  time  t)  by  recording  values 
of  V  along  the  track  in  the  u-v  plane  which  is  given  by  the  locus  of 
points  [u(t),  v(t)]. 
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We  can  now  see  why  it  was  important  to  have  integration 
times  relatively  short  compared  to  the  time  required  for  l&xsQ|  to 
vary  appreciably.  In  the  limit  of  zero  integration  time,  the 
interferometer  measures  V(u,v)  along  the  relevant  track  in  the  u-v 
plane.  However,  when  we  integrate,  the  interferometer  records  the 
average  value  of  V(u,v)  over  that  segment  of  the  track  which  is 
traversed  during  the  integration  time  period  tQ.  Because  B(x,y) 
is  the  Fourier  transform  of  V(u,v),  it  is  immediately  discernible  that 
the  brightness  distribution  function  is  distorted  when  observations 
with  to^0  are  made.  If  tQ  is  non-zero  but  small,  then  the  segment 
traversed  in  time  t  is  small  and  we  are  almost  measuring  V(u,v), 
since 

VCuItEWKI*'  -UJ-V/OvO 

This  case  generally  applies  to  earth-based  VLBI  systems,  and  averaging 
effects  have  previously  been  ignored.  Fomalont  (1973)  states  the 
criterion  for  short  averaging  intervals  as  being  those  cases  in  which 
Ifixs  |  changes  by  an  amount  less  than  the  antenna  radius  during  the 
integration  time  period.  However,  a  satellite  in  a  typical  orbit 
about  the  earth  moves  roughly  ten  times  faster  than  a  point  on  the 
surface  of  the  earth,  and  therefore  ^Sxs  |  undergoes  a  greater  change 
in  the  time  period  than  it  would  for  an  earth-based  VLBI  system. 

Although  there  is  no  physical  limitation  on  the  aperture 
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size  for  a  satellite  element,  practical  limitations  will  require  that 
the  antenna  aperture  sizes  be  considerably  smaller  (at  least  for 
initial  experiments)  than  earth-based  elements.  Consequently,  we 
realize  that  integration  times  will  have  to  be  correspondingly 
increased  to  achieve  an  acceptable  SNR,  which  has  the  adverse  effect 
of  making  the  averaging  considerations  for  V(u,v)  discussed  above  a 
more  significant  problem. 

To  get  an  order-of-magnitude  "feel"  for  the  Impact  of  this 
effect,  let  us  suppose  that,  all  other  factors  being  equal,  the 
satellite's  aperture  is  1%  of  that  of  a  typical  earth-based  antenna. 
The  interferometer's  aperture  Is  the  geometric  mean  of  the  two 
Individual  antenna  apertures,  or  ^A( .01A)  =  .1  A,  where  A  is  the  area 
of  the  antenna  on  the  earth.  The  signal  gain  from  the  interferometer 
therefore  decreases  by  a  factor  of  10.  The  interferometer  system 
sensitivity  is  (ideally)  related  to  the  Integration  time  tQ  by 
[Kraus  (1966)] 

sensitivity  oC/t^  . 

If  t  is  the  integration  time  period  for  an  earth-based  VLB  I  system 
and  tg  the  corresponding  period  for  our  satellite  system,  then  to 
increase  the  sensitivity  by  a  factor  of  10  to  compensate  for  the 
decreased  antenna  gain  Implies  ts=100te.  In  practice,  one  would  not 
approach  the  sensitivity  limit  of  the  interferometer  because  the 
corresponding  integration  time  period  could  well  be  of  the  order  of 
an  orbital  period  for  the  satellite.  Thus,  only  relatively  strong 
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sources  would  be  chosen  for  observation. 

Another  problem  which  we  have  not  yet  addressed  concerns 
itself  with  the  fact  that  practical  limitations  prevent  the  measure¬ 
ment  of  V(u,v)  throughtout  the  entire  u-v  plane.  This  problem  is 
independent  of  the  averaging  problem  and  is  present  in  earth-based 
systems  also.  Now,  from  the  basic  uncertainty  relation  for  conjugate 
functions  B(x)  and  V(u)  in  a  one  dimensional  Fourier  transform  pair, 
we  know  that  as  we  decrease  the  domain  of  the  function  B(x),  the 
corresponding  domain  of  the  transform  function  V(u)  necessarily  increases. 
The  implication  here  is  that  as  we  probe  into  smaller  regions  of  the 
x-y  plane,  i.e.,  as  we  examine  smaller  sources  Cor  extended 
sources  in  greater  detail)  by  using  our  satellite  VLBI  system  with 
large  |£l,  then  the  associated  visibility  function  V(u,v)  "spreads- 
out"  in  the  u-v  plane.  Limitations  on  the  extend  of  the  u-v  plane 
examined  during  the  VLBI  experiment  then  present  more  acute  problems, 
because  we  are  now  truncating  V(u,v)into  a  restricted  domain  which 
is  significantly  less  than  its  true  domain  of  definition.  To  see 
this  more  tangibly,  we  need  only  consider  that  a  point  source  has  a 
Fourier  transform  of  infinite  extent,  and  hence  any  truncation  of 
V(u,v)  due  to  practical  limitations  will  distort  the  brightness 
function  deduced  from  the  measured  V(u,v).  This  Is  equivalent  to 
saying  that  infinite  resolution  is  required  to  resolve  a  point  source. 

Since  all  celestial  sources  are  of  a  finite,  "region- 
limited"  extent  in  the  plane  of  the  sky,  a  generalization  of  the 
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sampling  (or  Nyquist)  theorem  to  two  dimensions  is  beneficial  to 
consider.  We  will  later  see  that  this  theorem  allows  us  to  merely 
sample  V(u,v)  at  distinct  points  in  a  2-D  lattice  extending 
throughout  the  u-v  plane,  and  nevertheless  reproduce  B (x,y)  exactly. 
The  Sampling  Theorem  is  the  basis  for  the  technique  of  aperture 
synthesis,  whereby  it  is  possible  to  synthetically  produce  information 
equivalent  to  that  obtainable  from  an  aperture  of  very  large 
extent  by  using  smaller  physical  antennas.  Essentially,  the  concept 
here  is  that  if  the  small  antennas  are  located  at  properly  spaced 
points  on  a  2-D  lattice,  then  their  correlated  information  is 
equivalent  to  one  hugh  antenna  whose  physical  aperture  would  be  the 
same  as  the  area  encompassed  by  the  2-D  lattice  of  antennas.  In  fact, 
it  is  sufficient  to  successively  move  a  few  antennas  throughout  the 
lattice  and  subsequently  combine  individual  recorded  results  at  a 
later  time  to  produce  the  aperture  synthesis  result  (if  the  signal 
from  the  source  is  not  time  dependent).  Earth-rotation  aperture 
synthesis  uses  array  elements  at  fixed  points  on  the  earth,  so  that 
the  earth's  rotation  causes  the  array  elements  to  move  in  some  fashion 
with  respect  to  the  source.  If  seme  of  the  proDerly  spaced  lattice 
points  are  not  occupied  by  an  array  element  during  the  observation 
period,  an  approximation  to  the  brightness  function  can  be  obtained 
from  the  data  obtained  along  the  tracks  in  the  u-v  plane  generated  by 
the  element  pairs.  Interferometer  arrays  having  enough  elements  to 
provide  good  coverage  of  the  u-v  plane  can  therefore  be  used  to 
obtain  reasonable  approximations  to  the  brightness  function.  If  an 
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array  includes  elements  on-  satellites,  then  the  averaging  problems 
discussed  above  must  be  considered.  However,  the  sampling  theorem 
alone  cannot  fully  compensate  for  the  effects  of  averaging  and 
truncation  (which  in  realistic  experiments  means  not  only  restricting  the 
extent  of  the  u-v  plane  examined,  but  also  restricting  measurements  of  V(u,v) 
to  tracks  in  the  u-v  plane  which  do  not  for  a  simple  2-D  lattice). 

1 .2  Thesis  Outline 

In  Chapter  2,  a  sample  calculation  for  a  track  in  the  u-v 
plane  for  a  two-element,  earth-satellite  interferometer  will  be 
presented.  Basic  VLBI  concepts  for  a  two  element  interferometer  will 
be  reviewed.  Basic  Fourier  transform  theorems  in  2-D  will  be 
reviewed,  and  some  examples  of  exact  transforms  given.  Truncation 
effects  are  discussed  and  summarized  by  a  theorem.  The  representation 
of  an  averaged  function  is  considered,  and  the  standard  Fourier 
transform  for  averaged  functions  is  presented.  We  then  examine  the 
sampling  theorem  and  its  relevance  to  the  2-D  visibility-brightness 
function  pair.  Aliasing  effects  and  the  effects  of  the  epoch  of 
sampling  intervals  are  considered.  In  Chapter  3,  we  apply  the 
standard  transform  method  to  a  typical  model  of  a  truncated  and 
averaged  visibility  function.  In  Chapter  4,  we  derive  theorems  for 
truncation  and  averaging  effects  and  develop  a  matrix  algorithm  for 
the  implementation  of  the  basic  theorem  for  these  effects.  The 
matrix  algorithm  is  applied  in  Chapter  5  to  several  one-dimensional 
test  cases,  and  these  results  are  compared  with  the  results  obtained 
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from  the  standard  Fourier  -inversion  method  applied  to  the  same  data 
sets  of  average  values  of  the  visibility  function.  Chapter  6 
presents  a  summary  and  the  conclusions  of  the  study.  A  definite 
conclusion  as  to  whether  the  standard  Fourier  inversion  method  or 
the  matrix  algorithm  is  the  better  technique  to  obtain  the  highest 
resolution  from  a  given  set  of  data  requires  further  analysis 
than  has  been  possible  to  perform  here.  Therefore,  the  conclusions 
presented  here  concerning  the  comparison  between  the  two  techniques 
should  be  considered  preliminary  conclusions,  subject  to  more 
extensive  investigation. 


CHAPTER  2 


BACKGROUND  AND  PRELIMINARY  CONCEPTS 

In  this  chapter,  we  will  present  a  brief  background 
discussion  on  VLBI  which  includes  a  derivation  for  the  baseline 
vector  of  an  earth-station  satellite-station  two  element  interferometer. 
We  then  consider  the  two  dimensional  Fourier  transform  and  prove 
some  fundamental  theorems  relevant  to  our  work  here.  The  standard 
Fourier  transform  method  for  an  averaged  function  is  presented.  The 
Sampling  Theorem  is  derived  for  the  two-dimensional  case  and  its 
implications  are  discussed. 

2.1  Baseline  Vector  for  the  Earth-Satellite  Two-Element  Interferometer 

The  baseline  vector  for  the  earth-satellite  interferometer 
will  be  taken  as  the  vector  from  the  earth  station  to  the  satellite 
station.  We  calculate  the  position  vector  as  a  function  of  time 
for  each  station  separately,  and  take  their  vector  difference  to 
obtain  the  baseline  vector. 

We  choose  a  rectangular  frame  of  reference  based  on  the 

standard  equatorial  coordinate  system.  The  origin  of  coordinates  is 
at  the  center  of  the  earth  with  the  z  axis  along  the  spin  axis  of 
the  earth,  the  x  axis  in  the  direction  of  the  true  vernal  equinox  of 
date,  and  the  y  axis  completing  a  right-handed  triad.  We  will 
consider  this  earth-centered  system  as  being  inertial  over  the  time 
periods  of  interest  to  us. 
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Let  the  earth  station  have  latitude  <f>E  and  longitude  ^ 
and  have  a  radial  distance  from  the  center  of  the  earth  of  rE> 
Define  the  vector  G  as  a  unit  vector  along  the  line  connecting  the 
center  of  the  earth  to  the  point  where  the  Greenwich  meridian 

A  A 

intersects  the  x-y  plane.  Let  gQ  be  the  angle  between  G  and  x  at 
the  time  t=0.  From  Figure  2.1,  it  is  easy  to  see  that  the  position 
vector  RE(t)  for  the  earth  station  at  time  t  is  given  by 

$E(t)  =  i  |rE  cos  <j>E |  cos(wt+g0+ij>E) 

+  j  lr£  cos  $E|  cos(wt+g0+^E) 

+  k  rE  sin  <t>E 

where  <d  is  the  angular  velocity  of  the  earth. 

We  consider  the  position  vector  for  the  satellite.  Viewed 
simplistically  (i.e.,  neglecting  effects  of  perturbing  bodies, 
non-sphericity  of  the  earth,  tidal  effects,  effects  of  general 
relativity,  etc.)  the  satellite  travels  about  the  earth  in  a  plane 
Keplerian  orbit.  Define  a  primed  coordinate  system  with  origin 
at  the  center  of  the  earth  and  with  z'  normal  to  the  plane  of  the 
orbit,  x'  along  the  line  from  the  center  of  the  earth  to  perogee, 
and  y'  completing  a  right-handed  triad.  (See  Figure  2.2.)  Let 
the  direction  cosines  of  the  primed  axes  in  the  unprimed  system 
be  ( av  i  v  9  Bv  i  ..  9  yv  i  7) »  (ciy  i  v  ,  6U i  y  »  Yw i  7) »  and 

^  x  jjr  a  ^  y  ,a  j  j 

(“z’x  ’  Bz',y  ’V.z^  for  x'»  y’s  and  z’’  resPectively-  Gl'ven  an 
initial  position  and  an  initial  velocity  for  the  satellite  attimet=0, 

one  may  compute  the  position  vector  ^’^(t)  for  the  satellite  at  any 

time  to.  (This  is  a  standard  wel 1 -documented  problem  in  astro- 
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FIGURE  2.2  —  Satellite  Orbit  Plane 
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dynamics  and  we  will  not  present  the  full  general  solution  here. 
The  Fortran  algorithm  which  was  used  to  solve  for  on  the 

computer  was  based  on  a  similar  algorithm  written  by  T.  Herring.) 
The  conversion  of  ^'$(t)  to  our  unprimed  system  is  given  by 


+  -5  (ft*  ft.;.  + ft*  k,)  +  V?(Rx-V«:>  +■  E«  £») 


where  ft's(t)  =  [Rx , ,  Ryl,  Rz , ] . 

The  baseline  vector  as  a  function  of  time  is  then  given 

by 

S(t)  =  ts{t)  -  T5E(t) 


2.2  Basic  VLBI  Concepts 

A  detailed  description  and  analysis  of  a  VLBI  system 

would  involve  a  far  more  extensive  study  than  is  possible  here  and, 
in  any  event,  there  are  extensive  references  in  the  literature 
discussing  the  intricacies  of  VLBI.  We,  therefore,  present  a 
simplified  view  of  VLBI  which  will  be  sufficient  as  background  and 
motivation  for  the  specific  problem  studied  in  this  thesis. 
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The  two  element  interferometer,  which  is  the  building  block 
for  multi -element  arrays,  consists  of  two  receiving  antennas  and 
equipment  to  process  and  correlate  their  received  signals.  If  a  real¬ 
time  link  such  as  a  cable  is  used  to  connect  the  two  elements, 
then  real-time  data  analysis  is  possible  and,  more  significantly, 
if  the  electrical  path  lengths  are  properly  set  then  the  relative 
phase  information  for  the  received  signals  is  directly  input  to  the 
correlator.  In  VLBI,  a  real-time  link  between  elements  is  usually 
not  possible  due  to  the  very  large  baseline  distance.  Therefore, 
real-time  data  analysis  is  not  carried  out,  and  instead  the  data 
are  recorded  for  later  processing.  To  preserve  the  relative  phase 
information,  which  is  essential  for  correlation,  to  an  acceptable 
level  of  accuracy  conmensurate  with  the  desired  level  of  accuracy 
expected  for  the  parameters  of  interest  now  requires  highly  accurate 
timekeeping  during  the  conduct  of  the  independent  observations.  The 
advent  of  hydrogen  maser  frequency  standards  provided  the  time¬ 
keeping  accuracy  required  to  enable  VLBI  to  exceed  the  accuracy 
levels  attainable  by  connected  element  interferometry  [Whitney  (1974)]. 
With  accurate  time-keeping  at  each  element,  the  received  signals 
can  be  recorded  on  magnetic  tape  for  later  processing.  In  fact, 
when  the  data  are  correlated,  it  is  possible  to  ascertain  the 
error  in  original  clock  synchronization  between  the  two  elements 
and  thereby  provide  a  technique  for  unprecedented  accuracy  (to 
about  .1  nanosecond)  in  clock  synchronization  over  intercontinental 
distances  [Shapiro  and  Knight  (1970)].  The  actual  reduction  and 
processing  of  data  from  a  VLBI  experiment  is  an  extensive  subject 


in  itself  and  cannot  be  discussed  here.  There  are  several  different 
algorithmic  approaches  to  the  problem  of  data  analysis,  and  ample 
references  may  be  found  in  the  literature  [see,  for  Instance, 

Whitney  (1974),  and  Moran  (1976)]. 

The  basic  geometry  for  a  two-element  interferometer  is 
illustrated  in  Figure  2.3.  We  assume  an  infinitely  distant  point 
source  so  that  plane  waves  are  received  by  the  interferometer 
elements,  which  therefore  point  in  the  same  direction  at  angle  4 > 
with  respect  to  the  baseline  vector  S  (i.e.,  we  neglect  parallax 
effects  which  would  have  to  be  accounted  for  if  the  source  were 
sufficiently  close  to  the  interferometer  to  cause  different  antenna 
orientations).  The  baseline  vector  t  points  from  the  reference 
station  to  the  remote  station,  where  the  reference  station  receives 
any  given  wavefront  earlier  than  the  remote  station  by  a  time 
factor  known  as  the  group  delay  t  *J$*s  .  The  projection  of  the 
baseline  normal  to  the  source  is  just  jbxs0|,  and  will  be  seen  to 
be  related  to  the  interferometric  resolution. 

Our  aim  at  this  stage  is  to  demonstrate  that  the  response 
of  the  interferometer  is  the  visibility  function,  which  is  the 
Fourier  transform  of  the  source  brightness  distribution  function. 

A  detailed  derivation  requires  the  application  of  coherence 
theory  to  the  two  element  interferometer  [see,  for  instance,  Swenson 
and  Mathur  (1968)].  However,  a  less  sophisticated  development 
will  suffice  for  our  purposes  here.  We  will  base  our  discussion 
on  the  treatments  presented  by  Moran  (1976),  Fomalont  (1973),  and 
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Whitney  (1974).  Partial  coherence  theory  establishes  that  all 
characteristics  of  a  source's  radiation  power  idstribution  (i.e., 
angular  distribution,  frequency  distribution,  and  polarization)  are 
related  to  spatial,  temporal,  and  polarization  cross  correlations 
of  the  received  signals  at  the  two  interferometer  elements.  We 
will  primarily  deal  only  with  the  spatial  correlation  aspects  here, 
which  assumes  monochromatic  plane  polarized  radiation  from  our 
infinitely  distant  point  source. 


Let  the  point  source  emit  a  sinusoidal  signal  at  frequency 
vQ  with  intensity  BQ.  Denoting  e-j  and  e£  as  the  received  voltages 
at  the  reference  and  remote  stations,  respectively,  we  have,  using 
complex  phasor  notation, 

iJLTT&t 

c. 


e- 


uirMt-Vi) 


The  interferometer  correlates  these  two  received  voltages.  [See 
Rogers  (1976)  and  Kraus  (1966)  for  discussions  of  various  techniques 
of  detection  and  concomitant  correlation  method.  We  will  assume  a 
simple  phase-switched,  or  multiplying,  interferometer.] 

Thecross  correlation  of  two  real  complex  functions  in 
the  Interval  0,tQ]  is  defined  by  [Bracewell  (1978)]: 
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For  our  interferometer,  then,  the  cross-correlation  function  for  the 
signals  received  in  the  time  period  [0,tQ]  is 

Ct)  kj)  B. 

where  superscript  "*"  denotes  complex  conjugation.  The  time  period 
tQ  is  referred  to  as  the  integration  time  period.  R(-rg)  is  the 
response  of  the  interferometer  for  the  point  source.  Now, 

AV  -S.  =  V)  co»Cplt> 

where  we  have  made  the  angle  $  between  &  and  sQ  a  time  dependent 
parameter  because  the  relative  orientations  of  S  and  sQ  are  changing 
as  a  result  of  the  earth's  rotation  (and  satellite  motion  for  a 
satellite  element  in  the  two  element  interferometer).  We  can  then 
write  ,  _ 

We  see  that  R( <{>)  is  an  oscillatory  quantity  with  a  mean  of  zero, 
amplitude  of  B  ,  and  time-varying  phase  factor  of  2ir~cos<t>(t) . 

An  alternative  description  for  the  slow  oscillations  in  R(<^>)  derives 
the  cross-correlation  function  in  terms  of  beating  between  the  two 
signal  frequencies  from  the  two  elements,  whose  difference  (in 
general  from  vQ  and  from  each  other)  arises  from  the  differential 
Doppler  shift  in  the  received  frequencies  due  to  the  relative 
motion  of  the  elements  with  respect  to  the  source.  [See  Rogers  (1976) 
for  this  approach.] 
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In  terms  of  4>,  we  may  think  of  R(<j>)  as  a  spatial 
correlation  function.  The  spatial  correlation  of  the  incident  radia¬ 
tion  field  is  determined  by  the  correlation  of  the  signals  from  the 
two  elements  located  at  different  points  in  space,  and  we  see  from 
Figure  2.3  and  from  the  expression  for  R(<j>)  that  the  two  signal 
phases  will  be  mutually  reinforcing  (or  partially  reinforcing)  or 
cancelling  depending  upon  the  magnitude  of  cos<t>,  i.e.,  depending 
upon  the  difference  in  phase  path  length  for  each  wavefront  reaching 
the  two  elements.  (We  have  tacitly  assumed  the  limit  T-*0  here  so 
that  the  spatial  correlation  is  performed  point-wise  in  space  as 
opposed  to  averaging  over  segments  of  space  for  non-zero  integration 
time  periods.)  This  phenomenon  is  directly  analogous  to  the 
fringe  patterns  of  optical  Interferometry,  and  R(<t>)  is  known  as  the 
fringe  pattern. 


The  phase  factor  2  ^  cos<f>  in  R(<<>)  is  known  as  the  fringe 
phase  $(t) .  R( <j>)  assumes  its  maximum  value  Bq  when  ^  cos<j>  is 
equal  to  an  Integer.  The  interval  between  successive  maxima  in 
R(<t>)  denoted  A4> ,  is  known  as  the  fringe  spacing.  Let  <j>-|  and 
be  two  values  of  $(t)  which  produce  adjacent  maxima  in  R($>).  Then 

cos ^  *-2’JT(~n+i)-2T~n  *  27T 

^  <f>a  -  -=-•»  <?, )  *  (-Z  wt  )  =  | 


<f>2+lh 

Assuming  that  we  write  — ^ .<f>2  and 

sin  -  *  ^2~  -  ^  •  Then,  -2^-  si n<}>  =  1,  which  gives  the 

magnitude  of  the  fringe  spacing  as  a<j>  =  j^yn—  But,  bsl.p(J)  =  jbxsj 


which  is  the  component  of  b  normal  to  the  souce  vector  s0- 
We  therefore  see  that  the  fringephase  provides  information  about  the 


source  location  s0. 


Having  obtained  some  understanding  of  an  idealistic 

interferometer  response  for  an  infinitely  distance  point  source 

emitting  monochromatic  radiation,  we  now  want  to  look  at  the  more 

realistic  case  where  an  extended  source  with  centroid  approximately 

at  s„  is  observed  in  the  bandwidth  Av  centc.cu  on  \>  by  an 
o  o 

interferometer  with  frequency  characteristic  F(v)  and  primary  power 
pattern  G(s-sQ)  (defined  as  the  product  of  the  two  elements' 

A  A  A  A 

voltage  patterns).  Heretofore,  we  had  assumed  G(s-s0)=6(s-s0) 
and  F(v)=s( v-vQ) . 


First,  we  note  that  if  @s  is  the  maximum  angular  source 
extent,  then  we  must  have 


A<f 


T— A —  >  0 
b  s  i  n<j>  s 


since  a  source  larger  than  the  fringe  spaceing  will  have  simultaneous 
reinforcements  and  cancellations  from  its  various  parts  and 
therefore  the  meaningful  fringe  pattern  is  lost.  The  amplitude  of 
the  fringes  for  an  extended  source  will,  in  general,  be  different 
from  the  amplitude  for  a  point  source  of  the  same  strength,  and  the 
ratio  of  these  amplitudes  is  known  as  the  fringe amplitude. 

A 

The  reason  for  a  specified  reference  direction  sQ  for  an 
extended  source  is  seen  in  the  following  analysis.  For  observation 
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of  a  point  source  over  bandwidth  Av  [but  assuming  an  infinitely 
narrow  beam  for  the  primary  power  pattern,  G(s-sQ)  =  6(s-sQ)]  we 

^  ,  /•«•+£  -b  i* 


-  t?  ^7T-fer  V$> 

“  tt  Jfc.  ik>  ...» <t>  e 


where  we  have  assumed  a  flat  frequency  response 

f  |  ^  c  L  C V.-  40/a  >  ( vV»*  4v>)/aJ 


/SM- 


and  included^  as  a  normalization  factor  to  compare  this  result  with 
the  monochromatic  result  above.  We  see  that  this  result  is 
identical  with  the  monochromatic  result  except  for  the  factor 

-2jp-  )  ?  e  -rr-b  «f  s  7rAV’/tj- 

which  is  a  taper inq  factor  which  reduces  the  fringe  pattern  amplitude. 

This  presents  the  conflicting  goals  of  using  Tg«^  (i.e.,  the 

signal  correlation  time  must  be  much  greater  than  the  group 

delay  time)  in  order  that  the  fringe  amplitude  be  sufficiently 

large,  and  the  other  goal  of  using  as  large  a  bandwidth  as  possible 

to  be  able  to  detect  weak  sources  or  source  components  without 

excessively  long  integration  times  (since  sensitivity  a  /==£=) . 

'Av 

We  can  achieve  both  of  these  goals  if  we  can  keep  Tg  very  small. 

A 

To  do  this,  we  define  a  reference  direction  sQ  to  be  roughly  in  the 
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direction  of  the  point  source  or  the  direction  of  the  centroid  of  an 
extended  source.  A  delay  factor  is  then  incorporated  into  the 
reference  station  arm  of  the  interferometer  by  using,  say  an  extra 
length  of  cable,  such  that  Tg=0  when  the  interferometer  elements 
are  pointed  in  the  direction  sQ.  Then,  Tg  will  be  small  for  all 
s  near  sQ,  which  is  usually  the  case  for  observed  celestial  sources. 
In  VLBI,  this  delay  factor  can  be  introduced  by  an  appropriate 
delay  synchronization  between  the  magnetic  tapes  from  the  two 
elements  when  they  are  brought  together  and  played  back  for  correla¬ 
tion,  or  equivalently  by  introducing  a  relative  shift  in  the  bit 
strings  of  data  during  cross  correlation  for  digital  systems.  With 

A 

the  delay  factor,  the  fringe  phase  is  set  to  zero  at  s0,  so  the 
fringe  amplitude  is  maximum  there  and  falls  off  very  slowly  due  to 
the  factor  for  s  slightly  offset  from  sQ.  We  note  that  the 

delay  factor  can  also  be  tracked  to  compensate  for  changing  Xg 
and  thereby  allow  <5.v  to  be  even  larger  than  for  fixed  delay. 


The  appropriate  generalization  for  the  correlation  function 
giving  the  interferometer  response  for  a  point  source  observed  in 

bandwidth  Av  with  accurate  delay  tracking  is 


'R(/,syBe,  e 


General izaing  this  to  the  case  of  an  extended  source  with  brightness 
distribution  function  B(s-sQ),  we  have 


I 
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R  = 


ultf  c 
e 


u- 

iZTT  cT  ® 


=  e 


Vet) 


where  the  function  V(S)  is  defined  to  be 


V(F)  = 


=  X  -  #  Fm  &cs-k)  e 


^TT^C^S-S^ 


We  can  interpret  this  result  as  follows.  The  factor  e 
is  just  the  interferometer  response  to  a  monochromatic  point  source 
with  unity  intensity.  This  factor  has  a  relatively  fast  oscillation 
rate  corresponding  to  closely  spaced  fringes.  The  complex  function 
V(F)  gives  the  amplitude  and  phase  offset  needed  to  properly 
characterize  the  extended,  non-monochromatic  source  which  is 
actually  observed.  Oscillations  in  V(F)  are  large  compared  to  the 
fringe  spacing  (or,  equivalently  in  terms  of  a  temporal  parameter, 
»(8)  oscillates  slower  than  do  the  fringes  as  the  baseline  rotates 
relative  to  sQ). 

We  can  express  the  result  in  a  more  convenient  form  in 

terms  of  a  suitable  astrometric  coordinate  system.  First,  we  note 

that  the  vector  s  can  be  expressed  by  an  angle  x  in  right  ascension 

and  y  in  declination  relative  to  the  reference  direction  s„. 

o 

Alternatively,  but  equivalently,  we  can  think  in  terms  of  a 
rectangular  coordinate  system  on  the  plane  of  the  sky,  with  origin 
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at  the  point  where  sQ  intersects  the  plane  of  the  sky,  and  with 
unit  vectors  x  and  y  with  ground  projections  in  the  east  and  north 
directions,  respectively.  We  will  adopt  the  latter  terminology 
here.  The  appropriate  units  for  x  and  y  are  radians,  since  the 
angular  description  is  the  operationally  correct  one,  although  the 
planar  description  makes  it  easier  to  conceptualize  brightness 
distribution  functions  on  the  sky.  (We  will  not  be  concerned  with  the 
longitudinal  component  here,  which  is  in  the  direction  sq=x  x  y.) 

With  this  notation,  G(s-sQ)  can  be  written  G(x,y),  B(s-§0)  as  B(x,y), 
ds  as  (dx-dy),  and  (s-so)  as  xx  +  yy.  For  monochromatic  radiation, 
we  may  now  write 


v(S)  = 


s-w,y) 


e 


Considerations  of  diffraction  limiting  effects  in  one 
dimension  show  that  the  maximum  resolution  of  a  two-element  inter¬ 
ferometer  is  -J (to  within  a  factor  of  order  unity),  which  occurs 
when  the  baseline  is  normal  to  sQ.  Clearly,  infinite  resolution 
is  possible  only  for  an  infinitely  long  baseline.  Infinite 
resolution  corresponds  to  zero  fringe  spacing,  and  the  resolution 
decreases  as  the  spacing  between  fringes  increases.  In  the  two- 
dimensional  problem  wich  we  are  considering  here,  we  must  deal  with 
two  orthogonal  components  of  resolution.  We  use  the  9  and  y  directions, 
and  we  define  unit  vectors  u=x  and  v=y  in  a  u-v  plane  which  is 
parallel  to  the  x-y  plane.  The  u  and  v  axes  then  give  the  orthogonal 
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resolution  components  in  the  east-west  and  north-south  directions, 
respectively.  The  units  of  u  and  v  are  the  same  and  are  reciprocals 
of  the  units  of  x  and  y;  for  x  and  y  in  radians,  u  and  v  are  in 


fringes  per  radian. 


t  . 


We  can  express  —  in  terms  of  the  resolution  components: 


?-  =  ux  +  vy  +  (Mj 


We  will  neglect  the  last  term  which  gives  the  longitudinal  component, 
and  consider  only  the  transverse  component  Sp  =  ux  +  vy.  In  terms 
of  the  standard  earth-based  coordinates  of  right  ascension  a  and 
declination  6,  we  write  (aQ,  <50)  for  the  angular  coordinates  of 
sQ  and  (afa,<5b)  for  b.  The  transverse  component  of  b,  which  is 
normal  to  the  source  direction  sQ,  can  then  be  written  as 

t  =  C  = 


+  9  -  Sn*Ss  Co*Sk  e*s(PV‘°0] 


We  identify 


u  =  W°V°0 


*  -fc  [*"*4.  °*s4 _  ~s£k  ces  K-*J] 


f 
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Now,  expressing  V(S)  in  terms  of  u  and  v,  and  noting  that 

lit  (KH+'W 


(s-sQ)p  =  xx  +  yy,  we  can  write 


V( 


u,v)  =  _/<U.Jiy  Gr  c^y)  e 


The  function  V(u,v)  is  know  as  the  visibility  function,  and  for  an 
isotropic  primary  power  pattern,  we  have 


V(u,v) 


*7 


Thus,  the  visibility  function  and  the  source  brightness  distribution 
function  are  a  Fourier  transform  pair.  Fomalont  (1973)  discusses 
the  major  assumptions  and  simplifications  which  are  incorporated 
into  this  result  and  its  inherent  limitations,  as  well  as  the  effect 
of  a  large  bandwidth  Av. 

We  note  that  u  and  v  are  time-varying  parameters  due  to 
the  relative  motion  between  £  and  sQ.  For  an  interferometer  with 
two  earth-based  elements,  the  locus  of  points  in  the  u-v  plane 

trans versed  by  u  and  v  will  lie  on  an  elliptical  tract,  since 

ul  +  Co--  op*  _  . 
of  5^  '  1 

where 

«  5  C“  <»  ct>s  Sfa  *"£3 

a"d  »o  s  -^7  <Sfc  “S  <Cs 
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However,  when  one  or  both  elements  of  an  interferometer  are  on 
satellites,  the  periods  of  the  individual  elements  (which  for  two 
earth-based  elements  are  both  24  hours)  are  no  longer  necessarily 
commensurate,  and  it  can  be  shown  that  the  track  in  the  u-v  plane 
is  not  a  closed  figure,  but  rather  is  "open-ended"  and  undergoes  a 
precessional-type  advance  through  the  u-v  plane.  Thus,  a  two 
element  interferometer  whose  elements  have  non-commensurate  periodic 
motions  can  (in  principle,  with  infinite  observation  time)  provide 
observational  coverage  over  the  entire  u-v  plane  inside  the  trunca¬ 
tion  limits  (which  are  the  maximun  resolution  limits  along  the  u 
and  v  axes) . 

Figure  2.4  illustrates  a  portion  of  the  track  which  is 
traced  out  in  the  u-v  plane  for  the  case  of  an  earth  station  located 
at  43°  latitude  and  70°  longitude  and  a  satellite  element  in  a 
circular  orbit  over  the  poles  at  distance  8-10-9 cm  from  the 
center  of  the  earth  (1.62-10^km  altitude).  The  source  coordinates 
used  were  those  for  3C273,  right  ascension  12H  26m  33s  and  declination 
2°20'  [Kraus  (1966)],  and  the  observation  wavelength  was  taken 
as  3cm.  The  only  portions  of  the  u-v  plane  which  may  actually 
be  observed  are  those  portions  where  the  source  is  visible  to 
both  elements.  From  the  figure,  we  can  get  an  appreciation  for  the 
role  of  the  integration  time  period  in  terms  of  the  extent  to  which 
the  visibility  function  is  averaged  as  opposed  being  sampled  (i.e., 
measured  point-wise).  The  integration  time  t  was  taken  as  1000 
seconds  for  the  example  case  here.  As  discussed  previously,  the 


FIGURE  2.4  —  Track  in  the  u-v  plane. 


Note:  The  inteval  between  markers  corresponds  to  a  1000  second 
integration  time.  The  start  of  this  portion  of  the  track 
is  the  end  near  the  origin. 
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integration  time  period  is  of  the  order  of  100  times  greater  for  an 
earth-satellite  interferometer  as  compared  to  an  earth-based  system. 

The  effects  of  averaging  must  therefore  be  accounted  for  when  a  satellite 
is  used  as  an  interferometer  element.  A  useful  simplification  arises 
if  we  shift  axes  so  that  vQ=0  or  if  the  data  are  "gridded"  (see  discus¬ 
sion  below)  onto  a  aregular  lattice.  In  this  case,  then,  the  arc  or  cell 
midpoint  at  (un,vm)  has  the  same  projection  magnitudes  on  the  u  and  v 
axes  as  does  the  arc  or  cell  with  midpoint  at  (u  „,v  „) ,  where  u  =-u„, 
v_m--vm-  If  we  can  call  these  projections  Tu  n  and  Ty  m,  respectively, 

then  we  can  write  T  =T  _  and  T  =T ,  This  simplification  will 
u,n  u,-n  v,m  v,-m  r 

be  used  later  in  our  discussion  of  averaging  effects. 

Essentially,  we  have  found  that  the  interferometer  measures 
the  Fourier  components  of  the  source  brightness  distribution.  As  we 
will  prove  later,  the  Sampling  Theorem  states  that  knowing  the 
Fourier  components  of  a  function  at  discrete,  properly  spaced  lattice 
points  in  the  u-v  plane  enables  one  to  completely  specify  the 
brightness  function.  This  concept  is  the  basis  for  aperture 
synethesis  techniques,  in  which  B(x,y)  is  determined  from  incomplete 
sampling  of  the  visibility  function  in  the  u-v  plane.  Since  the 
coverage  of  the  u-v  plane  is  confined  to  the  tracks  made  as  the 
relative  orientation  of  b  and  sQ  changes,  it  is  desirable  to  use 
many-element  arrays  with  various  baseline  vectors  so  that  ample 
coverage  of  the  u-v  plane  results  from  an  observation.  Each  pair 
of  elements  In  this  array  produces  a  track  in  the  u-v  plane.  If  the 
signal  from  the  source  is  not  a  rapidly  time-varying  phenomenon 


(such  as  the  solar  brightness  function),  or  is  approximately  a  steady- 
state  phenomenon  over  the  period  of  all  observations,  then  the 
various  observations  to  be  incorporated  into  the  aperture  synthesis 
map  need  not  be  made  simultaneously.  This  fact  presents  the 
added  advantage  of  allowing  one  to  use  a  certain  number  of  elements 
for  a  series  of  observations,  where  the  baseline  vectors  are 
changed  (say  by  moving  elements  on  railroad  platforms)  from  one 
observation  period  to  the  next.  In  this  manner,  it  has  been  possible 
to  obtain  aperture  synthesis  maps  which  result  from  excellent 
coverage  of  the  u-v  plane.  Although  the  density  of  observation 
points  may  be  reasonably  high,  there  still  remains  the  problem 
that  the  tracks  in  the  u-v  plane  do  not  (in  general)  pass  through 
the  properly  spaced  lattice  points  required  by  the  Sampling 
Theorem  to  yield  the  brightness  function  exactly  and  uniquely. 

Various  methods  are  in  use  for  extrapolating  (or  "grldding")  the 
known  data  to  the  lattice  points,  including  convolution  and 
cell  averaging  techniques  [see  Thompson  and  Bracewell  (1974) 
for  a  discussion  of  these  methods].  An  additional  problem  concerns 
the  fact  that,  although  good  coverage  of  the  u-v  plane  may  be 
obtained,  full  coverage  (i.e.,  infinite  resolution  from  infinite 
baselines)  is  not  possible  and  effectively  we  are  dealing  with  a 
truncated  version  of  the  true  visibility  function. 

Aperture  synthesis  theory  applies  to  VLBI.  In  a  particular 
case,  if  it  is  not  possible  to  obtain  sufficient  coverage  of  the 
u-v  plane  to  enable  a  high  resolution  brightness  distribution  map 


to  be  constructed,  then  a  simple  model  source,  commensurate  with 
the  known  data,  may  be  constructed  to  represent  the  true  source 


[see,  for  example,  Knight,  et  al  (1971)].  With  a  sufficient  number 
of  ground  based  stations,  and  by  using  satellites  as  Interferometer 
elements,  it  should  be  possible  to  obtain  good  coverage  of  the  u-v 
plane  to  enable  very  high  resolution  aperture  synthesis  maps  to 
be  produced. 

Having  demonstrated  that  the  relationship  between  the 
visibility  function  and  the  source  brightness  function  is  a 
Fourier  transformation,  we  now  review  the  basics  of  Fourier 
transform  theory,  and  subsequently,  we  will  study  the  effect  on 
the  brightness  function  due  to  truncating  and  averaging  the 
visibility  function. 

2.3  The  Two-Dimensional  Fourier  Transform 

2.3.1  Usefulness  of  the  Fourier  Transform  as  an  Integral  Transform 

Of  the  different  possible  integral  transforms,  the  Fourier 
transform,  because  of  its  basic  properties.  Is  one  of  the  most  useful 
and  widely  employed.  These  basic  properties  are  linearity,  shift 
invariance,  and  orthogonality  of  the  integral  transform's  kernel. 

If  we  denote  the  Fourier  transform  operation  by  the 
operator  symbol  &  ,  then 

■?  {&*>}  =  N/(u) 


I 


42 

means  that  the  Fourier  transform  of  a  function  B  of  variable  x  gives 
a  function  V  of  the  conjugate  variable  u.  In  terms  of  units,  u 
and  x  are  reciprocals;  for  example,  if  x  is  in  centimeters,  then 
u  is  in  inverse  centimeters. 

The  property  of  linearity  means  that 

0}  *0<N/,CU)  4-  £ 


where 

Shift  invariance  implies,  for  a  constant  c, 

-ftcWc*) 

where  :f --  \j^) 


and  f(c)  is  a  phase  factor,  depending  upon  constant  c  and  u,  but 
independent  of  x.  The  orthogonality  property  enables  the  inverse 
Fourier  transform,  denoted  by  the  operator*^  -1 ,  to  be  defined. 
Thus  if,  ^f[dcx)~j  = 


then 


In  one  dimension,  the  Fourier  kernel  is  written  as  e^*ux,  so  that 


the  Integral  transform  appears  explicitly  as 


VCO*  e^*T5oo  4* 


I 
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The  orthogonality  of  the  Fourier  kernel  then  lets  us  write 


.  .Lvn**  . 

4m.  C  Vtw'i 


'V-.Ju6  U— *  "Botxlxj 

-r  t a  r° ,  uttu.(k*x') 

J-«,cU  du.  6. 

=/ld*  BoO  “  "BCx') 


We,  therefore,  see  that  the  orthogonality  property  lelfs  define  the 
inverse  transform  kernel  as  e~^1TUX,  and  the  inverse  transformation 
appears  explicity  as 

I3oo*  dn  e.  Vcm) 


It  is  the  importance  of  linearity,  shift  invariance,  and 
inverse  transformation  operations  in  practical  applications  which 
accounts  for  the  diverse  and  widespread  use  of  Fourier  transforms. 


2.3.2  The  2-0  Fourier  Transform 


Multi -dimensional  Fourier  transforms  are  defined 
analogously  to  the  one-dimensional  form,  with  the  number  of 
conjugate  variable  pairs  being  equal  to  the  dimension.  In  particular, 


the  2-D  Fourier  transforms  are  given  by 

■Bu.v)*T’K«-4=  £X£d*-  e 


~L*irCuJIL+'ry) 


U.ir) 


-Li,  .'“’“jCa. 
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We  will  now  consider  a  few  examples  of  the  2-D  Fourier  transform 
which  will  prove  relevant  in  our  discussion  below. 

A  point  source  at  the  origin  is  represented  by  B(x,y)=<$(x,y), 
where  the  2-D  Dirac  delta  function  Is  defined  by  <5(x,y)=a(x)6(y) . 

The  corresponding  visibility  function  would  be 

\/ci*,u-)=J_«dxJ^.cla  e.  6(x,y)  =  J- 


Thus,  a  point  source  has  an  associated  visibility  function  of 
constant  amplitude  over  the  entire  u-v  plane.  If  the  visibility 
function  were  of  constant  amplitude  over  only  a  finite  rectangular 


area  of  the  u-v  plane 

N/Kc). 


1  Fa *  tr VT] 


ELSE  WMfrtT 


the  correspondinn  source  brightness  function  would  be 

•a. . -  f~.  r.  _-i*irCnx+uyK/, 


V(«,vr) 


If  _£«nc«f  -12-nr.y 


e - ,  e. 

S)M  2lT0o\ 


iZTTV. 


iZTH 


The  sine  function  divided  by  its  argument  will  appear  frequently 
in  our  work,  and  will  be  called  a  sine  function  following  Bracewell's 
convention:  sine  x  =  s- Thus,  the  previous  brightness 
function  may  be  written  as 


As  a  final  example,  we  will  consider  a  Gaussian  visibility 
function,  _ 


— TT  (  U,  +  v  ) 


r*  r°  T(u*-+V%) 

Then  B  CX,  y)  ^  C. 

I 


Consider  the  u  integral.  Completing  the  square  of  the  exponent, 

7T^V(^7nc)U.=  (>J1Tu  +  l  firx^TT*’ 

-tfrU.,'  +  C«TTU.Xj 


let  us  wri 


«x,jc; 


tin.  £* 


-7TX 


=  c  -jGk  e 


-[*ru.+lJir<r 


Define 


{  S  -fifK  +  iV5F*  ,  So 

if*  Y5r  dn  . 

c  &  iE‘= 


o'™ 


(l^) 


-7Txa 


Then, 


46 


So  we  can  see  that  the  corresponding  brightness  function  to  a 
Gaussian  visibility  function  is 

.  -irCx. 

e. 

It  should  be  noted  that  the  Gaussian  is  the  only  non-general i zed 
function  which  preserves  its  function  form  under  Fourier 
transformation.  (The  generalized  function  III,  discussed  later, 
will  also  be  seen  to  have  this  property.) 

2.3.3  Fundamental  Theorems 

We  have  seen  that  even  a  simple  Gaussian  requires 
some  effort  to  perform  the  Fourier  transformation  by  direct 
evaluation  of  the  integral.  Fortunately,  there  are  a  number  of 
fundamental  theorems  which  not  only  ease  the  calculation  of 
transforms,  but  also  provide  further  insight  into  properties  of 
Fourier  transforms.  Since  we  will  extensively  rely  on  these 
theorems  in  what  follows,  we  will  prove  a  number  of  them  here. 
[Most  of  these  theorems  are  stated  without  proof  for  the  2-D 
Fourier  transform  by  Bracewell  (1978).]  The  symbol  "3"  will 
be  used  to  denote  a  Fourier  transform  pair,  so  that 

N/ch,^  D  B  (.*,*3)  /W* 

Note  that  one  dimensional  versions  of  these  theorems  are  easily 
realized  by  setting  y  and  v  to  constants  in  the  equations.  (Of 
course,  these  arbitrary  constants  do  not  affect  the  results.) 


I 
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Separable  Product  Theorem:  If  B-j(x)DVi(u)  and  B2(y) 
3V2(v),  then  B1(x)B2(y)  3  V-,(u)V2(v). 

Proof:  Let  V(u,v)  =  V-j(u)V2(v)  and  B(x,y)  =  B-|(x)B2(y). 

Then  {Vcu.ir)]  e,'**"*  V.CuWxC^) 

.[&*  e“”“v.<4£>  e  a~VJ 

B^ly) 


2-D  Addition  Theorem:  If  B^x.y)  D  V-|(u,v)  and 
B2(x,y)  DV2(u,v),  then 

"B.  u.<fl  +  Ei  C*.y)  o  M  ,w.o-)  +Mz  m.*) 

r-  r°  .ia-irc ,  r  r,  .i«r<*un-'ry> 

=hh)J*  «■  V,cu.«0  +X2‘i/->r  e 

=•  B,  Cx.y)  + 

It  is  readily  apparent  that 

=  V,Cw,vO  +  V2  i 


2-D  Shift  Theorem:  If  B(x,y)^  V(u,v) ,  then 

t^TTCcuCt  +  bvr) 


Bc*.-a.,y-W)  =>  e 


Vdi.ir) 


and  \/cu..<l,0--\o)  Z)  e 


Proof: 


_,r  >  r  -t.z-rrCu.iL+'J'^ 

V  (tl-euy-lo) 

r  r  -iZ-rTCa*x+^)  .  Z 27T E^- ^  + Cvr.  U) y] 

=JLs* Cu^dEo-w)  e  <2.  \ 


V(U-<L,'^'*0 


Similarly,  we  have 


-£HS<n> 


^7TCa.u.  +  Vj^) 

=  e 


V  fa,ir) 


2~D  Similarity  Theorem:  If  B(x,y)  O  V(u,v),  then 

N/(a.u,W)  ^  UV»\  "6(  a:>ib") 

and 

»f  Proof:  1  f?  n  ; 

T  [ v<ccu,  v,*)}  c  N  ^sV) 


“  \Owl  |V»1 


4!  Vav^j  c  N  co-u, 


-  ■£&  ^U:>$) 
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Similarly,  we  haver-  r-  *£TT(ux+OrU) 

{bCo-x.W]  =JW^  e. 

=  6w>  e"**'***®,^ 


■&(*•*,  Vl) 


2-0  Definite  Integral  Theorem: 

j£*/«js  T5U,4)  =  \|(«)°) 


Pro  f:  tr  r.  ia.ir<UX+<»il) 


.Ml] 


iA.*a 
u*«  o 


[Wi]u-o  =  N/<o.o) 


Similarly, 


jQu.jQo-N/c^o'^  -  BCo.o) 


2-D  Differentiation  Theorem: 


OkT  (fy)  BU,^  =  (-lziruF(-iZirop  Me*.*) 


Proof: 


ftdB'aw  Buaj- =£*>  *  <w 


•J&  e^tfsTBu.,) 
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Integrate  by  parts  once  for  the  y  integral: 


The  integrated  part  vanishes  if 


Li*  B(x,y)  =  0 

a-** 


(as  we  will  be  concerned  with).  Perform  (m-1)  further  intergrations 
by  parts  for  the  y  integral,  and  then  n  integrations  by  parts  for 
the  x  integral;  the  result  is 


*(Cfcrt*T  (-w)  e  ^<*y) 


=  (,-UV uf”  (-  LZTTu-y*  T) 


Similarly, 

Bcx,y)  ■ 


Hermiticity  Theorem:  If  B(x,y)  is  a  real  function,  then 
V(u,v)  must  be  Hermitian,  i.e., 

V(u,v)  =  V*(-u,-v) 


Proof: 


.  .  VW  .  .  go  p*0 


4*7T(UX+  ^ 

e  BlM) 


Since  B  is  real,  we  have  C©5[2'n'(uX**'r^))  3CX.*Ji 

Vic  5,4  l2’7r^x+vr^]  tx,y) 

where  »V^(‘*.u'V -v  iVt-cn,^  . 


We 


clearly  see  that  NjyCU,^  -Vjv(.u  -u-V  ,  \/i(M.U')-~V<t-«.-4r) 


Hence  ^14,  U“)  -  V*(*M;”'r)  . 


Hermiticity  Corollary  1  --  If  B(x,y)  is  a  real  function 
and  is  even  in  its  joint  argument,  B(x,y)  =  B(-x,-y),  then  V(u,v) 

Is  real  and  even  in  its  joint  argument,  V(u,v)  =  Vr(u,v)  =  Vr(-u,-v). 

Proof:  From  the  equations  for  Vr(u,v)  and  V.(u,v)  in  the 
proof  of  the  Hermiticity  Theorem,  we  need  only  note  that  if  B  is 
even  In  x  and  y,  then  the  integrals  in  the  equation  vanish,  and 
the  Vr  equation  yields  Vr(u,v)  =  Vr(-u,-v). 

Hermiticity  Corollary  2  —  If  B(x,y)  Is  a  real  function 
and  odd  in  its  joint  argument,  B(x,y)  =  -B(-xt-y),  then  V(u,v)  is 
imaginary  and  odd  in  its  joint  argument. 

Proof:  From  the  proof  of  the  Hermiticity  Theorem,  this 
case  causes  the  Vf  equation  integrals  to  vanish  while  yielding 
Vj(u,v)  =~Vj(-u,-v) . 


From  these  theorems,  we  deduce  the  concepts  that 
reality  of  one  function  in  a  fourier  transform  pa1>  implies 
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Hermiti city  of  its  conjugate  function,  and  that  a  real  even 
function  has  a  real  even  conjugate  while  a  real  odd  function  has 
a  pure  imaginary  odd  conjugate. 

The  convolution  of  two  functions,  B^(x,y)  and  Bgfx.y), 
is  denoted  by  and  is  defined  by 

We  now  prove  a  very  useful  theorem  relating  the  Fourier  transform 
of  the  convolution  of  two  functions  to  a  product  of  their  separate 
transforms . 

2-D  Convolution  Theorem:  If  B-j  3  V]  and  B2  3  Vg.  then 


Proof: 


r 1 uu 1 •  w 


;Z7r<.ux+\rq) 


W)[M&  crm — 


=  M  6«’r0lR‘4’ri,^cc'J3 


•y- 


•  k'*Mi*r*'nuuum^tag 


53 


We  similarly  obtain  the  converse  of  this  theorem  as: 

V ,  iw.vr)  *  V»Cu ,<r)  ■=>  Bi  u.s)  •  Ba  «,>i)  , 


2.3.4  Truncation  Effects 


We  define  the  2-D  box  function  in  terms  of  1-D  rectangle 
functions,  following  Bracewell's  notation  (1978): 


a  rr  c«)  •  mu-) 


where 


l«l=t 

iw\>i 


Rectangular  truncation  of  a  visibility  function  In  the 
u-v  plane  is  equivalent  to  multiplying  that  visibility  function  by 
a  2-0  box  function  which  has  zero  amplitude  outside  the  specified 
rectangular  region.  Since  truncation  can  be  expressed  as  a 
product  of  functions  in  u-v  space,  then  the  corresponding  brightness 
function  can  be  determined  by  the  Convolution  Theorem.  To  consider 
the  general  case  of  rectangular  truncation,  we  must  know  the 
transform  of  -r  t»  /  li-C*  IT-  C.%*. 


mU- Cm  LT-  C.-  \ 

—  >—£7) 


which  specifies  a  box  function  with  unit  height  Inside  the  rectangular 


region  with  edge  lengths  and  by,  centered  on  cu  along  the  u  axis 
and  cv  along  the  v  axis,  respectively.  Since  this  transform  will 
be  important  later,  we  now  prove  a  theorem  for  the  general  case. 


Box  Function  Transform  Theorem- 


Proof: 


«*&•*&) 

-i*' rK(Cn^  -iMTXCCu- V)  ,**7T/ CCr*^  -WlPfCC^jg) 

_ *  ««  6>  &  —  €» 

-Aairx 

__  ^  ^iair^cuxtc>r*j)  j-#»4  ttWx  s/n 

irx  “  TpJ  ■ 

The  effect  on  the  brightness  function  resulting  from  the  truncation 
of  the  visibility  function  is  expressed  by  the  following  theorem. 


Truncation  Effect  Theorem:  If  V  B,  then  if  V 
represents  the  truncated  visibility  function, 

=.nk'fl«'  ^-i2-7r(c“xVc^') t,M7ra,/ 

try 

Proof:  For  2-D  rectangular  truncation  of  V(u,v),  we 


can  write 
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We  will  use  the  curly  overbar  'W  to  denote  a  truncated  visibility 
function  and  its  associated  transform.  Then,  by  the  Convolution 
Theorem 

-  T'[\7(w,vr)] 

-  e“4,air^CHX+t,r^ ajKirWX.  arwTWvr j 

TTX  -7T2J 

J—  «/-«■>  Tra' 

•  ©Ox-xMJ-W, 

Qualitatively,  truncation  of  V(u,v)  gives  rise  to  oscilla- 
tions,  or  ringing,  in  the  function  B(x,y)  as  compared  to  B(x,y). 

The  origin  of  the  ringing  is  seen  In  the  convolution  of  B(x,y) 
with  an  oscillatory  factor  e“'*2lT(cux+cvy)  ancj  With  the  sine 
functions,  which  have  decaying  oscillations.  The  effect  of 
truncation  can  also  be  qualitatively  described  in  terms  of  a  loss 
of  resolution  in  the  sharp  features  of  the  conjugate  function. 

This  is  especially  clear  for  the  conjugate  functions  of  Inter¬ 
ferometry,  since  truncating  V(u,v)  means  that  we  are  not  using  the 
high  resolution  information  for  large  u  and  v  outside  of  the 
truncation  limits,  and  therefore  B(x,y)  is  "blurred"  or  "smoothed" 
by  being  convoluted  with  the  sine  factors  to  produce  B(x,y). 

2.3.5  Representation  of  an  Averaged  Function 

ConHder  a  visibility  function  which  is  averaged  in 


rectangular  cells  over  the  u-v  plane.  Let  one  cell  be  centered  on 


the  point  (0,0).  Then  we  can  specify  any  cell  by  double  integer 
Indices  (n,m),  where  the  "n"  specifies  the  n^  cell  along  the  u 
axis  and  "m"  the  mth  cell  along  the  v  axis.  We  let  n  and  m  range 
from  -®  to  +«,  where  the  +/-  sign  will  denote  the  relevant  half¬ 
axis  which  the  cell  is  referenced  to.  The  (n,m)^  average  value 
of  V(u,v)  is  then  given  by 

-r. «,  , 

r  -  ,  i  ru**’5r  r**  a 

2  * 

where  Tu  n  and  Ty  m  denote  the  lengths  of  the  (n.m)^  averaging 
cell  along  the  u  and  v  axes,  respectively,  and  up  and  vm  are  the 
midpoint  coordinates  for  the  (n,m)th  cell.  The  overbar  and 
brackets  will  be  used  to  denote  an  average  value  of  the  visibility 
function. 

To  represent  the  average  value  of  V(u,v)  In  the  (n,m)th 
cell  as  a  function  of  u  and  v,  we  will  use  a  2-D  delta  function  at 
the  cell's  midpoint  whose  strength  is  given  by  the  average  value 
over  that  cell.  We  therefore  write 


Note  that  the  area  of  the  (n,m)th  cell,  (Tu  n'Ty  m) ,  has  been  included 
as  a  factor  multiplying  the  delta  function  in  order  that  the  delta 
function  properly  represent  the  average  value  over  the  enti re  cell. 

The  V(u,v)  function  averaged  over  its  entire  domain  can  now  be 
written  as  a  sum  over  these  average  value  terms,  and  will  be 
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designated  by  7*(u,v): 

Z"^Vm  [V-n.'-n,]  £  (H’My, 

A  few  consents  are  in  order  regarding  the  specification 
of  the  (n=0,m=0)  cell  as  being  centered  at  the  origin  (u=0,v=0), 
and  the  particular  choice  of  representation  for  the  average  value 
in  a  cell  as  a  function  of  u  and  v.  Specifying  that  the  (n=0,m=0) 
cell  be  centered  on  the  origin  is  based  upon  a  natural  generalization 
of  the  2-D  Sampling  Theorem  to  the  case  of  averaging.  The  Sampling 
Theorem  specifies  that  sampling  be  performed  in  a  regular  fashion 
with  one  sampling  point  located  at  the  origin.  The  effect  of  the 
epoch  of  sampling  will  be  further  discussed  below  in  our  considera¬ 
tions  of  the  Sampling  Theorem.  The  use  of  delta  functions  to 
represent  the  average  values  is  also  a  natural  generalization  from 
the  Sampling  Theorem.  The  other  reasonable  way  to  represent  the 
average  values  would  be  by  rectangular  boxes  whose  height  is  the 
relevant  average  value  In  the  particular  cell.  Now,  a  box  representa 
tion  for  a  particular  cell  can  be  expressed  as  the  convolution  of 
a  delta  function,  which  is  located  at  the  cell  midpoint  and  whose 
amplitude  is  the  average  value  of  V(u,v)  over  that  cell,  with  a  box 
of  unit  height  and  identical  cross  section  to  the  relevant  averaging 
cell;  . 

*7T(£;- 
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where  V°m  denotes  the  box  representation  function.  Consider  the 

case  of  constant  cell  sizes,  so  T  „  and  T  _  are  constants  for 

u,n  v,m 

all  n  and  m.  From  the  convolution  theorem,  we  immediately  see 
that,  if  66(x,y)  represents  the  transform  of  the  delta  function 
representation  for  the  averaged  visibility  function: 

then  B  (x,y),  which  represents  the  transform  of  the  box  representa- 

O 

tion  for  the  averaged  function,  can  be  written  as 

,  (  S1*' Sm‘ TT*73»^ 

-rz  tv  V  its  7r y 


Si"“  -rt-c  =  XX- 

we  see  that,  as  the  averaging  cell  sizes  go  to  zero,  B  and  B 

□  6 

become  identical.  However  for  large  Tu  and  Ty  ,  as  we  will  be 
concerned  with,  the  box  representation  results  in  preferentially 
decreasing  the  amplitude  of  &5  as  we  move  away  from  the  origin,  by 
the  factor 

(TU  sine  TiA*)(Tr  <  A.  poR- 


In  fact,  the  box  representation  produces  zero  brightness  at  all 
zeros  of  sin(irTux)  and  s1n(TrTyy).  We  would  like  Bfi  and  to  match 
the  true  brightness  function  B(x,y)  as  closely  as  possible.  However, 

A 

B  (x,y)  incorporates  a  sine  modulation  factor  which  must  have  zeros 

□ 


i 
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at  definite  points  in  the  x-y  plane,  regardless  of  the  amplitude 
of  B(x,y)  at  those  points.  This  consideration,  as  well  as  that 
of  the  natural  generalization  of  the  Sampling  Theorem,  makes  the 
delta  function  representation  the  preferrable  method  to  use  to 
represent  the  averaged  visibility  function.  A  comparison  of  these 
two  representations  for  a  number  of  one-dimensional  test  functions 
was  conducted,  and  those  results  support  the  choice  in  favor  of  the 
delta  function  representation. 


2.3.6  The  Standard  Fourier  Transform  for  an  Averaged  Function 


The  averaged  visibility  function  is  specified  by 

X  [  Vm.-w]  £(n,-Hv,,vr-i£^  . 

n*-*  m«-«* 

The  Fourier  transform  of  this  function,  designated  fi6(x,y)  will  be 
an  approximation  to  the  true  brightness  function  B(  ,y) .  Taking 
the  Fourier  transform  of  •  gives  ^  . 

\  ^  r-  -i  f r  r 

(  *“  *l)  ^  Tu,-n*73-,  ^ C  £ 

•mi-*  y 


If  V(u,v)  has  been  truncated,  then  we  will  have  only  a  finite  number 
of  terms  in  the  sum,  which  can  then  be  written  as 

A  P  -r-  -T-  rr,  l 


T*.,X.49,.,J  a" 


*w  ^ 


where  we  have  assumed  symmetric  truncation  about  the  origin  along 
the  u  and  v  axes. 

Now,  B(x,y)  must  be  a  real ,  non-negative  function  since 

it  represents  the  radiation  intensity  from  a  physical  body.  The 

visibility  is,  in  general,  a  complex  function.  We  can  therefore 

write  [v.  „]  =  [V„  m3  +i[V„  „].  ,  where  the  subscripts  refer  to 
n,m  n,m  r  n,m  i 

real  and  imaginary  parts.  It  appears  at  first  sight  that  may 
be  a  complex  function,  and  would  therefore  not  be  a  good  approxima¬ 
tion  to  B(x,y).  However,  we  have  shown  that  the  Fourier  transform 
of  a  real  function  must  satisfy  the  Hermiticity  property,  and  this 
property  is  invariant  under  averaging  operations.  We  will  see 
that,  for  an  untruncated  or  symmetrically  truncated  visibility 
function,  B,  will  always  be  real,  although  it  may  possibly  assume 
negative  values  contrary  to  the  nature  of  the  true  brightness 
function. 


Clearly,  since 


2  Z 

the  average  values  [Vn  m]  will  satisfy  a  "discrete  Hermiticity 


relation"  if  V  is  Hermitian: 


V Cu,vr\  *  V  (-  U ,-vr)  =5=t>[  V*n,Tn]  “  [  V->nr"*3 


This  statement  assumes  that  the  epoch  of  averaging  is  at  the  origin 
[(ns0,m=0Mu=0,v=0)].  If  we  write 
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then  for  symmetric  truncation  (or  no  truncation  if  N-*»,  we  have 

h  <mwE  £  x.„-c J [[ v,. J  +[v  j 

.  Sl«ZTT{*y>\+#^  r  sirf^Cu*  c«s2TT(a„  t+'Vjjj 

Since  [Vn  mlr  is  even  in  n  and  m  while  sin  ^(u^+v^)  Is  odd  in 
(un,vm),  then  their  product  is  odd  and  vanishes  in  the  double 
summation  process.  Similarly,  the  second  imaginary  term  vanishes 
in  the  double  summation.  Both  of  the  real  terms  are  even  in  (n,m) 
and  (un,v_),  so  we  can  write 

&£tX,iO=I  X.n,X;of[v^C'32^A+[^Ji3,N2^M 

Zl  t  f[ +[?„, J w<ai$w«4 

We  have  shown  that  §5  is  indeed  real  for  symmetric  truncation  of 
V(u,v)  (and  for  no  truncation  at  all).  Furthermore,  the  last 
equation  implies  that  a  knowledge  of  the  average  values  of  V(u,v) 
over  only  half  the  u-v  plane  is  required  to  specify  B^.  The  two 
quadrants  in  the  chosen  half  plane  must  be  adjacent,  not  diagonal 
to  one  another,  in  order  that  the  Hermiticity  property  be  used 
to  deduce  the  average  values  in  the  other  half  of  the  u-v  plane. 

If  the  visibility  function  Is  truncated  asymmetrically, 
then  Bfi  may  be  complex  valued.  The  Hermiticity  property  satisfied 
by  V(u,v)  for  real  B(x,y)  enables  us  to  relate  the  values  of  the 
visibility  function  in  two  adjacent  quadrants  of  the  u-v  plane  to 


its  values  in  the  other  two  quadrants.  Therefore,  in  cases  of 
asynsnetrlc  truncation.  It  may  be  possible  to  ascertain  the 
average  values  of  V(u,v)  needed  to  symmetrize  the  truncation.  If 
this  is  not  possible,  then  the  set  of  average  values  of  V(u,v) 
should  be  decreased  to  result  in  a  data  set  reflecting  symmetric 
truncation  in  order  to  produce  a  real-valued  B6(x,y).  One  dimen¬ 
sional  test  functions  were  examined  for  the  case  of  asymmetric 
truncation,  and  the  results  corroborate  the  assertion  that 

a 

asymmetric  trunction  produces  a  complex  B3(x,y)  with  the  imaginary 
portion  contributing  no  information  useful  in  discerning  B(x,y) 
or  a  best  approximation  to  it. 


2 . 4  The  Sampling  Theorem  for  Fourier  Transforms 

2.4.1  Statement  of  the  Sampling  Theorem  in  1-D 

The  Sampling,  or  Nyquist,  theorem  for  Fourier  transforms 
In  1-D  may  be  stated  as  follows.  Consider  a  function  B(x)  which 
is  non-zero  only  in  a  finite  portion  of  the  x-axis,  from  -Lx  to 
Lx.  Such  a  function  is  usually  referred  to  as  being  "band-limited". 
Let  V(u)  be  the  Fourier  transform  of  B(x),  and  let  V(u)  be  period¬ 
ically  sampled  with  the  epoch  of  sampling  at  the  origin  (i.e., 
a  sample  is  taken  at  u=0  and  periodically  thereafter  in  both 
directions).  If  the  sampling  interval  is  less  than  or  equal  to 
ir-  ,  then  it  is  possible  to  reconstruct  the  function  B(x)  exactly 

£S< 

[and  hence  of  course  also  V(u)],  The  largest  sampling  Interval 
which  can  be  used  and  still  reproduce  B(x)  undistorted  is  known  as 
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the  Nyquist  interval,  which  we  will  denote  by  T^. 

In  Its  one-dimensional  form,  the  Sampling  Theorem  is  of 
use  in  signal  processing  applications.  A  suitable  generalization 
to  2-D  makes  this  theorem  of  use  in  VLSI  data  analysis.  Since  the 
proof  of  the  Sampling  Theorem  in  2-D  is  entirely  analogous  to  the 
one-dimensional  proof,  we  will  prove  this  theorem  below  only  for 
the  2-D  case  of  interest. 

2.4.2  Relevance  to  VLBI  Data  Analysis 

Since  all  celestial  emitters  of  radio  waves  have  a 
finite  physical  size,  their  brightness  distribution  functions  have 
non-zero  amplitude  only  over  a  finite,  bounded  region  on  the 
plane  of  the  sky.  In  two  dimensions,  we  will  refer  to  such  a 
function  which  is  non-zero  only  in  a  certain  portion  of  the  plane 
of  definition  as  being  "region-limited". 

As  discussed  above,  an  interferometer  measures  the 
visibility  function  which  is  the  Fourier  transform  of  the  brightness 
function.  But  since  all  brightness  functions  are  region-limited, 
then  a  generalization  of  the  Sampling  Theorem  to  2-D  would  imply 
that  sampling  V(u,v)  in  some  periodic  lattice  fashion  over  the  u-v 
plane  would  provide  sufficient  information  to  reconstruct  B(x,y) 
exactly,  provided  that  the  sampling  cells  were  acceptably  small. 

The  ramification  of  this  theorem  for  radio  interferometric  observa¬ 
tions  using  aperture  synthesis  is  obvious,  and  with  this  motivation 
we  now  consider  the  2-D  Sampling  Theorem  in  di'ail. 
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2.4.3  Derivation  and  Discussion  of  2-D  Sampling  Theorem 


Consider  a  region-limited  brightness  function: 
|'B(x,y)  |x|<Lx,  |y|*Ly 

0  otherwise 


B(x,y)  * 


Let  e(x,y)  be  the  periodic  extension  of  B(x,y),  such  that 
B(x+pLx,  y+qLy)  »  B(x,y) ,  if  p,c 
8(x,y)  in  a  2-D  Fourier  series: 


8(x+pL  ,  y+qL  )  =  B(x,y),  if  p,q  are  integers..  We  can  expand 
x  y 


I  „  e 


Note that  2L  and  2L  give  the  full  source  extent  along  the  x  and  y 
x  y 

axes,  respectively.  The  Fourier  coefficients  are  given  by 

c  . 


21* 


fL*  jm*  f !  f *  iffi! £ 

Li*  «■  {rJ*  «•  ^  pa,*) 

r1*  f  w‘ i  .  /  a.tr-t\x  .  4Twy\ 

-  fiCx.MN  £  tC  ' 


£ 

' - M,*W1 

r1-*  rw>i 

■=  £ 

In  the  region  of  interest,  8(x,y)  B(x,y),  so  we  may  write 

£  I  e"iQ  >  **  M>K 


o 


Let  V(u,v)  be  the  visibility  function  which  is  the  Fourier  transform 

of  the  brightness  function  B(x,y): 

r  r°  .i.z-rr(M.x+vr^  . 

e  V(u,vr) 


Since  B(x,y)=0  outside  the  region  |x|.<LY,  |yjsLv,  we  have  for  the 

V  y 

in  erse  transformation: 
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&  B(x,y). 

Consider  the  particular  values  of  V(u,v)  at  the  regularly 
spaced  lattice  points 

/ ,  --yv  .  <-»*  \ 

(  '  {J’" 

Note  that  the  point  (0,0)  is  included  in  this  set.  At  these  points, 
we  have 


Comparing  this  equation  with  the  previous  one  for  C  „  shows  that 

n  ,m 


~  4L*Lvj 

Thus,  we  see  that  a  knowledge  of  periodically  sampled  values  of 
the  visibility  function  in  the  u-v  plane  provides  all  necessary 


Fourier  series  coefficients  C_  _  to  completely  determine  the 

n  9  »n 

function  B(x,y)  exactly.  Hence  we  can  write 


:  /  ZTT'n\  .  iJTPtnyA 

Biyjj.r  r .  c,.„  c  *  -scf) 


I  -  '.C  WW  - aC 


I  ft  a.  .  hv.  s  + 

-4tT;Z  Z  v(-2n--2L7)e.  ZL‘  **■> 

Of  course,  exact  knowledge  of  B(x,y)  implies  that  V(u,v)  is  also 


completely  specified  from  the  Fourier  transform  relations. 


We  therefore  conclude  that  for  a  region-limited  brightness 
function,  knowledge  of  periodically  sampled  values  of  the  associated 
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visibility  function  completely  specifies  B(x,y)  and  V(u,v).  The 
proper  lattice  constants  for  the  periodic  sampling  lattice  are 

AH  *  -  X,  *x  f  Atr  *  S.  TXx 

where  Tu  ^  and  Ty  ^  refer  to  the  "Nyquist  interval"  lengths 
along  the  u  and  v  axes,  repectively. 


Because  B  is  a  real  function  and  V  is  therefore 
Hermitian,  we  can  further  reduce  the  sampling  requirement.  Only 
two  adjacent  quadrants  need  be  periodically  sampled  In  order  to 
reproduce  B(x,y)  and  V(u,v)  exactly,  since  samples  of  V(u,v)  in 
these  two  quadrants  are  related  to  the  sample  values  (on  a  regular 
lattice)  In  the  other  two  quadrants  by  the  Hermiticlty  property. 

If  we  write 


and  expand  the  complex  exponential  in  terms  of  sines  and  cosines, 
we  can  simplify  to  get 


H  tit,'  °)«» +>k<3ti  °) 


Note  the  similarity  of  this  result  to  that  obtained  for  Bg(x,y); 

If  the  averaging  intervals  T  n  and  T  were  fixed  at  the  Nyquist 
intervals,  and  we  replaced  the  average  values  by  sample  values, 
then  the  results  would  be  identical  (since  un=^-  and  vm=^-  now). 


./ 
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2.4.4  Aliasing 

Conceptually,  the  aliasing  effect  in  1-D  is  identical 
to  the  effect  in  2-0.  Therefore,  we  will  base  our  discussion  and 
work  in  one  dimension,  and  state  the  very  obvious  generalizations 
to  two  dimensions. 

In  one  dimension,  the  Sampling  Theorem  may  be  written 
for  a  real  function  B(x)  which  is  band-limited  in  |x|sLx,  as: 

"Tm*  N/r  (®)  4*  £  "TniE  j^v  X 

+  VifVl«)  %,H  ZTTyC^x.  Xj 

•j  "* 

where  j  =  ,  and  2Lx  is  the  full  1-D  source  extent.  The 

expression  written  for  B(x)  Is  recognized  as  just  a  standard 
Fourier  series  expansion  for  that  function  in  Its  domain  of 
definition  [-Lx,l_x].  Note  that  the  periods,  Xp  Nj,  of  both 

trig  functions  are  defined  by 

Y  -»  — 1 -  -  2Lk_ 

Af,  MX 

The  fundamental  period  occurs  for  n=l ,  with  all  other  values  of  n 
producing  shorter  harmonic  periods.  Since  all  terms  In  the  sum 
have  the  common  beat  period  21  ,  we  see  that  the  superposition  of 
the  various  harmonic  terms  produces  a  function  6(x)  which  Is 
identical  to  B(x)  in  the  interval  [-LX,LX]  and  reproduces  B(x) 
periodically  along  the  entire  x  axis  with  period  Xp_^  j^j=2Lx  In  the 
bands  [rl  ,(r+2)Lv],  r=odd  Integer.  This  is,  of  course,  the 
standard  periodic  extension  of  a  function  when  represented  by  its 
Fourier  series.  What  is  important  to  recognize  for  our  considera- 
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tions  is  that  the  particular  choice  of  TNI  for  the  sampling 
Interval  produces  a  Fourier  expansion  which  properly  "fits"  Into  the 
true  domain  of  definition  for  B(x),  [-LX,LX].  e(x)  reproduces  B(x) 
in  each  band  without  interference  between  bands  if  we  sample  at 
intervals  T^j  because  all  terms  In  e(x)  have  the  coitmon  beat 
period  2LX>  Figure  2.5a  illustrates  this  case. 


Suppose  now  that  we  consider  sampling  V(u)  at  intervals 


'NI 


result 

then  is 

:  f\u 

ml 

hfTH)  '-•**hr^r\T<  x 

+  s/H 

1 

-L. 

v  JL 

*r> 

>  -TSx- 

The  period  is,  now 

^f,<  * 

and  the  common  beat  period  is 


B<(x)  is  the  periodic  extension  for  this  Fourier  series  which 
matches  B(x)  in  the  interval  , 

l-k'ix) 


and  reproduces 

with  r  an  odd 
period  here  is 
between  bands. 
B<(x)  which  Is 
this  interval. 


B(x)  an  all  intervals 

T  JC, 

L  ’  %T<  * 

integer.  We  Immediately  see  that,  since  the 
greater  than  Xp_-|  NI  we  still  prevent  Interference 
Note  that  we  may  define  a  band-limited  function 
Identical  with  B(x)  In  [-LX,LX]  and  zero  outside 
but  whose  basic  domain  of  definition  is  chosen  as 
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[-L  -6  ,  L  +6  ],6  >0.  If  y  (u)  B  (x)  then  the  sampling  theorem 

XX  A  A  A  ^  ^ 

properly  applies  here  with  T<,Ni=2fir+$T'  But  clearly  V(u)sV<(u) 
since  the  Fourier  transform  integral  is  performed  over  the  entire 
one-dimensional  space.  As  illustrated  In  Figure  2.5b,  B<(x)  is 
identical  with  B(x)  in  [-Lx »LX3  but  its  periodic  bands  are  now 
spaced  a  distance  apart  instead  of  being  juxtaposed.  Since  we  are 
really  interested  only  in  the  function  B(x),  we  conclude  that  sampling 
at  T  <Tni  Is  capable  of  reproducing  B(x)  exactly  in  the  region  of 
Interest. 


Now  consider  T^T^.  as  the  sampling  interval.  The  component 
periods  in  8>(x)  now  are 

^r,>  “  *vr>  < 

and  the  common  beat  period  for  all  terms  is 

Ty  <  TjJ  * 

We  have  a  situation  where  other  bands  have  a  non-harmonic  overlap 
into  the  interval  [-LX,LX],  and  B(x)  is  not  exactly  reproduced  in 
its  basic  domain  of  definition.  Thus,  sampling  at  T>TNj  causes 
distortion  of  B(x).  We  will  presently  show  by  an  argument  using 
the  Convolution  Tneorem  [based  on  a  similar  discussion  by 
Bracewell  (1978)]  that  the  effect  of  using  T>  is  to  simply  cause 
^(x)  to  be  a  superposition  of  overlapping  bands  wherein  e^x) 
would  match  B(x)  except  for  the  overlap  effects.  It  Is  not 
possible  to  deduce  the  true  form  of  B(x)  by  compensating  for  the 
overlap  in  some  way,  and  thus  sampling  V(u)  at  ™NI  results  In  an 
unrecoverable  loss  of  Information  on  B(x). 


•V 


We  can  concisely  prove  these  observations  by  using  the 


Convolution  Theorem.  We  define  the  Sampling  Function  by 

3ZE(t)2T-J L.5U-*t) 

where  T  is  the  sampling  interval.  Sampling  V(u)  at  intervals  T 
is  equivalent  to  multiplying  V(u)  by  III(ij!-).  Denoting  the  sampled 
visibility  function  by  Vs(u),  we  write  Vs(u)»V(u) -III (^) -  Let  b(x) 
be  the  Fourier  transform  of  Vs(u).  Then  the  Convolution  Theorem 
gives  us 

(3U)=  T'{ N/Siu^  =  *  T  ‘[mute)] 

=  bu)  *  #*'  [30:  (t)}. 

The  Fourier  transform  of  a  sampling  function  is  another  sampling 
function,  III(u)  III(x)  [see  Bracewell  (1978)  for  a  discussion 
of  the  proof  of  this  transform  pair].  By  the  Similarity  Theorem, 
we  then  have 

jf ‘'{^Or)]  =  f  -rzrL-r*). 

5° '  =  Bex)  *  -r2U-r%) 

.rv Bw'iif 

niw 

****-«• 


This  result  clearly  reveals  the  effect  of  the  size  of  the 
sampling  inteval .  For  B(x)  band-limited  in  [-LX,LX],  then  clearly 
TNI=2T~  causes  e(x)  to  be  a  non-overlapping  periodic  extension  of 
B(x).  Similarly,  for  T<T^j,  we  simply  widen  the  bands,  as  was 
illustrated  in  Figure  2.5b.  For  T>T^j,  then  =p<2Lx  and  we  see  that 
the  bands  now  overlap,  producing  a  distored  version  of  B(x)  in  the 
Interval  of  Interest  [-LX,LX].  The  present  statement  of  this  result 
lucidly  shows  that  the  distortion  results  from  superposition  of 
different  bands,  each  of  which  may  be  conceived  as  containing  one 
undistorted  period  of  B(x).  {See  Figure  2.5c) 

The  term  "aliasing"  has  been  used  to  describe  this  effect 
because  sampling  at  Intervals  greater  than  the  Nyquist  causes 
spatial  frequencies  In  different  bands  to  mix  due  to  the  overlapping 
of  the  bands.  It  is  therefore  said  that  higher  spatial  frequencies 
are  posing  as  lower  frequencies,  and  hence  the  higher  frequencies 
have  lower  frequency  aliases  because  of  the  overlapping. 

The  generalization  to  two  dimensions  Is  quite  straight¬ 
forward.  Sampling  In  2-0  is  performed  on  a  rectangular  lattice  in 
the  u-v  plane,  and  sampling  along  any  line  of  lattice  points  in  the 
u  or  v  directions  is  independent  of  behavior  in  the  other  direction 
and  hence  exactly  equivalent  to  the  1-D  case.  Me  can  therefore 
think  of  the  2-D  lattice  as  being  a  direct  product  of  two  1-D 
lattices.  We  then  have 

B(x.'i)  (tl*, To. 4) 


where 


M  (-r^T.'i)  s  UT(ruK)^ILCXr‘i) 

•llJz  s(r^2 

The  result  is: 

We  immediately  see  that  all  of  the  previously  discussed 
cases  for  one-dimension  apply  to  each  argument  of  B  independently. 

In  other  words,  for  true  recovery  of  B(x,y)  in  the  region  of 

interest  {[-lx,Lx],[-Ly,Ly]} ,  we  must  have  Tu  <  Tu  NI  =  and 

1  ’x 

Tv  <  Tv  |^j  =  .  If  the  sampling  interval  along  u  or  v  exceeds 

Its  relevant  Nyquist  length,  then  aliasing  effects  occur  and 

B(x,y)  will  be  distorted  along  the  respective  direction  (i.e., 

if  TU<TU  but  TV>TV  NJ,  distortion  of  B(x,y)  occurs  only 

parallel  to  the  y  axis). 

We  will  call  a  rectangular  cell  with  edge  lengths  Ty  ^ 
and  Ty  ^  a  Nyquist  cell.  When  we  say  that  a  given  cell  is  smaller 
than  the  Nyquist  cell  size,  we  mean  that  T  <T„  KIT  and  T  <T  MT 
for  the  given  cell.  A  cell  which  is  referred  to  as  being  larger 
than  a  Nyquist  cell  will  have  TU>TU  NI  or  Tv>Ty  Nj,  or  both. 

2.4.5  Effect  of  the  Epoch  of  Sampling 

Heretofore,  we  have  used  the  Sampling  Theorem  with  the 
epoch  of  sampling  at  the  origin,  which  Is  the  conventional  statement 
of  the  theorem.  We  now  investigate  the  effect  of  an  epoch  of  sampling 


i 

.1 


not  based  on  the  origin.  Again,  we  will  work  in  1-D  and  state  the 
generalization  to  2-D. 


Suppose  that  we  sample  V(u)  at  regular  intervals  T,  but 
we  do  not  center  the  samples  on  u=0;  this  Is  equivalent  to  shifting 
the  Sampling  Function  by  some  Increment  a,  with  0<a<T.  Then 

ur($)  ($-*-). 

The  sampled  visibility  function  V$(u)  is  now 

V,(*)=  vm- nr  (>-«•). 

From  the  Shift  Theorem,  we  have 

The  Convolution  Theorem  then  requires  that  the  Fourier  transform 

°f  Vs<u)  fHcu^BGO *  e  “~t  arcr*) 

-  -rfr 

~n»—c  ' 

As  before,  B(x-^)  is  a  periodic  extension  of  the  band-limited 
function  B(x).  Now,  however,  3  (x)  has  a  complex  factor  in  its 

terms.  Note  that  y<l ,  and  if  is  irrational,  then  e”i27rrY 
Is  non-periodic.  If  is  rational,  then  e’^2irn^  is  periodic 
in  n.  Suppose  that  n  Is  the  Integer  period  such  that  na=T  and 
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for  any  integer  r.  Then  the  complex  function  6ft(x)  will  have  a. 

period  ^  In  particular,  for  T=TNj,  Ba(x)  has  a  period 
2TU 

n2Lx  =  — —  >  2LX.  Hence,  If  T<T^,  no  overlapping  occurs,  and  we 
see  that  for  n=0,  we  do  In  fact  recover  B(x)  in  the  Interval  of 
interest  [-LX,LX]. 

Consider  the  case  T>T^j  now.  Clearly,  we  will  have  over 
lapping  as  discussed  previously,  but  the  distortion  now  will  be 
far  worse  than  when  the  epoch  of  sampling  was  at  the  origin, 
because  the  adjacent  bands  which  mix  into  the  region  of  interest 
have  complex  weighting  factors  in  them. 


The  generalization  to  2-D  incorporates  a  translational 
displacement  of  the  epoch  of  sampling  for  the  lattice  from 
(u=0,v=0)  to  (a,b).  The  result  is 


f  ©(x-T?: .  £*iZir^e 

*  •  *»)»■'•  W|>«0  ' 


-il  7T-£ 


The  cases  discussed  above  for  1-D  apply  separately  now  to  each 
coordinate. 

As  a  simple  illustration  of  the  effect  of  the  epoch  of 
sampling  In  one-dimension,  let  us  consider  a  case  where  a=^-  for 
the  conjugate  pair  V(u)  =  sine  u  l-|x|=B(x)  for  |x|<l.  We  have 
TNi  *  a  *5.  The  complex  modulating  factor  is  now  e"1lrn,  with 
period  n=2.  But 
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n  even 
n  odd 


The  resulting  function  6a(x)  for  T=TNi  is  pictured  as: 


as  compared  with  the  standard  Sampling  Theorem  result  for  sampling 


Based  on  the  considerations  of  this  section,  it  is 
considered  preferable  to  use  an  epoch  of  sampling  at  the  origin. 

In  our  discussions,  we  will  want  to  consider  averaging  cells  with 
dimensions  both  smaller  and  larger  than  the  Nyquist  cell.  When  we 
represent  the  average  value  of  V(u,v)  over  a  cell  by  a  delta 
function  at  the  cell  midpoint,  we  are  imitating  the  sampling  process, 
and  hence  the  observations  on  the  effects  of  cell  size  versus  Nyquist 
cell  size  are  applicable.  If  we  average  over  the  entire  visibility 
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function  with  cell  sizes  equal  to  or  smaller  than  the  Nyquist 
cell  size,  then  aliasing  effects  will  not  be  present  and  the 
transform  of  the  averaged  function  will  represent  B(x,y)  except 
for  possible  distortion  due  to  the  fact  that  average  values  were 
used  instead  of  true  sample  values.  Averaging  with  cell  sizes 
greater  than  the  Nyquist  cell  size  produces  greater  distortion  due 
to  aliasing  -^fects.  As  we  have  seen,  an  epoch  at  the  origin 
produces  less  distortion  than  a  shifted  epoch  in  cases  where  the 
Nyquist  size  is  exceeded,  and  we  would  therefore  prefer  to  use 
an  averaging  epoch  at  the  origin  (u=0,v=0)  for  these  cases.  Since 
we  will  be  considering  various  averaging  cell  sizes,  and  since  the 
epoch  does  not  affect  the  result  in  the  region  of  interest  for 
cells  smaller  than  the  Nyquist  cell  whereas  It  may  adversely  affect 
the  result  for  cell  sizes  greater  than  the  Nyquist  cell  size,  we 
see  that  fixing  the  epoch  of  averaging  at  the  origin  for  all  cases 
of  cell  size  (including  mixed  cell  sizes)  is  the  best  prescription. 
It.  should  be  noted  that  in  an  actual  experiment,  it  may  not  be 
possible  to  choose  an  epoch  of  sampling  at  the  origin.  However, 
although  the  elliptical  tracks  may  not  conform  with  a  sampling  or 
averaging  epoch  at  the  origin,  if  the  data  are  gridded  onto  a 
rectangular  lattice,  then  it  would  be  preferable  to  use  a  lattice 
with  epoch  at  the  origin  if  such  an  extrapolation  is  feasible. 
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CHAPTER  3 


THE  STANDARD  FOURIER  TRANSFORM  METHOD 


3 . 1  Statement  of  the  Method 

In  an  actual  experiment,  V(u,v)  can  be  examined  over  only 
a  finite  portion  of  the  u-v  plane,  and  hence  the  visibility  function 
is  effectively  truncated.  We  will  assume  that  the  truncation  is 
symmetric  in  u  and  v  (or  has  been  made  symmetric  as  previously 
discussed).  We  have  shown  above  that  the  Fourier  transform  of 
the  truncated  and  averaged  visibility  function  is 

'BJx.'i)  =1  *iru.x 

HX.M  Vn»| 

+[9^Ji  »«  2Tr6ux+*.y)j 

where  T  ,T  are  the  averaging  cell  lengths  in  the  u  and  v 
u  ,n  v  »ni 

directions,  respectively;  \,vm)  is  the  midpoint  of  the  (n,m)th 
cell  in  which  the  average  value  is 


For  completeness,  we  state  the  result  for  asymmetrical 


truncation: 


For  asymmetrical  truncation,  £,(x,y)  will  be  complex. 


We  can  easily  show  that,  for  a  given  truncation  of  the 
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visibility  function,  the  value  of  Bg(0,0)  is  independent  of  the 
particular  averaging  lattice  which  is  used,  and  is  in  fact  equal 
to  B (0,0),  the  origin  value  of  the  function  resulting  from  trunca¬ 
tion  of  V(u,v)  alone,  without  averaging.  We  have  shown  that  the 
effect  of  truncating  V(u,v)  produces 

t 


where:  c..  =  i  coordinate  value  of  the  midpoint  of  the  non-truncated 
portion  of  V(u,v) 

b.  =  length  along  the  i  axis  of  the  non-truncated  part  of 
V(u,v) . 

By  the  Convolution  Theorem, 

TTl ^)] 


Since  the  integrations  are  independent  of  x  and  y,  we  have,  using 
the  definition  of  rJ  »  ^^5')  j 


We  can  write  £gfx,y)  in  a  form  which  applies  to  either  symnetric  or 
asymmetric  truncation, 


fegbAyt  f>TwRvJ 

i-M, 


t 
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Thus,  §6(0,0)  =  B(0,0)  =  Volume  under  the  visibility  function  inside  the 
truncation  limits. 


3.2  Illustrative  Case:  The  Two  Component  Gaussian  Source 

We  now  want  to  consider  a  practical  example  to  examine 
the  effects  on  resolution  due  to  truncating  and  averaging  the 
visibility  function.  As  our  model  source,  we  chose  a  two  component 
Gaussian  brightness  distribution  function.  This  model  source  has 
been  chosen  not  only  because  of  its  relevance  to  actual  observed 
sources  which  hav°  two  or  more  peak  components  in  the  brightness 
function  which  we  wish  to  resolve,  but  also  because  it  provides 
an  excellent  case  to  study  changes  in  symmetries,  the  vector 
between  the  peak  brightness  points,  and  the  width  of  the  components 
as  functions  of  truncation  limits  and  averaging  cell  sizes. 

The  basic  form  of  the  brightness  function  for  this  model 


source  is 


_Hu£U-oW] 

BlX,1!)*  e  +  e. 


where  w  Is  a  parameter  characterizing  the  width  of  each  component. 
The  Gaussian  components  have  peak  brightness  of  unity  at  the  points 
(-1,0)  and  (1,0).  Figure  3.1  illustrates  the  cross-section  of  the 
model  source  for  four  different  values  of  w.  Figure  3.2,  a-d, 
presents  three  dimensional  views  of  the  model  source  for  the  four 
values  of  w,  as  well  as  their  associated  visibility  functions. 

(Note  that  these  source  and  visibility  functions  are  smooth 
functions;  any  sharp  edges  in  the  3-D  plots  are  due  to  the  mesh 
size  used  in  the  plotting  algorithm.) 


3.2.1  Exact  Brightness  Distribution  function  and  Its  Transform 


The  Fourier  transform  of  the  brightness  function 


_ .  .irw[(x-i)z+alJ  -irufU+4  +*J  J 

+e 


is  derived  as  follows.  Let 


Note  that 


Vu,vr)=  t  [bim)} 

-irw[xV|*J 


By  the  Similarity  Theorem  and  the  Separable  Product  Theorem,  we 
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immediately  have 

*{ 


u3  6. 


and  by  the  Shift  and  Addition  Theorems  we  have 

~f  { s(K-o  +  £U+o}  =  z  c. •*  . 


Hence,  by  the  Convolution  Theorem,  we  get 


SlCuy)^  2*71*4.  & 


Note  the  ease  with  which  this  transform  is  obtained  using  the  theorems 
as  opposed  to  a  direct  calculation. 


Since  the  double-Gausslan  model  that  we  are  using  is  not 
a  truly  region-limited  function,  one  cannot  properly  define  a 
Nyquist  cell  size.  However,  the  Gaussians  fall -off  sufficiently 
fast  so  that  we  may  define  effective  limits  for  the  extent  of  B(x,y) 
in  the  x  and  y  directions.  This  effective  length  will  be  specified 
as  twice  the  value  xc  for  which  B(x,y)  has  an  amplitude  less  than 
or  of  the  order  of  .001  of  its  peak  amplitude  for  all  x>xc,  and 
similarly  for  yc>  The  effective  Nyquist  interval  lengths  are 
then 

1  1 
Tu,NI  =  and  Tv,NI  s  Tq  • 

Table  3.1  shows  the  effective  cut-off  lengths  and  corresponding 
Nyquist  intervals,  as  well  as  the  peak  values  of  B(x,y),  the 


location  of  these  peaks  in  the  brightness  function  and  the  half¬ 
widths  [defined  as  half  the  peak  value  in  B(x,y)]  for  the  four 
width  values  w  used  in  this  study. 

It  was  desired  to  examine  cases  of  both  severe  and 
negligible  truncation  of  the  visibility  function.  Since  V(u,v)  is 
sharply  peaked  and  falls-off  very  fast  for  w=.25,  whereas  it  decays 
quite  slowly  for  w=10,  It  was  decided  that  a  common  truncation 
limit  of  |u|=l  and  | v | =1  for  all  cases  studied  would  achieve  the 
goal  of  examing  both  radical  and  minimal  truncations. 

Three  different  sets  of  averaging  intervals  commensurate 
with  these  truncation  limits  were  used.  One  set,  designated  A, 
used  averaging  cells  which  were  all  smaller  than  the  effective 
Nyquist  cells  for  the  w=10  and  w=l  cases  (but  were  larger  than 
the  Nyquist  cells  for  the  other  two  w  cases).  Set  B  used  cell 
sizes  which  were  all  greater  than  the  effective  Nyquist  cell 
sizes  for  all  four  w  cases.  Set  C  used  cell  sizes  which  were  a 
mixture  of  cells  both  smaller  than  the  smallest  effective  Nyquist 
cell  and  larger  than  the  largest  effective  Nyquest  cell  in  the 
four  w  cases.  All  cells  were  taken  to  be  souares,  i.e.,  the  u 
and  v  edge  lengths  of  any  given  cell  were  the  same.  The  three 
different  averaging  schemes  are  summarized  below,  where  It  Is 
understood  that  the  Intervals  along  the  v  axis  are  the  same  as 
the  Intervals  along  the  u  axis  which  are  illustrated. 


Averaging  Scheme  A  i  •*  -Z  •  <*.  ■»* .  •*.  ,tf.  ♦**.  •£  .  j 


Averaging  Scheme  B  i  .4>25*  .  ,  .3ff  ,  At  ,  .42S* 


Averaging  Scheme  C 


.5“  .  .iff  .2  .1  .2  .ZS 


3.2.2  Method  of  Analysis 


All  of  the  program  references  in  this  section  refer  to 


Program  A  of  the  Appendix. 


r 

t 


The  function  B(x,y)  was  obtained  by  performing  the 
relevant  convolution  integral  numerically.  Since  the  magnitude  of 
B (x,y)  is  less  than  .001  for  all  x  and/or  y  values  greater  than 
3.5  for  the  three  cases  w=10,  1,  .5  (the  corresponding  value  along 
the  x  axis  for  the  case  w=.25  is  3.96),  the  highest  x  and  y  values 
for  which  B(x,y)  can  be  caluclated  must  be  3.5  less  than  the 
magnitude  of  the  x'  and  y'  limits  on  the  respective  integrals.  In 
other  words,  since  we  must  put  a  practical  bound  on  how  far  out  we 
go  to  Integrate  over  x'  and  y'  in  the  convolution  integral,  we 
must  Insure  that  the  chosen  limits  are  sufficiently  far  out  so 
that  effectively  all  of  the  volume  under  the  shifted  function 
8(x-x',  y-y ' )  In  the  intergrand  is  included  in  the  numerical 
integration.  The  falling  sine  factors,  which  are  less  than  unity 
except  at  the  origin,  actually  help  to  taper  the  integrand, 
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sinc(tu,nx)  sinc(Ty  j^y)  B(x-x ' ,y-y;) ,  so  that  the  cut-offs  at  3.5 
are  somewhat  conservative. 

The  expense  and  CPU  time  involved  in  performing  the 
convolution  integrals  to  a  reasonable  level  of  accuracy  (about 
.05  of  true  values)  limited  the  range  of  x  and  y  values  from  0 
to  2.5  in  increments  of  .25.  Furthermore,  B(x,y)  was  actually 
evaluated  only  at  lattice  points  in  one  quadrant.  This  was 
considered  acceptable  because  the  integrand  arguments  are  even  in 
x  and  y,  and  truncation  of  V(u,v)  was  also  performed  symmetrically, 
so  that  the  values  of  B(x,y)  in  one  quadrant  should  be  symmetrically 
related  to  its  values  in  the  other  three  quadrants. 

The  effective  limits  used  for  both  the  x'  and  y'  integra¬ 
tions  were  -6  to  6,  in  accordance  with  the  comments  above.  A  mesh 
with  step  size  .1  in  both  x'  and  y'  was  used  to  numerically 
evaluate  the  integrals.  The  simplest  integration  technique  of 
summing  over  cell  volumes  (value  of  integrand  at  cell  midpoint 
times  the  cell's  base  area)  was  used  for  speed;  improved  accuracy 
would  involve  at  least  a  fourfold  increase  in  computer  time.  To 
generate  one  value  of  B(x,y)  required  14, 400  iterations  of  the 
sunmation  loop  to  evaluate  the  integrals  numerically,,  so  that 
even  the  very  sparce  x-y  lattice  containing  100  points  required 
considerable  computer  time.  Attempts  to  Increase  the  density  or 
extent  of  the  lattice  points  examined  were  considered  cost- 
ineffective  for  our  purposes  here. 


Figure  3-3,  a-d,  illustrates  the  visibility  function  after 


truncation.  Truncation  is  negligible  for  the  cases  w=.5,  .25,  and 
therefore  the  long  convolution  integral  calculations  were  omitted 
for  those  cases. 


To  obtain  an  overall  comparison  of  B(x,y)  with  E(x,y), 
the  following  parameters  were  defined: 


Normalized  Deviation  = 


Normalized  Root  Deviation  Squared 


A _ 

[BIX,*#- £1*,'$*’ 


where  the  integrals  are  performed  over  a  square  region  of  the  x-y 
plane  whose  area  is  A.  The  division  by  area  was  included  to 
normalize  the  result  so  that  comparisons  could  be  made  with  other 
cases  where  the  section  of  the  x-y  plane  which  was  examined  had  a 
different  area.  Since  we  evaluated  B(x,y)  only  over  the  region  0 
to  2.5  along  the  x  and  y  axes,  the  integrals  above  had  to  be 
numerically  performed  over  this  single  quadrant.  The  mesh  step  size 
was  .25  along  the  x  and  y  axes  (the  available  data  values  for  5(x,y) 
have  this  increment  step),  which  is  quite  large.  Therefore,  the 
calculated  parameters  should  not  be  viewed  as  very  accurate 
quantities  here,  but  rather  as  quantities  which  are  only  sufficiently 
accurate  for  comparitive  purposes  with  other  similarly  calculated 
parameters . 
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Other  tests  (such  as  peak  locations  and  halfwidths)  were 
not  performed  for  B(x,y)  due  to  the  prohibitively  long  calculations 
needed  at  each  step  of  the  iterative  procedures  required  for  such 
tests. 


Numerically  performing  the  averaging  process  for  the 
visibility  function  was  considerably  faster  and  more  accurate, 
using  a  mesh  with  step  .025  along  both  the  u  and  v  axes.  The 
visibility  function  was  numerically  integrated  over  an  averaging 
cell  and  the  average  value  calculated  by  dividing  by  the  averaging 
cell's  area.  Only  the  upper  half-plane  was  used  (by  virtue  of 
the  Hermiticity  relation)  and  the  values  of  §5(x,y)  were  calculated 
using  the  standard  discrete  Fourier  inversion  method  discussed 
in  Chapter  2. 

Several  tests  were  then  performed  to  study  the  function 

A  A 

B  (x,y)  and  its  relation  to  B(x,y)  and  B _(x,y).  The  two  main  peaks 

y  o 

A 

in  B.(x,y)  were  located  using  an  iterative  search  procedure 

o 

incorporating  the  subrotine  VA06A  of  the  Harwel  Library.6  The 
accuracy  control  for  this  calculation  was  specified  by  the  condition 

H  b\*=  p_to)2+  [2Api]*<.oi 

for  the  iteration  to  terminate. 

The  magnitude  and  angle  (with  respect  to  the  x  axis)  of  the 
vector  connecting  the  two  peaks  was  then  calculated.  Since  the 
half-widths  of  the  separate  Gaussian  components  are  not  well- 
defined  in  the  x  direction  due  to  overlapping  of  the  two  components. 


\ 


I 
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the  determination  of  the  half-width  for  each  component  In  B.(x,y) 

0 

was  performed  in  the  y  direction  where  half-width  is  well- 
defined  (i.e.,  unaffected  by  overlap  effects).  The  half-width 
was  defined  as  the  distance  between  the  y  value  at  a  component 
peak  and  the  y  value  where  the  brightness  magnitude  is  half  the 
peak  magnitude.  The  half-widths  were  determined  by  iterative 
search,  and  the  error  bound  criterion  specified  that  the  magnitude 
of  Bfi(x,y)  be  within  .01  of  the  true  half-width  magnitude  for 
the  iterative  search  to  terminate. 


Overall  comparative  parameters  were  defined  for  Bfi(x,y) 
in  an  analogous  fashion  to  those  used  for  B(x,y).  For  comparison 
with  the  exact  brightness  function,  the  parameters  were: 


Normalized  Deviation  = 


_ 


VM-is 

Normalized  Root  Deviation  Squared  =  - 2 - 


The  square  region  (with  area  64)  from  -4  to  4  along  both  the  x  and 
y  axes  was  examined,  and  the  mesh  step  size  used  to  perform  the 
Integrations  numerically  was  .2  along  both  the  x  and  y  axes.  Again, 
the  results  here  should  be  viewed  as  sufficiently  accurate  for 
comparative  purposes  but  not  highly  accurate  in  themselves. 

For  overall  comparisons  with  B(x,y),  the  following  parameters 


were  used: 
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LL  [Su.'D-'B.Cx.'i)] 

Normalized  Deviation  =  - - ^ - 

Normalized  Root  Deviation  Squared  = - ^ - - - 

the  same  first  quadrant  region  and  mesh  size  as  used  for  the 
comparison  of  El(x,y)  with  B(x,y)  was  used  here. 

B.(x,y)  was  also  examined  to  see  if  the  reflection 
symmetries  through  the  x  and  y  axes  which  are  present  In  B(x,y)  are 
preserved  in  B5(x,y).  The  region  from  -6  to  6  along  both  the  x 
and  y  axes  was  examined  using  a  mesh  with  step  size  .1  along  each 
direction.  The  criterion  for  symmetry  between  two  reflection 
symmetric  points  was  that  the  two  values  of  §5(x,y)  agree  with 
each  other  to  within  .001,  and  the  condition  for  an  overall  judge¬ 
ment  that  B  (x,y)  was  symmetric  under  reflection  through  the  x  and 
y  axis  was  that  no  more  than  10  pair  of  points,  of  the  7200  pairs 
examined  for  each  reflection  synmetry  case,  be  unsymmetric.  This 
tolerance  of  about  .12  was  incorporated  to  allow  for  any  possible 
rounding  errors  or  other  inaccuracies  in  the  computer's  calculations 
which  may  have  been  present  and  incorrectly  imply  an  asymmetry. 


\  .y 
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3.2.3  Results 

Table  3.2  summarizes  the  results  of  the  model  study. 
Every  case  was  found  to  preserve  the  reflection  symmetries  through 
x  and  y  axes,  and  therefore  these  results  were  not  included  in 
Table  3.2  Similarly,  the  component  peaks  were  symmetric  about  the 
origin  In  all  cases  studied,  so  that  only  the  second  component 
peak  (at  positive  x  value)  was  included  in  the  Table. 
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3.2.4  Observations  and  Discussion 


Preservation  of  the  reflection  symmetries  of  B(x,y)  in 

B.(x,y)  is  an  expected  result,  since  V(u,v)  was  symmetrically 

<5 

truncated  and  the  averaging  cells  were  symmetric  in  the  sense  that 


T  =T  and  T  =T  _ . 
u,n  u,-n  v,m  v,-m 


The  component  peaks  in  B  (x,y)  are  generally  not  coinci- 

0 

dent  with  the  peaks  in  B(x,y).  In  some  cases,  the  peaks  are  seen 
to  shift  inward  along  the  x  axis  toward  the  origin,  in  other  cases 
the  shift  is  outward;  however,  in  all  cases,  there  is  no  shift  in 
peak  location  along  the  y  direction.  Peak  shift  is  most  severe  for 
averaging  scheme  B,  where  the  cells  are  all  larger  than  the 
relevant  Nyquist  cell  size  (for  all  four  w  cases).  In  an  extensive 
study  of  about  300  one-dimensional  sample  cases,  it  was  found  that 
peak  shifts  result  from  truncation  effects  alone  as  well  as  from 
the  combined  effects  of  truncation  and  averaging.  Qualitatively, 
we  can  account  for  peak  shifting  by  realizing  that  a  peak  in  B(x,y) 
results  from  the  overall  cumulative  reinforcement  of  the  component 
sinusoids  of  which  B(x,y)  is  composed  in  its  Fourier  expansion,  and 
therefore  omitting  some  of  these  sinusoids  [due  to  truncating  V(u,v)], 
altering  their  relative  amplitudes  [due  to  averaging  V(u,v)], 
and  using  a  limited  set  of  "non-Nyquist"  expansion  "frequencies" 
(specified  by  the  aveaglng  cell  midpoints)  can  give  rise  to  a  rein¬ 
forcement  peak  at  a  point  which  Is  displaced  from  the  peak  In 

B(x,y).  Three  cases  in  the  model  survey  show  a  peak  In  B.(x,y) 

0 

coincident  with  the  peak  in  B(x,y),  to  within  the  accuracy  limits 
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of  the  calculation  as  discussed  above.  Two  of  these  cases  are  for 
averaging  scheme  A,  where  the  cells  are  smaller  than  the  respective 
Nyquist  cell  size.  The  truncation  of  V(u,v)  was  not  severe  for 
these  cases.  A  very  small  shift  (i.e.,  undetectable  to  our 
accuracy  limit)  is  therefore  not  surprising  for  these  cases.  The 
third  case  involved  radical  truncation  of  the  visibility  function 
and  used  averaging  scheme  C  which  incorporates  mixed  cell  sizes; 
apparently,  the  peak  shifting  Influences  of  truncating  and 
averaging  V(u,v)  tended  to  cancel  In  this  case  producing  a  very 
small  (i.e.,  undetectable)  net  peak  shift. 

The  peak  magnitudes  of  B(x,y)  are  seen  to  differ  from 

6 

the  peak  magnitudes  of  B(x,y).  In  most  cases,  the  peak  magnitude 

of  B  (x,y)  was  less  than  the  peak  magnitude  of  the  corresponding 
6 

B(x,y);  however,  in  two  cases  the  converse  is  seen  to  hold.  For 
w=10,  all  three  averaging  schemes  show  peaks  in  B  (x,y)  which  are 
significantly  less  than  the  peak  In  B(x,y).  The  visibility  function 
was  severely  truncated  in  this  case,  and  therefore  we  would  expect 
that  the  omitted  Fourier  components  had  significant  amplitudes, 
so  that  a  peak  in  B(x,y)  results  from  a  greater  number  of  reinforcing 
Fourier  components  with  significant  amplitude  than  are  available  to 
form  a  peak  in  Bs(x,y).  (The  relation  B.(0,0)=iT(0,0)  for  any 

o 

averaging  scheme  holds  here.  Note,  however,  that  different  averaging 
schemes  with  the  same  truncation  limits  can  lead  to  different  peak 
brightness  magnitudes  here  because  the  peaks  are  not  at  the  origin.) 
The  effect  of  truncation  is  only  one  factor  accounting  for  the  differ- 
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ence  in  peak  magnitudes  between  B„(x,y)  and  B(x,y).  The  relative 
amplitudes  and  "frequencies"  of  the  component  sinusoids  in  B5(x,y) 
are  determined  by  the  particular  averaging  scheme  which  is  used, 
and  these  factors  also  influence  the  peak  brightness.  For  w=10, 
the  effect  of  truncation  is  dominant,  whereas  for  the  other  three 
cases,  the  peak  magnitudes  are  determined  by  a  more  balanced 
contribution  of  the  above  cited  factors.  In  particular,  the  two 
cases  where  the  peak  magnitude  of  B5(x,y)  exceeds  the  peak  magnitude 
of  B(x,y)  can  be  qualitatively  viewed  as  cases  in  which  the  combined 
factors  result  in  greater  reinforcement  magnitude  for  than  for  B. 

The  halfwidth  of  B5(x,y)  is  similarly  influenced  by 

truncation  and  averaging  effects.  From  the  previously  discussed 

theorem  for  the  effect  of  truncation,  we  know  that  truncating 

V(u,v)  leads  to  a  broader  and  less  sharply  resolved  function  B(x,y) 

as  compared  with  the  true  source  brightness  function  B(x,y).  For 

the  case  w=10  where  truncation  was  severe,  we  clearly  see  that  the 

truncation  effect  is  dominant  and  the  halfwidths  are  all  larger 

(by  about  a  factor  of  2)  than  the  halfwidth  of  B(x,y).  In  the 

other  cases  where  the  truncation  and  averaging  influences  are  more 

balanced,  the  differences  between  the  halfwidths  of  B.(x,y)  and 

0 

B(x,y)  is  less  severe  than  the  case  w=10.  (These  differences  are 
about  10%  as  compared  with  the  factor  of  2  difference  for  w=10.) 

The  parameters  defined  to  enable  overall  comparisons  to 
be  made  between  the  various  cases  confirm  the  expected  results  that, 
generally,  truncation  without  averaging  yields  a  better  approximation 
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to  the  brightness  function  than  does  truncation  with  averaging,  and 
the  smaller  the  averaging  cells  the  better  the  overall  approximation 
of  B6(x,y)  to  6(x,y).  Again,  these  are  general  trends  and  it  is 
possible  (as  occurred  in  several  cases)  for  B5(x,y)  to  give  a 
better  approximation  to  B(x,y)  than  does  B(x,y)  for  some 
particular  choi  es  of  truncation  limits  and  averaging  scheme.  This 
observation  is  supported  by  sample  cases  among  about  300  one- 
dimensional  test  cases  where  B(x)  and  B^(x)  were  plotted  and 
visually  inspected  for  goodness  of  fit  to  B(x).  The  entries  in 
Table  3.2  clearly  show  that  averaging  scheme  B  (cell  sizes  all 
greater  than  the  Nyquist  cell  size)  was  always  significantly  worse 
in  fitting  B(x,y)  than  were  the  other  two  averaging  schemes  where 
some  or  all  of  the  cells  were  smaller  than  the  Nyquist  cell  size. 

Having  obtained  a  qualitative  understanding  from  the 
model  study  of  the  effects  on  the  brightness  function  resulting  from 
truncating  and  averaging  the  visibility  function,  we  now  analyze 
these  effects  quantitatively  in  Chapter  4  to  further  elucidate 
the  nature  of  these  effects. 
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CHAPTER  4 

THEORY  OF  TRUNCATION  AND  AVERAGING  EFFECTS 

Having  seen  the  effects  of  truncation  and  averaging 
illustrated  in  Chapter  III,  we  now  take  an  analytical  look  at  the 
causes  of  these  effects. 

4.1  Periodicity  Considerations 

Our  discussion  on  the  periodic  nature  of  the  brightness 
function  resulting  from  sampling  or  averaging  the  visibility 
function  has  thus  far  been  concerned  with  a  fixed  sampling  or 
averaging  interval.  We  now  discuss  the  periodicity  phenomenon  for 
irregular  intervals,  i.e.,  intervals  of  varying  sizes.  We  will 
work  in  one  dimension  and  then  generalize  to  2-D. 

Theorem  4.1:  Consider  a  visibility  function  V(u)  and 
let  it  be  averaged,  with  averaging  intervals  in  the  set  {Tn>.  Then 
B5(x)  is  a  periodic  function  in  x,  with  period 

x„  -  LCM[|l-|]n,0 

where  LCM  =  Least  Common  Multiple. 

Proof:  For  the  general  case  of  truncation  and  averaging 
(including  asymmetric  truncation)  we  have 


--.r  M  YT’W..El_ 


■arum* 


®*O0=-I  T.  Cv,] 


e 


=  f  T,  [[t*,],  C03  3/tf  2tTU^ 

*n-  -W« 

-*•  [.KX  S'w^TTM-  -  [v^  «•»  £-r«„x]j 


For  a  given  n,  all  trig  functions  have  the  common  period  Xp  n  = 
The  number  of  distinct  periods  is  just  (Ng+N-j+l).  The  period  of 
B  will  therefore  be  the  lowest  beat  period  of  these  "component" 


peri ods : 

X,»UCM[XtJ  * 


We  specify  n/0  because  uQ=0,  for  which  the  trig  functions  reduce  to 
constants.  Therefore,  the  n=0  term  is  automatically  commensurate 
with  the  common  period  of  all  the  other  components.  Now,  un  is  the 
midpoint  of  the  nth  averaging  interval,  and  if  umi-n  is  the  lower 
truncation  limit,  we  can  write 

+  ~T~  . 

If  we  use  a  constant  averaging  interval  T  and  truncate 

symmetrically,  then 

'~V\  €. 

a  ’Hv.t,  +ZX  +•  ■%- 
*  H*,  +  (■ti-/+w+i)t  +  ^ 

=  +  (.  +-*  +  +1T 


and 


106 


Then, 


=  /.<w[) 


11 


But, 


+(  a  ^  +  2^4-07-  -J^^o  • 

'  *0  —  -  (a.M+i)  t  ;  Sc 

Xp-  U-CN\  [ )  -r  '  ^  +-  2  M  f  ;u»  +»  |]  ,  =  LCA^[|-^-  i|] 


-m*o 


This  agrees  with  our  previous  results  for  fixed  T,  including  the 
limiting  case  where  V  is  not  truncated  so  that  N-*»  (provided  that 
for  then  jjpO) . 


In  two  dimensions,  we  have 

^6  ^  x>  ^  "  Z  Z 

Clearly,  if  we  fix  one  variable,  say  y,  then  B.(x,yJ  Is  a  periodic 

0  0 

function  in  x  with  a  period 

*n^o 

exactly  analogous  to  the  1-D  case  considered  above.  Thus,  both 
orthogonal  components  of  the  two  dimensional  brightness  function 
satisfy  the  1-D  periodicity  relation  independently. 


r-  -;*7TG4„X-»-OUJi) 


In  a  real  experiment  many  different  cell  sizes  occur, 

and  their  midpoint  coordinates  generally  do  not  form  a  set  whose 

inverses  have  a  low  Least  Common  Multiple.  For  example,  if  two  values 

of  |^-|  in  the  set  are  1  and  .617,  then  their  LCM  is  617,  without 
un 

even  considering  other  midpoint  values.  This  observation  leads  to 


jK- 
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the  conclusion  that  the  x  and  y  periods  of  B  (x,y)  will  be  very  large 
for  a  realistic  case,  and  hence  aliasing  effects  should  not  be 
important  (unless  the  data  is  gridded  onto  a  regular  lattice  with 
cell  sizes  greater  than  the  Nyqulst  cell  size). 

4.2  Effect  of  Truncation  Alone 

We  have  previously  shown  that  truncation  of  the  visibility 
function  V(u,v)  into  a  rectangular  region  with  sides  bu  and  by,  and 
center  at  (cu,cv),  results  in  a  brightness  function  B(x,y)  given  by 

^  U  )  TTbs*  5C  3/*i  TTUt*! 

where  B(x,y)  V(u,v).  B(  x,y)  exhibits  ringing,  i.e.,  it  resembles 
B(x,y)  (to  an  extent  depending  upon  the  severity  of  truncation),  but 
with  truncation- induced  oscillations  as  a  modulating  factor,  and  B 
is  smoother  than  B(x,Y). 

4.3  Effect  of  Averaging  Only 

In  the  next  several  sections,  we  investigate  the  effects 

of  averaging  over  the  entire  visibility  function  and  averaging  a 
truncated  visibility  function.  Initially,  we  will  work  In  one  dimen¬ 
sion  with  a  fixed  averaging  interval  size.  We  then  generalize  to 
the  case  of  irregular  averaging  intervals,  and  finally  we  will 
generalize  the  results  to  two  dimensions. 
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The  following  theorem  will  be  considered  our  basic  theorem 
relating  an  average  value  of  the  visibility  function  to  the  true 
source  brightness  function.  We  will  later  formulate  an  algorithm 
which  uses  the  information  on  the  effect  of  averaging  conveyed  by 
this  theorem  in  an  attempt  to  improve  the  brightness  function 
resolution  by  accounting,  to  some  extent,  for  the  fact  that  V(u,v) 
has  been  averaged.  The  standard  Fourier  inversion  method  does  not 
incorporate  this  information. 

Theorem  4.2:  For  the  case  of  averaging  with  the  epoch 
of  averaging  at  the  origin,  each  average  value  of  the  function  V(u) 
can  be  written  in  terms  of  a  definite  integral  over  x  which  Includes 
B(x),  the  true  conjugate  function  for  V(u),  in  the  integrand. 

Proof:  Let  V(u)  be  averaged  with  Integer  index  ne[-®,»] 
denoting  the  nth  averaging  Interval  with  epoch  (n=0)  at  the  origin, 
such  that  the  sign  of  n  denotes  the  relevant  half-axis  on  which  a 
particular  interval  is  located.  The  ntfl  average  value  of  V(u)  is 
given  by 

[  V-J  *  W-u,.  T>  ^  ' 

We  can  replace  V(u)  in  this  expression  by  its  Fourier  transform 
relation:  -j- 

m  l  /_>  B<*)] . 

z 


We  assume  that  the  order  of  integrations  can  be  reversed  (which  is 
true  for  analytic  functions,  and  we  will  justify  this  even  for  the 
case  where  B(x)  is  a  generalized  function,  say  representing  a  point 


source).  Then, 

[<}  -  tjj 

_  ^  T 


■t&«  iwr? 


*  (  I***-*  e,-™  ) 

.  x  fdx  b(k)  ei2,n4'x 

7TX. 

=  jfd*  bu)  =«.t;x  e*2’™'* 

=  [ ftlx  BCx)  *"■«■"£ X  c.anu-xl  =  f  f B(.xl  s->Xx) 

«'■"*  vuuM  *■  Jv. 

Note  that  if  we  average  with  a  fixed  interval  T,  then  un=nT.  We  also 

note  that  if  B(x)  is  a  band-limited  function  in  the  interval  [-LX,LX], 
then  the  definite  integral's  limits  need  be  only  from  -Lx  to  Lx. 

We  must  now  prove  several  intermediate  results  which 
will  be  needed  in  our  development. 


Theorem  4.3: 


r-»  - iZir^nTy,  i  *  s'  ,  w\  \ 

I  <=■  =7 I  SCx-^J 


Proof: 


rE  S(x.-^)=  m  (tx) 


'f  {3mCTx]}  =  ^3r(f)=Z. 


Now, 


no 


8ut'  **’{*«*)}-£ .HI ea7(, M„ 

s  _  -iZirvTX 

*■£  C 


Hence, 


*2  «*-*)  =2  t 


-i2  n'+T*. 


~n  «-« 


Qualitatively,  we  see  that  this  result  is  reasonable  because  if  x=y, 
then 


<£ 


so  the  sum  is  infinite.  However,  for  x^j,  the  infinite  sum  over 
sines  and  cosines  with  varying  phases  will  vanish.  ■ 


Corollary  4.3A: 

-tv*  ^.iarn-r*  ,lftZ7rMT  (x-^) 

where  N  and  Blare  integers. 


Proof: 

w+*\ 


X  £(u-  ‘n T)  - 


T 


3zr^r)-'ir(- 


u-tt-T 

2  NT 


)■ 


Then  analogous  to  the  previous  proof. 


w+d 
fin l 


-IZ’W'ti "Tx 
C 


But,  by  the  Convolution  Theorem, 


Ill 


*l+n &)  •  tt^)}=  m-rx)  *[*"*”"*&&] 

*  +£  SU-X'-^r) 

*  w  X  T|r-€ 

-i27r<J0T-(*- 


55  -i- 


rl  e 


s**Ztt  fjr(  x-^-) 


Comparing  results,  we  find 


'FI*  e - 

-H=-^+Jrt 


-  ?  P  -ix.Tr'tnr*.  ^^irfH^Cx-  ¥) 

ttU- 


Corollary  4.3B: 

^-iZirnTx  a,w^TA/T(x-^-) 

w,*— »  "^C  X  -  ^) 

Proof:  Setting^l=0  in  Corollary  A,  this  result  is 

obvious . 


M 


We  can  now  apply  these  intermediate  results  to  prove 
several  theorems  on  averaging  effects. 

Theorem  4.4:  The  effect  of  averaging  V(u)  over  its  entire 
domain  with  constant  averaging  interval  T  is  to  form  the  periodic 
extension  of  B(x)  multiplied  by  sinc(Tx). 

Proof:  As  previously  defined, 

'r[yJ 


I 


11 2 


Using  Theorem  4.2  for  [Vn],  we  write 

fiM- rf  l/jx'bcv)  «T«'  e ",n,r,‘j  {(*.„„) 


"Bs  (X)  = 


&■  eT'^rS  jJTjx'  boo  wr*e~ 


T£  e 


J.tO 


For  averaging  with  the  epoch  at  the  origin,  u  =nT,  so  we  can  write 


~  -  rr - - -  -  ’  z>-'  *  -p  - - -  ■  "  ■  ’ 

&(.*)  =  t£  Mx'  B(x')  *^TX'  e"*  1W) 

•  -nc>«  “'■‘® 

__  ,  £  -t2T^T(X-x’) 

~  T  /dx  B 60  3>*c I  e 

=  /V  B(X')  swc. 7x'f  <ffrx-x')->] 


where  we  used  Theorem  4.3.  So, 


Btu)=  £ 


*■**  -«p 


This  is  a  periodic  extension  of  the  true  brightness  function  modulated 
by  the  sine  factor.  Note  the  implications  for  aliasing  if  B(x)  is 

band-limited  in  [-LV,LV]  but  Tv^i— ,  or  if  B  were  not  band-limited 

XX  b 

at  all . 
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We  now  consider  the  important  case  in  which  the  visibility 
function  has  been  symmetrically  truncated  and  averaged  with  epoch 
at  the  origin. 


Theorem  4.5:  The  effect  of  averaging  a  function  V(u) 
which  has  been  symmetrically  truncated  is  expressed  by  a  convolution 
integral  which  contains  B(x),  the  true  conjugate  function  to  V(u), 
in  its  integrand. 


we  have 


Proof:  If  V(u)  is  symmetrically  truncated  and  averaged, 

r<u.)=  TZ  [vj  £(u.-v,t) 


where  un=nT  for  epoch  at  the  origin.  Proceeding  exactly  as  we  did 
In  the  previous  theorem,  we  get 

Bs(x)=  Tjjj*'  B(x') 

Since  the  summation  is  no  longer  from  -<*>,+«,,  as  It  was  In  the 
previous  theorem,  we  cannot  simplify  the  convolution  Integral  any 


further.  Noting  that 

U  & 


•  \  +  2  I  c  OS  [2TT-r)T (  X-’X')]  f 

*nti 


we  can  give  an  alternate  statement  of  the  result: 

=  T JJ*  t3(x')  "7 *'  [l  +  ZZ  cos ZtTrtT C X-X ')] 

The  presence  of  the  cosine  factor  accounts  for  the  ringing  expected 
when  the  visibility  function  is  truncated.  A  t.lr'rd  alternate  form 
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can  be  written  using  Corollary  4.3B: 

Aoo-f  fTx'  Be*')  *~7*' 

4.4  Examination  of  the  Theorems 

We  now  examine  the  theorems  in  further  detail.  We  will 
check  for  consistency  with  known  results  for  the  limiting  case  as 
T-+0,  and  then  discu  s  the  case  where  V(u)  Is  a  linear  function. 

This  case  is  important  because,  if  V(u)  is  linear,  then  averaging 
and  sampling  produce  identical  results,  and  thus  the  distortions 
predicted  by  our  theorems  should  vanish  for  t=tni- 

4.4.1  Limiting  Cases 

If  V(u)  is  averaged  over  its  entire  domain,  we  found 

BSU)*L  BCx-'t') 

-yt»  -•® 


The  limiting  case  is  then 

Lin A  "Br(*) 

T->0  6 


-rr-(x-  r, 


"T-*0  *"»»*■  "**e 


Using  the  theorem  for  the  limit  of  a  product  (Kaplan,  1952),  which 
states 


bw,  [  <UV  ^00]  =  [j.«*  <JW)]( 


we  can  write 


BAA 

•p+o  ® 


Z  [u 
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Now, 


£■* 

T-*  o 


"Bco 

BC*— ) 


IF  <-K*  O 

IF  o 


where  we  understand  B(x-“)  to  mean  the  function  B(x)  displaced  to  », 
and  hence  aliasing  effects  are  no  longer  relevant  (assuming 

J-i*  BCX)  =  © 


or  that  B(x)  is  in  fact  a  band-limited  function).  Also, 

5//r  7T7X  TTkc  -  7fTx  St*7Tr\ 

T-+P  !T~r CX-—J  T"+o  ^  T'TLk-'^) 

5W  TTX_  f  © 

1  1  iF“*.*o 


i*,  ~  -  t,«  f- 

TT(K-~)  T~o  i 


} 


—  LiM  C^vSTti 

T-*o  TT  A-TTfS. 


Hence,  we  have 

LiM  [b]  -  BU) 

-TV*© 


where 


'BCX') 

BOC-tf) 


and  the  Knonecker  delta  is 


IF  S*.=  O 
IF  St*0 


I 


/ 

I 


We  have  obtained  B(x)  as  expected  for  T-K5,  in  which  case 
we  are  effectively  sampling  every  point  in  V(u).  Note  that  B(x), 
not  its  periodic  extension,  is  obtained. 


If  we  symmetrically  truncate  V (u)  and  average  with  intervals 
T  with  epoch  at  the  origin,  we  found  that  -  ~ 

ri  -lZTT^T  (x-x') 


=  tJj*'  "-cT*Xe 


Then,  .  *  ,  ,  -r-.r  I  -r-  -lzTT-aTC*'*  J 

Lit*  BfCx)  •*  Lm  T  «ix  B(,x. )  sttte  7x  2.  ^ 

-pv«  k  r-*»*  <J~*c 

*  fix'  Bu1)  u.  [r  **.  -rv£ 

T*ou  w»s-N  J  ‘ 

Note  that,  in  this  case,  we  require  N-^°  as  T-+0,  in  order  to  keep 
the  truncation  limits* (N+^)T  finite  and  fixed  in  value  as  we  take 

the  limit.  Consider  the  limiting  term: 

■  I  ■  /n  'iZTPflf  (X“X  ) 

L  s  Lm  T s/ncTx  £  e 

T-*0  Wls-M 

r«  -T-  1 1  r.  __ «  -Ai27r*7)T”C x-x )  i 

as  Um  S/KC  rx  J  |  L»a»  TT  c,  J 

lT-*o  J  Lr-*e  ^5?-^  J 


s»  [l].[u  T*£  6 

L  J  L  tv* 


-t^TS iT(X-X’) 


Um  5  a*  2TrNT(K-K,-'*r) 
JZ2>  **  —  7T(X-X'-  !£.) 

A..  * 


NT  Wire 

where  we  have  used  Corollary  4.38.  So, 

,  a  ,  jw  aTl/rQu-x1-^)  « 

L _ir-' 

NT  Ritcrt 


S/M  *T&*T)<X-xV  ....  . 

- TTU-^ -  ,r*-° 


5/H  g7T(A/T)CX.-X-a0) 
TCx-x'-oo) 


|FM.*C> 


For  n?<0,  the  term  is  understood  to  represent  the  sine  function  shifted 
to  «,  and  since  the  magnitude  of  the  sine  function  goes  to  zero  as 
(x  -x')-»»,  we  see  that  in  any  region  of  Interest  for  a  band-limited 
brightness  function,  (x-x')<®,  and  then 

si*  2 

i r(x-x'-o©) 
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if  effectively  zero.  Our  result  then  is 

Lm  -RCxi./j*'  Bex')  f^jCaA/rKx-x'J 

T-vo  5  TT4X.-X'j 

But,  for  symmetric  truncation, 

^=&  gcxo 

where  bu  is  the  length  of  the  non-truncated  portion  of  V(u), 

bu=umax”umin=^N+^"!'=^  ^or  w^e  NT  remains  finite.  Hence, 

Aim  B  ( x) =B ( x ) , 

T+0 

as  expected.  It  is  straightforward  to  show  that  this  result  is 
also  obtained  if  V(u)  were  asymmetrically  truncated,  but  we  will  not 
go  through  the  details  here  since  the  asymmetric  truncation  case 
is  not  of  interest  for  our  purposes. 

4 .4£  Linear  Visibility  Function 

If  V(u)  is  a  linear  function  with  a  real  fourier  trans¬ 
form,  then  V(u)=iautb,  and  B(x)  must  then  be  a  generalized  function: 

?(.%)=  T '‘{v(w)}  —  §£  %  Sm  +  \.s«l 

where  we  have  used  the  derivative  theorem. 

We  first  show  that  the  previously  discussed  reversal 
of  Integration  orders  is  valid  if  B(x)  is  a  generalized  function. 
Suppose  B(x)=d(x)  and  V(ul=l.  Then 
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If  we  reverse  integration  orders,  we  get: 


"kJj*  ^u)[X,.T^du  £‘Z7rux]' -luLi* 


=Jl*  £«)  =- 1 


Thus,  reversing  integration  orders  is  acceptable  if  B(x)=$(x). 
Suppose  now  that  V(u)=i2T*und  B(x)=^[<$(x)] .  Then, 


iZiru*. 


&Cx)J 


a  -  i-ZlT 


^[(U.+  ^f-Cu.-l^)2']  *  ~^[2lL« Ty,] 
-  -iziru^- 


If  we  reverse  integration  orders,  we  get 

,x,+'E 


e*  J 

ii*  e.[£u)] 


=  i. 


‘■'a’rM”x  S«T7XX 


_  _  «i  r  _  1 

cU  Le  svwcT^xJ^ 

*  -  [ea7riu^L[^i;x|^ 


=  -  LZTTU^ 


I 


Thus,  if  V(u)  is  an  arbitrary  linear  function  such  that  B ( x )  is 
impulsive,  then  reversal  of  integration  orders,  which  was  essential 
in  deriving  the  theorems,  is  valid. 


Consider  now  the  case  where  V( u ) =1  so  that  B(x)=5(x). 

The  theorem  for  averaging  over  the  entire  visibility  function  with 
fixed  interval  T  then  says 

BrU)=t  BU-*) 

=f  7r"n>XV^)- 

TT'T 

To  analyze  this  result,  consider  the  following  integral.  Let  f(x) 
be  some  function  of  x. 

I j&[ v 

=£ 


J£*' J  T  <*)  • 


Therefore,  the  generalized  function 

t  £(-/.-»  — 

behaves  as  T 

and  so  we  may  write 

(*•)=£ 
d  *ni-* 

which  is  the  periodic  extension  of  B(x). 


I 
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Suppose  now  that  V(u)=i2Tru  so  that  B(x)=^{6(x )].  Then 


we  have 


_  V  «*t*-*)  <L  _ Tr?  v  *34 S'/At  7TT(x-^.) 

-ir-'-ir-  .«*-•*) LbL*  t;j  tttlx-'?) 


To  analyze  this  result,  consider  the  following  behavior,  where  f(x) 
is  some  function  of  x: 


=  -E  [f  <*'♦»&■  5^V 


* 


*w»-«o 


=  -£  m) 

-A[£4rUc*-*)l  4£*)] . 

HZ-.0 


We  therefore  see  that  the  generalized  function 


'~h= -*> 


behaves  as 


S/Af  ~7T  T  (  x.-  go 

-n-T(x-^) 


2jfe 

Hence,  we  see  that 


b*=2  it€<x-^-j 

■is***  ' 

the  periodic  extension  of  the  brightness  function  P(x). 
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We  conclude  that,  if  we  average  a  visibility  function 
V(u)=iau+b  with  Fourier  transform  B(x),  then  we  get  for  the  Fourier 
transform  of  the  averaged  function  a  periodic  extension  of  B(x) 
with  period  y.  Note  the  similarity  with  the  Sampling  Theorem  result 
for  point-wise  sampling  at  intervals  T,  which  is  the  expected  result 
here  because,  for  a  linear  V(u), 

1  f(n+?)T 

Average  Value  =  [V  1  =  r  ,  V(u)du=  V(nT)=V(u_)=Sample  Value 
n  1  j(n4)T  n 

Note  also  that  If  T>T^j  we  will  get  aliasing  in  B^  since  we  have 
effectively  point-wise  sampled  V  at  intervals  T>TNj. 

We  have  tested  the  theorems  on  averaging  for  limiting 
cases  as  T*0  and  for  the  case  of  a  linear  visibility  function  with 
constant-interval  averaging  over  its  entire  domain.  All  tests  have 
indicated  that  these  theorems  are  consistent  with  the  general 
properties  of  Fourier  transforms.  We  will  not  examine  the  case  of 
a  symmetrically  truncated  and  averaged  visibility  function  here  due 
to  the  algebraic  complexity  involved  in  that  case.  However,  it  is 
reasonable  to  expect  that,  based  on  the  other  successful  consistency 
tests,  this  case  also  will  yield  the  required  result  that  averaging 
a  linear  visibility  function  is  equivalent  to  sampling  that  function. 

4.5  Peak  Shift  Effects 

One  of  the  most  important  pieces  of  Information  to  be 
obtained  from  interferometric  observations  is  the  location  of  a  point 
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source  (or  source  below  the  resolution  limits)  or  the  peak 
brightness  for  an  extended  source.  As  we  saw  in  Chapter  3,  it  is 
possible  for  truncation  and  averaging  effects  to  distort  the  bright¬ 
ness  function  so  as  to  cause  shifts  in  the  location  of  the  peak 
brightness  away  from  the  true  source  peaks.  This  case  can  easily  be 
demonstrated  for  the  case  of  point  sources. 


Consider  two  point  sources  in  one  dimension  with  bright¬ 
ness  function  _  ,  r  „  .  .T 

^  L  +  SCx.+i)] 


The  visibility  function  is  then  V(u)=cosiru.  Suppose  we  truncate 
V( u)  at  ] u | =1  and  average  with  intervals  of  T=.4.  Figure  4.1 
Illustrates  the  true  source  brightness  function  and  compares  it  with 
B(x)  and  Bs(x).  The  peak  shifts  are  quite  evident. 


In  this  section,  we  attempt  to  discern  the  origin  of 

peak  shifts  and  conditions  for  a  peak  in  B  (x)  to  be  coincident 

6 

with  the  corresponding  peak  in  B(x).  We  will  discuss  truncation 
and  averaging  effects  separately,  and  then  their  combined  effect 
will  be  examined.  Although  we  are  most  interested  in  the  peaks,  the 
discussion  will  actually  deal  with  relative  extremum  points.  For 
the  case  of  averaging  over  the  entire  visibility  function,  we  can 

express  B,(x)  directly  In  terms  of  B(x):  however,  for  B(x)  and 

0 

for  the  case  of  averaging  a  truncated  visibility  function,  the  expres¬ 
sions  Involving  B(x)  are  convolution  integrals  and  will  be  difficult 
to  work  with.  For  this  reason,  we  will  analyze  the  case  of  averaging 
over  all  of  V  in  terms  of  B(x),  whereas  we  will  analyze  the  other  two 
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cases  in  terms  of  the  visibility  function. 


^B(x)  as 


The  Derivative  and  Convolution  Theorems  let  us  write 


h  £(.*)  [-UTU  v<«)  7T(-S^)} 

-LZttux 


-L 


c'-% 


dn.  LZiru.  €. 


\J(  u) 


We  will  set  Cu=0  so  that  the  truncation  is  symmetric.  Let 
V(u)=Vr(u)  +  1V.j (u)-  we  ^can  then  write 

cos  & ruyt,  -uVr  2-jtu.X  'j 


—  L  (u  Yr  +  uVi  3</J/?7TUx)J 


The  imaginary  terms  are  odd  functions  of  u,  so  they  vanish  in  the 
symmetric  integration.  If  ^B(xQ)  [3^-  B(x)lXi,x  =0>  then 


cl 


_  ^  u 

,)  =  (uVt  2.7TUX.—  uVV  ^  2-rrux}j 


Similarly,  we  can  obtain 


«  feu 

~~i  SUol  = -4ir2 JZciu  fu2\l,rni>sZirux„ 

+i £Vi*«2iru>C). 

It  Is  evident  that,  in  general,  ^6(xo)=0  does  not  imply  g^6(xo)=0. 

Although  ^  B(xQ)  may  be  small.  It  will  consistently  vanish  only  for 
a  peak  at  the  origin  for  the  case  in  which  B(x)  Is  an  even  function. 

In  this  case,  V1-(u)=0,  and  xo=0  causes  the  term  in  Vr(u)  In  the  expres¬ 
sion  for  ^  8(x  )  to  vanish,  so  'B(xo)=0  and  B(xo)<0.  We  are 

UA 
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usually  assured  that  the  second  derivative  is  negative  because  Vf 
is  symmetric  with 

VCoH  N,-t0)=X-  J*  >  ° 

and 

J-K  Vr(*)  =  BCo)  >  o 


although  in  some  particular  cases  we  may  find  a  different  result 
possible. 


On  the  other  hand,  for  a  peak  in  B  at  xQ/0,  or  if  B(x) 
is  not  an  even  function,  then  the  general  case  will  be  that  the 
peak  in  B  will  occur  at  a  different  point  than  that  In  B(x),  although 
again  some  particular  cases  should  be  possible  where  the  peaks 
are  coincident.  One  particular  case,  in  which 


will  cause  the  peaks  to  be  coincident.  Also,  the  more  general  case 

where  r  ^  ,  f  _ 

1  ducV-N/;  c-ocj 2irux,  ?=■  I  4u.  "KVy  s/^TTax^, 

J-'f  J‘Tt 


can  likewise  lead  to  coincident  peaks. 


Unfortunately,  it  does  not  appear  possible  to  write  a 
concise  expression,  valid  for  a  general  visibility  function  with 
arbitrary  symmetric  truncation,  which  sets  an  upper  bound  for  the 
peak  shift  in  B(x)  resulting  from  truncation.  However,  If  we  assume 
that  the  peak  shift  is  small,  we  can  arrive  at  an  approximation  to 


the  peak  shift,  which  we  will  designate  by  Ax. 

Consider  the  case  where  ^6(xo+^=0,  g^B(xo)=0.  Then 

Ax)  a.  O  =■  juV{.c*s  2-TrutA*^  -ttVc  W 2iru  b„) J 

3  t*s  2-rruXp  cciZir^  -  uVi  3<k 

*  - 14 Vy.  a,/J  x,  Z7n*tLA  -  U.\t  c©s  ^7n<x,  ^^-iru.4A1 

If  we  make  the  assumption  that  ax  is  small,  we  can  expand  to  J 

first  order  in  Ax  to  get 

CjS  2-JTUXo  -  2TAAuf  Vi  ZlfUXc 

-M.Vr  S/w2tuX*  -'4?7T4AlLZVrC*s2XUXo] 

We  can  now  solve  for  Ax: 

U.  2-TTU.Xo  -  Vv- 

vf  [Vi  s/AiZmx^  +  vr 

As  an  example,  consider  the  case  V(u)=cosuu,  truncated 
at  |u|=l.  We  have  then,  for  the  point  source  at  x  =  5, 


4*  = 


A*3 


dk 


U.  Cos7TU.  StfTTU. 


Zvj  u.Z  CoyVic 


a*  JL- 


_  ~  27T  3/N  TTu_ 


o  *  i 


Jo  <lu>  ^  [< 


2 Tj ^  <iu.  u.2,  Cwa7TU. 


C03>  TTm.  —  S/At  71X.J 


_  JU  r  -cfe  '  rffa)s*f2iru.  - 


-  jj^S^TT^  1* 
±313:  * 

UCos^ttk. 

-Hr4*  '• 


=  .  c>LSeif7 
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The  exact  peak  shift  is  calculated  as  follows.  The  exact  form  for 

SU-4)+  SU+i)] 

=  5/*  Z(x.*t)  4-  Sjft*. 

Searching  for  the  peak  iteratively  gives  the  peak  location  at  x*. 56285 
(to  five  places),  and  thus  the  peak  shift  is  .06285,  which  agrees 
reasonably  well  (to  about  5%)  with  our  estimate  based  on  the  expansion 
to  first  order  in  the  peak  shift. 


We  now  consider  peak  shifts  which  are  due  to  averaging 
effects  only.  If  V(u)  Is  averaged  with  constant  intervals  over  its 
entire  domain,  then 

A 


TfcOO  =E  B(X-^>  Tt*-» 

**.*-•0 


Consider: 


o 


■K““  L  C*-^-)  U-*)  J 

Suppose  now  that  g^B(xo)=0.  Then, 

Ai  Bj(x.)  «  BUJ  [-^2Z2f  -  ]  +f ^  wc.7fv^) 

“v*o  L 

+  BU.  tH  <*.-*•)  (X.-  J 

The  terms  in  the  sum  are  seen  to  represent  possible  aliasing  effects 
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which,  of  course,  can  themselves  cause  peak  shifts.  In  particular, 
let  us  now  suppose  that  B(x)  is  a  band-limited  function  in  j x | ^Lx . 
Then,  if  we  consider  J xQf<Lx ,  all  terms  in  the  sum  over  non-zero  n 
are  seen  to  vanish  provided  that  TST^j,  since  if  we  average  at  inter¬ 
vals  as  small  as  or  smaller  than  the  Nyquist  interval,  then  there 
will  be  no  aliasing  effects  to  influence  peak  shifting. 

With  these  qualifications,  we  may  write 

&  =  'Be*.'*  [  ]  • 


In  this  case,  then  if  B5(x0)  is  to  be  a  relative  extremum  point 
coincident  with  a  relative  extremum  In  B(x),  we  must  satisfy  one  of 
the  two  following  conditions: 


1.  B(xo)  =  0 


2  Co5  7r'rXa 

X. 


Sir*  TTTX* 
TTXo2 


~Ti«  7TtXo 


TTTX.. 


Note  that  as  T-»0,  the  condition  tawrTX0=^TX0  becomes  satisfied,  being 
exact  in  the  limit  T=0,  as  required,  for  any  xQ.  Note  that  this 
transcendental  equation  has  an  infinite  number  of  possible  solutions 
for  xQ.  If  a  peak  in  B(x)  occurs  at  one  of  these  solutions,  for  a 

A 

given  value  of  T,  then  the  corresponding  peak  In  B.(x)  will  also 

0 

occur  at  xQ.  Note  that  this  result  predicts  that  the  peak  shift 
oscillates  in  such  a  way  as  to  be  zero  for  those  values  of  xQ  which 
satisfy  condition  2  above,  for  a  given  fixed  value  of  T. 


-j  - 
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If  xQ=0,  we  must  consider  the  following  limit: 

L.M  -4r  T— 7TT)L  -  5/w7TtX.  1 

*.-*«  **  L  -rrTX.  J 

Using  L'Hospital's  rule  gives 

i  ...  C-a3  ffTXo  -  g/ffaTx0  _  L  ^  -7TT  S/h  TTTXo  “^[s/Wc75uJ 

xl-i-o  Xo  X.**®  j  =  ° 

since  the  slope  of  the  sine  function  Is  zero  when  its  argument  is 
zero.  From  this  we  conclude  that  averaging  with  fixed  TiT^,  so 
that  aliasing  effects  do  not  occur,  will  not  produce  a  peak  shift 
in  a  peak  at  the  origin.  This  result  resembles  the  fact  that 
truncation  alone  does  not  produce  a  shift  In  a  peak  at  the  origin. 

For  x^O,  let  us  suppose  that 

i;B((X.+  i.]so  if  -jfe-BOO-o. 

Then 


which  implies 

corTTTXc  -  s/*'77Tx0  *'*7TfAA  s^r T\>  Cafirn^siMirr^ 
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C«3  71  l  Xo 


-  7TT“A*  s/*  7TTXe 

X.  +  ** 


=  Sift  ITT Xo  +  7TTA*  ITT"  Xc 

tttca/*  +  2x.Ax) 


Rearranging,  and  neglecting  terms  in  Ax  higher  than  first  order  gives 


Xo  St«  TTY.  -  TTTX^  c*3  TTTXo _ 

7TTXo  <^»s  7 TTy.  -  (7TATzXoZ  +  05///7r7^o 


As  an  example,  consider  the  transform  pair 

VW--  (i-v  e.'lL  z>  ('fjr  +  irx)  e"”  *  =  Bix}. 

-5 

Now,  V( 1 u | -3. 5)^10"  ,  so  we  effectively  are  averaging  over  the  full 
visibility  function  If  we  truncate  at  |u|=3.5.  Since  B( |x|*l)S10~4, 
the  "effective"  Nyquist  Interval  is  TNI=  .5.  We  choose  a  fairly  gross 
averaging  interval  T=1 .  The  exact  peak  in  B(x)  is  calculated  to  be 
x0=. 07879.  The  peak  in  B5(x)  can  be  determined  by  searching  on  the 
computer,  and  is  xQ=.0685,  to  four  places,  giving  Ax=-. 01029 
for  the  true  peak  shift.  Our  equation  predicts  A  =-. 01976,  giving 
the  peak  at  x=.059.  a  is  a  good  estimate  for  A  but  is  off  by  a 
factor  of  1.92  as  opposed  to  the  previous  case  for  shifts  due  to 
truncation  only,  which  agreed  to  within  5%  with  the  true  shift. 

We  now  consider  peak  shifts  resulting  from  symmetric 
truncation  and  averaging.  For  the  case  of  averaging  a  symmetrically 


1 


_ _ 


-  j 
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truncated  visibility  function,  we  have 


b£u)  =  ■?*'{  =f  t  [vj  c";‘27ru’'x 

•j  -  rti„+  15l. 

=1  e-i27rM'^  V(a) 

*n«-W  T 


Then, 


i  =£  (-tfJTUO 


Treating  the  value  un  as  constant  independent  of  the  u  integration. 


we  may  write 


...  ...  -  -  ^  ^  /*^+^  r 

d*  =  ca^U‘ SirfiJiru.,,* 

-  Vl  Cas  ^TTU^)  +  L  (  Vt  simZttu^X  +  ^  «u»j£7TUv,x)|. 

Noting  that  up  is  odd  in  n  with  uo=0,  Vr  is  even  in  u  and  is 
odd  in  u,  we  see  that  the  sum  over  integrals  can  be  written  as 

i  a  H  r*-*^ 

'ST  B£ £*)  =  - 4ttE iu. (Vr s«2iru,x-  Vt  **s£7nu»x). 

T 

As  before,  we  suppose  that  ^  B(xo)=0  and  that  B^(xq+ax)=0.  Then, 

\  K  f^+%*  / 

&B.(\,+Ax)  =  o  =  -4tt£  U-J,,  T  cMVv s'* ^TTOv, U+  A* ) 


-  tc^ZlTM^  (X.+  A*)) 


which  implies 


(^2muX.  coi2Tri^bK 
+  cos^TTi^Xo  ^£7^  A*)  '  Vl 

—*  sw  Z7TlLy\  Xo  >27T7^»v\ 
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Assuming  that  the  peak  shift-  in  small,  we  expand  to  first  order 
in  ax  and  simplify  to  get 


(Vr  -  Vi  cp5  ZlTU^ 

(Vl  s"<  2ttu^x.+  \4c**2tt2ux.J 


As  an  example  for  this  case,  we  consider 

..  .  -|ul  _ l _  .  I 

VCU)-  e.  c/»5irU‘D  ^  [2TT(.ir-x)]i  r  l  +  [27T6?+x)]Z 

We  want  to  truncate  a  fair  amount  of  V(u)  and  average  with  T>TNj. 

We  choose  truncation  limits  at  |u|=1.5,  and  use  T=.6,  since 
V( |u|*l  .5) - .2231  and  B( |x|>7.2)<.001  implies  TNI=.0694  as  the 
effective  Nyquist  interval.  The  true  peak  is  at  xQ=. 49939  and  the 
peak  in  B{(x)  is  at  xQ=. 49877,  both  found  by  iterative  search  to 
five  places.  Our  equation  for  peak  shift  gives  an  estimate  of 
Ax=-. 0002374,  whereas  the  exact  peak  shift. is  -.00063.  Agreement 
here  is  to  within  a  factor  of  2.6.  Our  equation  gives  an  order  of 
magnitude  estimate. 


In  a  real  experiment,  one  would  not  know  the  true 
visibility  and  brightness  functions,  but  would  rather  have  only  the 
average  values  of  the  visibility  function  as  data.  For  the  cases 
of  averaging  over  the  entire  visibility  function  and  averaging  a 
symmetrically  truncated  V(u),  our  equations  involve  only  Integrations 
over  the  visibility  function,  which  are  just  equal  to  the  average 


'Boo 
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value  (which  we  have  as  data)  times  the  knownvalue  of  the  averaging 
interval.  We  can  thus  estimate  a  margin  of  error  for  all  peaks 
in  the  brightness  function  generated  from  the  average  value  data. 
For  the  case  of  truncation  only,  we  cannot  determine  the  integrals 
(which  now  contain  other  factors  of  u)  from  the  available  data; 
however,  this  case  is  not  of  significant  importance  since  experi¬ 
mental  data  will  be  in  terms  of  average  values  of  the  visibility 
function. 


4.6  Irregular  Averaging  Intervals 

In  the  development  of  the  averaging  theorems  in  this 
chapter  and  in  the  discussion  of  peak  shifts,  we  have  used  the 
simplifying  assumption  that  T  was  a  fixed  averaging  interval  size. 
We  now  want  to  consider  the  generalization  to  the  case  where  T 
is  a  variable  parameter. 


For  the  case  where  V(u)  is  averaged  over  its  entire 
domain  with  irregular  averaging  intervals  with  epoch  at  the  origin, 
we  have 

x, [C/J 

‘n*-n 

where  the  midpoint  of  the  n^  interval  is  now  given  by 

-  _  **-/  ___ 


■U-K  = 


■nr* _  _ 

"t  x *f«-+ lb  *•  '*<o 

Va-I  *■' 


'•tt+l 


it'  ^>0 

n.-  o 
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Then, 


m£xf  «ix‘  Buo 

iu-<o  J  _eD 


-i£7T4„(X-X') 


Because  u  ^iTfor  irregular  intevals,  we  cannot  replace  the  present 
sum  by  a  sum  over  delta  functions  and  thereby  simplify  the  integral, 
as  was  possible  for  constant  T. 

For  the  case  where  V(u)  is  truncated  symmetrically  and 
averaged  with  irregular  intervals,  we  get 

•  * 

In  particular,  for  the  important  case  Tn=T_n,  this  can  be  simplified: 

B(x')  5/rfcTx*  +  /d*'BU')  5^Xx' 

°  nt>|  *'"• 

-CosZTTll^r) 

-  Xj  BU‘)  sMcrmj  -HZZ'n^cAs^muoc  f  Jx'BOO 

*W.» 

•nn 

r*  n  no  ru**1 

SU*  BU')  +1  Sg-OK,*' 

4  s<h  x*J  4-  s/v  Ziru-nX  ^~7  le-*  a*  -  WX»  *'lj 

where  ^ 

a  7TTv,+  2”ITU-h  ~  "TT  L"^e  2LTi-J 

r» 

*  7T~7^, “■  27TXlyy  -  -7r[X+  Tcr^ 


awgv ,  r; 


For  irregular  intervals,  we  have  seen  that  only  the  bounds 


on  the  index  n  distinguish  the  cases  of  averaging  over  the  entire 

V(u)  and  averaging  a  symmetrically  truncated  V(u).  Therefore,  we 

need  to  discuss  the  generalization  of  the  peak  shift  equation  only 

for  averaging  a  symmetrically  truncated  visibility  function.  We  will 

make  the  further  restriction  that  T  =T  „  which  will  often  be  the  relevant 

n  -n 

case  (along  both  the  u  and  v  axes)  for  the  radio  interferometer 
system  in  which  we  are  interested.  It  is  obvious  from  a  review 
of  the  method  used  to  derive  the  peak  shift  Ax  for  the  case  of 
averaging  with  fixed  T  over  a  symmetrically  truncated  visibility 
function  that  the  appropriate  generalization  for  irregular  intervals 
is  just  the  previous  result 


*'*£m*x.  -Vi  *** 

5  trJ  ZmU  X*  +Vv  ^2rrni^ 


with  un^nT  now. 


4 . 7  Generalizations  to  Two  Dimensions 

The  generalization  of  our  1-D  results  to  2-D  is  entirely 
straightforward.  We  will  illustrate  the  procedure  for  one  theorem 
and  then  simply  state  the  other  results.  The  epoch  of  averaging  Is 
always  at  the  origin.  We  will  assume  that  the  cell  sizes  are 
Irregular,  and  state  the  simplified  form  for  constant  cell  sizes  if 
one  exists. 
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Theorem  4.6:  Each  average  value  of  V(u,v)  can  be  expressed 
In  terms  of  a  definite  integral  over  the  true  source  brightness 
function: 


Proof:  The  average  value  of  V(u,v)  in  the  (n,m)th  cell 


i-  "i  \  f  M  zT  I 


TH.-bT^ih  J  7m, -w  tr -"TV^ 

*.  ~z 


do-  V(u,xr) 


where  the  cell  midpoint  coordinates  are  (un,vm) 
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For  fixed  cell  sizes,  un=nTu»  vm=m^v-  We  now  state  the  other 
generalizations  to  2-0. 


The  effect  of  averaging  V(u,v)  over  its  entire  domain 
with  fixed  cell  size  is  given  by 

For  irregular  cell  sizes,  the  result  is 


l^cx, 
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If  V(u,v)  Is  truncated  symmetrically  in  u  and  v  and 
averaged  with  irregular  cell  sizes,  we  simply  change  the  infinite 
limits  for  the  n  and  m  sumnations  to  (-N  to  N)  and  (-M  to  M), 
respectively.  (Note  that  averaging  the  truncated  V(u,v)  with 
fixed  Tu  and  Ty  does  not  simplify  this  result  except  to  remove  the 
cell  size  (T  *T  )  from  the  summation.)  For  the  case  where  T  =T 

WV  U)il  Uj"ll 


and  T  =T  we  have 
v,m  v,-m 
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We  may  express  the  peak  shift  in  two  dimensions  by  (a  ,a  ) 

a  y 

where  the  x  and  y  components  independently  satisfy  the  1-D  peak  shift 
equation  for  irregular  averaging  intervals  given  in  the  previous 
section.  Note  that  this  result  exDlains  why  there  was  no  peak 
shift  in  the  y-direction  in  the  model  source  study  since  both  peaks 
occurred  at  y=0  values,  and  we  found  that  peaks  at  the  origin  do 
not  displace. 


Development  of  an  Algorithm 


Having  examined  the  effects  on  the  brightness  function 


139 


of  truncating  and  averaging  the  visibility  function,  we  now  attempt 
to  develop  an  algorithm  which  uses  the  relation  between  the  average 
values  [Vn  m]  and  the  exact  brightness  function  B(x,y)  to  obtain 
a  better  estimate  of  B{x,y)  than  does  the  standard  Fourier 
inversion  method.  The  advantage  to  this  approach  lies  in  the  fact 
that,  although  the  fourier  inversion  method  disregards  the  fact 
that  the  visibility  function  has  been  averaged,  an  algorithm 
incorporating  this  information  might  be  capable  of  giving  higher 
resolution  than  one  which  disregards  it.  The  algorithm  which  we 
will  develop  here  is  by  no  means  the  only  approach  to  the  problem, 
nor  is  it  necessarily  the  best;  however,  it  may  be  a  useful 
alternative  to  the  standard  inversion  method,  and  as  such  deserves 
consideration.  We  will  develop  the  basic  concept  in  1-D  and  then 
generalize  to  the  two  dimensional  case  where  V(u,v)  is  averaged 
over  a  rectangular  grid. 


We  have  seen  that  the  average  value  of  V(u)  in  the  n^1 
interval  can  be  expressed  as 

iaTTlU* 

_ .  ; —  PL 
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Because  the  true  source  brightness  function,  B(x),  appears  in  the 
integrand  of  a  definite  integral,  we  cannot  Invert  the  equation  to 
solve  for  B(x).  An  approximation  is  therefore  necessary. 

Since  we  are  concerned  with  band-limited  functions,  let  us 


I 
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assume  that  B(x)  is  band-limited  to  the  interval  [-LX,LX].  We  can 
then  expand  B(x)  in  the  complex  Fourier  series  form 


where  the  Cp  may  be  complex,  and  B(x)  is  of  course  a  real  function. 

If  we  substitute  this  Fourier  series  for  B(x)  into  the  definite 
integral,  we  would  be  able  to  calculate  the  definite  integral's 
value,  with  the  integration  limits  now  set  at  -Lx  and  Lx.  Although 
we  can  now  compute  the  value  of  the  definite  integral,  we  are 
still  not  able  to  evaluate  the  infinite  number  of  Fourier  coefficients 
to  deduce  B(x)  exactly.  What  we  propose  to  do,  then,  is  to  convert 
the  (2N+1 )  complex  (in  general)  average  values  of  the  visibility 
function  into  the  first  (2N+1)  complex  Fourier  coefficients  In  the 
expansion  for  B(x). 

Let  B(x)  denote  the  approximation  to  B(x)  determined  by 
this  method.  Then 

c.p  a  B4X-) 

f*-*/ 

Inserting  B(x)  for  B(x)  in  the  definite  integral  gives 

= E  cJV  e  -iSTX  (V  u") 

We  have  (2N+1)  equations  for  the  (2N+1)  values  [V  n],  ne[-N,N].  Since 
the  definite  Integral  is  just  a  number,  we  call  it  I  .  We  therefore 
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can  wri te 

[v„>  E  cP  x,„ 

Pz-fJ 

We  have  a  system  of  (2N+1)  linear  equations  in  the  (2N+1)  unknowns 
Cp.  Hence,  we  have  the  matrix  problem 


and  we  can  obtain  the  required  Fourier  coefficients: 


The  matrix  elements  are  determined  as  follows: 

p,>  J-Wi  xt-hX  c 
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The  imaginary  portion  has  an  odd  integrand  integrated  over 
symmetric  limits  and  therefore  vanishes.  Using  a  trigonometric 


identity  gives 

*-r  -  _j _ 


where 


TT>+  2T  (TT,- U*)]  ,  pf  «  K-iT^-U.]] 
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The  integral  cannot  be  given  in  closed  form;  the  series  result  is 


*TTTn4*^  ( Z'r+MZ'r* 


The  matrix  elements  1  „  are  real  numbers,  and  we  can  therefore 

p,n 

separate  the  problem  in  complex  parameters  into  real  and  Imaginary 
parts: 


The  series  which  results  from  the  above  Integration  must 
be  checked  for  converaence.  (The  alternating  theorem  definition 
and  theorems  which  we  will  cite  here  are  taken  from  Marsden  and 
Weinstein,  1980.  B  the  algebraic  sum  rule,  we  can  treat  the 
composite  series  in  [(a_  _L„) 2r+1+( s_  J-V)2r+1]  as  two  separate 
series,  and  provided  each  converges  separately,  then  the  composite 
series  converges  to  the  sum  of  the  Individual  series.  Consider  the 
following  limit: 
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Both  of  these  limits  are  zero  (see  Hildebrand,  1976,  for  verifica¬ 
tion  of  the  second  limit).  The  same  result  is  obtained  for  the  series 
in  ( 6p  nLx)  .  Since  each  of  the  separate  series  alternates  in 
sign  and  the  summation  terms  in  both  have  the  limit  zero  as  r+°>, 
then  both  of  these  series  satisfy  the  definition  of  an  alternating 
series.  We  are  therefore  assured  of  convergence  by  a  theorem  which 
states  that  all  alternating  series  are  convergent.  Additionally, 
another  theorem  insures  that  if  we  use  only  a  partial  sum  to 
represent  the  series  (as  we  will  have  to  for  practical  calculations), 
then  the  error  resulting  from  using  the  partial  sum  will  be  no 
greater  than  the  first  omitted  term. 


We  must  now  address  the  problem  that  this  method,  which 
we  will  call  the  "matrix  method",  presupposed  a  knowledge  of  the 
true  source  extent  L  in  order  to  calculate  the  matrix  elements  I 

X  P  9  *  ' 

In  the  realistic  case,  lx  is  one  of  the  parameters  which  would  be 

determined  from  the  interferometric  observations.  To  compensate 
for  this  difficulty,  we  propose  that  a  zeroth  order  estimate  for  Lx 
be  used  initially  in  the  calculations.  This  estimate  may  be  based 
on  a  lower  resolution  observation  capable  of  yielding  a  close 
estimate  for  Lx  while  still  being  unable  to  resolve  the  source 
structure.  Using  the  zeroth  order  estimate  L  we  determine  the 
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Cp's  and  hence  deduce  B^(x).  Now  the  definite  integral  theorem 


states  that 


r 


*BCX)-  V(o). 


We  can  define  an  approximate  definite  integral  relation  by 

form  which  we  can  deduce  a  better  estimate  for  the  true  source 
extent,  L  Using  L  we  can  determine  B^(x)  from  the 

matrix  method.  This  process  may  be  continued  iteratively  until  we 
satisfy  the  approximate  definite  Integral  relation  to  within  some 
error  bound  (assuming,  of  course,  that  the  Iterative  scheme 
converges) . 


During  initial  1-D  trial  cases  using  the  matrix  method, 

it  was  found  that  cases  in  which  irregular  averaging  intervals 

were  used,  especially  when  the  Interval  sizes  differed  by  factors 

not  near  unity,  gave  totally  Incorrect  results.  It  is  believed  that 

the  cause  of  this  phenomenon  lies  in  the  fact  that  a  "frequency" 

factor  Jr-  was  assumed  in  the  Fourier  expansion  for  B(x).  Now, 
n  x 

{^r— }  are  the  natural  set  of  -Frequencies  required  to  expand  a 
CLx 

function  B(x)  which  is  band-limited  in  the  inteval  [-LX,LX].  However, 
we  are  required  to  sample  or  average  V(u)  at  the  Nyqulst  Interval, 
TNi=^J-  ,  In  order  to  produce  this  natural  set  of  frequencies  for 
the  expansion  of  the  band-limited  function  B(x).  This  Is  just  the 
logic  involved  In  the  proof  of  the  Sampling  Theorem.  By  using  a  set 
of  averaging  Intervals  of  different  sizes,  we  are,  in  essence. 


"confusing"  the  frequency  Information  from  which  we  are  attempting 

to  derive  the  fourier  series  coefficients,  C  . 

P 

To  circumvent  this  difficulty,  it  is  possible  to 
reformulate  the  problem  into  one  using  averaging  intevals  whose 
sizes  are  fixed  at  the  Nyquist  interval.  Thus,  given  a  set  of 
data  containing  average  values  of  V(u)  over  Irregular  intervals, 
we  can  convert  the  data  into  an  appropriate  form  for  use  in  the 
matrix  method. 

The  proposed  reformulation  algorithm  proceeds  as  follows. 
First,  we  examine  the  truncation  limits  to  see  if  an  integral  number 
of  Nyquist  intervals,  with  epoch  at  the  origin,  will  fit  in  between 
these  limits.  If  this  does  not  occur,  we  symmetrically  extend  the 
truncation  limits  so  as  to  satisfy  this  requirement.  In  the 
extension  region,  we  assign  a  zero  value  for  the  average  visibility. 
(Assigning  a  zero  value  for  the  average  vlxibility  in  the  extension 
region  is  similar  to  other  methods  where  zero  values  are  implicitly 
or  explicitly  assigned  to  portions  of  the  u-v  plane  for  which  data 
is  not  available;  however,  other  techniques,  such  as  assigning  a 
data  value  based  on  a  Gaussian  taper,  may  be  more  advantageous. 

These  possibilities  have  not  been  Investigated  here,  but  further 
study  should  include  considerations  of  alternative  data  extension 
techniques.)  We  will  say  that  the  original  data  set  Index  n  runs 
n«[-N,N],  where  It  Is  understood  that  N  accounts  for  the  domain 
extension.  If  it  was  necessary.  We  now  "reaverage"  the  average  values 
of  the  visibility  function  to  obtain  Nyquist  interval  average  values. 
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Consider  the  Nyquist  interval,  with  indexing  scheme 
in  k  similar  to  the  indexing  scheme  In  n  for  the  data  values.  Let 
there  be  0^  original  data  intervals  lying  wholly  or  partly  in  the 
kth  Nyquist  interval.  Denoting  as  the  average  value  in  the  kth 
Nyquist  interval  whose  midpoint  is  u^,  we  have 


where  n^  is  the  lowest  index  number  (in  the  possibly  extended  set 
nc[-N,N])for  an  original  averaging  interval  which  is  contained 
wholly  or  partly  in  the  Nyquist  Interval.  Using  the  original 
[V  ]  data  values  in  this  equation,  we  compute  the  values  for  the 
"reaveraged"  data  set  Vfc. 

Now,  we  can  also  write  our  expression  for  [Vn]  in  terms 
of  B(x)  for  the  terms  in  the  sum: 

V*  -  TV  fL\'  B«')  WT;*1  e/*™'*’. 

Let  there  be  (2K+1)  Nyquist  intervals  so  that  ke[-K,K].  The  Fourier 
expansion  for  B(x)  is  ^  .  27 rP* 

=  I  crel  2V-* 

and  we  can  obtain  (2K+1)  coefficients  from  the  (2K+1)  Nyquist 
Interval  average  values  of  V(u): 
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Using  this  expression  for  B(x)  in  the  integral  gives 
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where,  as  done  previously. 
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Now,  since  the  sum  over  n  here  Is  Independent  of  the  sum  over  p,  we 
can  write  *  j* 


__  K 
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where  we  define 


We  have  now  completely  reformulated  the  original  problem  Into  a  form 
in  which  the  frequency  confusion  problem  for  irregular  averaging 
Intervals  has  been  circumvented.  We  again  have  a  linear  system,  but 


now  with  (2K+1)  equations  in  the  (2K+1)  unknown  C  . 
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For  the  generalization  to  two  dimensions,  we  start  with 


the  previous  result 


>£ 


to  Jmm 


We  Initially  suppose  that  all  averaging  cells  are  fixed  at  the  Nyquist 
cell  size.  Let  there  be  (2N+1)  by  (2M+1)  cells  Inside  the 
truncation  limits.  The  appropriate  approximation  to  B(x,y)  is  now 

_  *  *  -i£*£t  -1203. 

Btx,4)  =  SE  cf,,e  *“  e  *4  «Btx,>p. 

Using  this  expression  for  B(x,y)  in  the  previous  integral  gives 

t^l-%  t  c ,,,[/> 

Each  integral  is  recognized  as  being  of  the  same  form  as  the 
previously  considered  1-D  integral.  We  therefore  define 
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We  can  formulate  a  matrix  equation  by  trating  (n,m)  and  (p,q)  as 
the  two  indices  [  ]: 


The  solution  to  this  2-D  problem  therefore  proceeds  in  a  manner 

directly  analogous  to  the  1-D  case,  including  Iterative  refinements 

for  L  and  L  . 

*  y 


The  reformulation  algorithm  also  proceeds  In  direct 
analogy  to  the  1-D  case.  Let  there  be  (2K+1) -(2L+1)  Nyquist  cells 
inside  the  truncation  limits.  The  truncation  limits  may  have  had 
to  have  been  extended  along  the  u  and/or  v  axes,  as  previously 
discussed  for  the  u  axis  in  the  1-D  case.  We  will  say  that  the 
original  data  set  indices  run  as  ne[-N,N],  m*[-M,M],  where  It  is 
understood  that  N  and  M  account  for  the  domain  extension.  If  it 
was  necessary.  Index  kc[-K,K]  Is  the  Nyquist  cell  reference  Index 
along  the  u  axis,  and  le[-L,L]  is  the  Nyquist  cell  reference  index 
along  the  v  axis.  Let  -j  denote  the  average  value  of  V(u.v)  in 
the  (k,l)th  Nyquist  cell.  We  "reaverage"  the  average  value  data 
into  Nyquist  cells: 


t  laJ’ga.  fg  1 
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Here,  nk  and  m-j  are  the  lowest  values  of  n  and  m  such  that  the 
(nk,m.j)th  cell  in  the  original  (possibly  extended)  data  set  lies 
wholly  or  partly  in  the  (k,l)^  Nyquist  cell.  Similarly,  Jk  and  G-j 
are  the  highest  values  of  n  and  m  such  that  the  (Jk,Gi)th  cell  in 
the  original  (possibly  extended)  data  set  lies  wholly  or  partly  in 
the  (kj)*^  Nyquist  cell.  Then,  using  the  analogous  procedure  as 
we  discussed  for  the  1-D  case,  we  have 
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We  formulate  a  matrix  problem  here  by  grouping  Indices: 
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where 


-  x'  x* 


Solving  this  linear  system  for  the  ^  vector  gives  us  B(x,y). 
Note  that  a  suitable  Index  contraction  scheme  is: 


(M)  =  (k-1)  K  +  i 

(p,q)  -  (P-1)  P  +  q  • 
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CHAPTER  5 

APPLICATION  OF  THE  MATRIX  ALGORITHM 

In  this  chapter,  we  apply  the  matrix  algorithm  to  several 
one-dimensional  test  cases.  Development  of  a  suitable  Fortran  source 
code  for  the  implementation  of  the  two-dimensional  algorithm  has  been 
in  progress,  but  as  yet  has  not  been  completely  debugged.  An 
advantage  to  the  presentation  of  some  one-dimensional  test  cases  is 
that  the  results  may  be  graphically  displayed  without  excessive 
computer  time  necessary  to  generate  a  sufficient  density  of  points 
for  plotting,  as  would  be  required  for  a  decent  plot  of  a  2-D 
test  case.  Additionally,  numerical  evaluation  of  the  convolution 
integrals  needed  to  study  truncation  effects  can  be  performed 
faster  and  to  higher  accuracy  for  one-dimensional  examples  than  for 
two-dimensional  cases.  After  presenting  the  test  cases,  the 
limitations  of  the  algorithm  will  be  pointed-out,  and  we  conclude 
with  a  brief  discussion  of  a  possible  direct  application  to 
averaging  along  curvilinear  tracks  In  the  u-v  plane. 

5.1  Application  to  Test  Cases 

The  Fortran  source  code  for  the  implementation  of  the  1-D 
matrix  algorithm  is  presented  as  Program  B  of  the  Appena1x.  Evalua¬ 
tion  of  B(x)  and  B{(x)  was  performed  in  a  manner  similar  to  that  used 
for  the  2-D  model  study  (see  Program  A  of  the  Appendix).  The  only 
significant  difference  was  that  the  numerical  Integration  required 


to  compute  'B*( x)  used  the  average  value  of  the  two  Integrand  values 
at  each  end  of  an  integration  mesh  interval  times  the  interval's 
length  instead  of  the  integrand's  value  at  the  interval's  mid¬ 
point  times  the  interval.  Since  the  1-D  program  to  determine 
B(x)  and  Bg(x)  is  considerably  simpler  than,  but  analogous  to  the 
2-D  case,  the  1-D  program  is  not  included  in  the  Appendix. 

Four  simple  brightness  functions  were  used  as  test 
functions,  and  each  was  studied  for  two  averaging  schemes  applied 
to  their  respective  visibility  functions.  The  averaging  schemes, 
designated  A  and  B,  were  as  follows: 

.1 

.  .S  .Mi--*—'  .is  .S  , 

Averaging  Scheme  A  J - 1 — 11  y  1 1  1  -  { 

-l  o  | 

r,  fc  .7 

Averaging  Scheme  B  |— : - 1 — ^ — 1 — - -| 

“»  ©  I 

Table  5.1  presents  the  parameters  of  interest  for  the  test 

functions.  Averaging  scheme  A  uses  Intervals  which  are  all  smaller 
than  the  Nyquist  Interval  for  test  functions  2  and  4.  Utilization 
of  the  matrix  algorithm  is  intended  primarily  for  cases  where  the 
effects  of  averaging  are  significant,  i.e.,  ..here  the  averaging 
Intervals  are  larger  than  the  Nyquist  interval;  averaging  scheme  B 
uses  Intervals  which  are  larger  than  the  Nyquist  interval  for  all  of 
the  4  test  functions. 


Figures  5.1  through  5.8  provide  comparative  illustrations 
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FIGURE  5.1  Case  1:  Averaging  Scheme  A 
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Coincident  with  B(x)  for  | x [ < . 25 


FIGURE  5.6  Case  2:  Averapinq  Scheme 
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of  B{x),  B(x),  B(x)  and  B^(x)  for  the  eight  test  cases  which  were 
examined.  Results  are  illustrated  only  within  the  bounds  of  the 
estimated  source  extent  Lx  which  was  assumed  for  each  model  source 
in  order  to  apply  the  matrix  algorithm. 

5 . 2  Observations  and  Discussion 

For  averaging  scheme  A,  it  Is  clearly  seen  from  the 
figures  that  B(x)  surpasses  §5(x)  in  estimating  the  true 
brightness  fuction  B (x) .  For  the  case  of  two  point  sources,  we  also 
see  that  the  peaks  in  B(x)  are  closer  to  the  true  point  source 
locations  in  B(x)  than  are  the  corresponding  peaks  in  B£(x); 

•  A 

thus  it  appears  that  B(x)  givesa  smaller  peak  shift  than  does  B5(x). 
We  should  note  that  the  matrix  algorithm  does  not  yield  a  non¬ 
negative  estimate  for  B(x). 

The  estimated  source  extent  Lx  was  chosen  as  1.25  for  all 
four  cases  using  averaging  scheme  A.  This  then  sets  the  Nyquist 
interval  TNI=.4,  and  a  domain  extension  was  not  required  to  apply 
the  matrix  algorithm  since  five  Nyquist  intervals  with  an  epoch 
of  averaging  at  the  origin  fit  ins'.dv.  v..  ;•  •  •  •' 

that  the  chosen  estimate  for  lx  is  not  the  true  source  extent  (or 
even  an  effective  source  extent  based  on  our  previously  specified 
criterion)  and  it  therefore  appears  that  the  matrix  algorithm  is  not 
very  sensitive  to  the  source  extent  estimation,  provided  that  d-~a in 
extension  is  not  required. 


An  estimated  source  extent  of  1.  was  chosen  for  use  with 
averaging  scheme  B.  This  value  for  Lx  set  the  Nyquist  Interval 
at  .5,  and  a  domain  extension  at  .25  was  then  required  at  both 
truncation  limits  in  order  to  have  the  epoch  of  averaging  at  the 
origin.  For  these  four  cases,  the  choice  as  to  whether  Bg(x)  or 
B(x)  gives  the  better  estimate  for  B(x)  is  not  as  clear  as  the 
previous  four  cases.  Averaging  scheme  B  used  averaging  intervals 
which  were  larger  than  the  effective  Nyquist  interval  for  all 
four  test  functions.  The  peak  shift  In  Bt(x)  for  the  case  of  the 
two  point  sources  is  quite  severe,  and  B(x)  yielded  a  smaller  peak 
shift  than  Bg(x).  However,  B(x)  has  flattened-out  in  comparison 
to  Bg(x)  and  in  this  way  B(x)  is  less  representative  of  the  true 
source  than  is  Bfi (x) .  Thus,  we  have  a  trade-off  in  this  case, 
with  B(x)  providing  a  better  estimate  to  the  locations  of  the  point 
sources  but  with  Bg(x)  providing  a  better  resemblance  to  the 
symmetric  nature  of  the  true  point  sources.  The  most  probable  cause 
of  the  flattening  in  B(x)  is  the  fairly  severe  domain  extension 
and  the  assignment  of  a  zero  value  for  the  average  visibility 
in  the  domain  extension  region,  since  the  required  extension  here 
was  25%  of  the  estimated  source  extent. 

p 

For  the  test  function  B(x)=l-|x|  and  B(x)=sinc  x,  the 
matrix  algorithm  clearly  gave  a  better  estimate  for  the  true  bright¬ 
ness  function  than  did  B.(x).  For  these  two  test  functions,  the 

Q 

corresponding  visibility  functions  have  small  or  zero  values  in  the 
domain  extension  regions,  and  the  assignment  of  a  zero  value  for  the 


average  visibility  in  the  domain  extension  regions  was  therefore 
not  a  bad  assumption.  Note  also  the  severe  effect  of  aliasing  in 
B;(x)  in  Figures  5.6  and  5.7.  B(x)  does  notshow  this  severe  aliasing 
effect. 

For  the  Gaussian  test  function,  we  see  that  B(x)  is 
clearly  worse  in  estimating  B(x)  than  is  Bg(x).  B(x)  has  broadened- 
out  in  comparison  to  Bg(x).  Again,  the  likely  cause  for  the  poor 
estimate  given  by  B(x)  is  the  fact  that  the  visibility  function 
has  non-negligible  values  in  the  domain  extension  regions,  and 
the  assignment  of  zero  average  visibility  in  the  domain  extension 
regions  is  therefore  a  bad  assumption  for  this  case. 

From  the  rather  limited  survey  of  test  functions  which 
has  thus  far  been  conducted,  we  can  make  the  following  tentative 
conclusions,  which  require  further  analysis  to  fully  justify  them. 
First,  it  appears  that  the  matrix  algorithm  reduces  the  size  of 
the  peak  shift  interval,  and  therefore  B(x)  provides  a  better 
estimate  for  the  locations  of  the  peaks  in  the  true  brightness 
function  than  does  Bjx).  The  matrix  algorithm  does  not  appear 
to  be  very  sensitive  to  the  estimation  of  Lx  if  such  an  estimate 
does  not  require  a  domain  extension.  However,  the  matrix  algorithm 
does  appear  to  be  sensitive  to  estimates  for  lx  which  require  domain 
extension  if  the  visibility  function  has  non-negligible  values  in 
the  extension  regions  and  a  zero  value  for  the  average  visibility  is 
assigned  in  the  extension  regions.  A  possible  remedy  for  this 
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problem  might  be  to  assign  a  value  for  the  average  visibility  in  the 
extension  region  based  on  a  Gaussian  taper  rather  than  the  arbitrary 
assignment  of  a  zero  value.  Finally,  we  have  seen  that  the  matrix 
algorithm  does  not  always  provide  a  non-negative  estimate  for  the 
value  of  the  brightness  function,  and  therefore  B(x)  as  well  as 

a 

B5(x)  sometimes  fail  to  yield  physically  meaningful  (l.e.,  non¬ 
negative)  values  for  the  brightness  of  a  real  source. 

5.3  Limitation  of  the  Matrix  Algorithm 

In  this  section,  we  will  discuss  problems  with  and 
possible  limitations  of  the  matrix  algorithm. 

From  the  test  functions  which  were  examined,  it  can  be 
seen  that  reasonable  estimates  for  Lx  appear  to  be  sufficient  for 
the  calculation  of  the  matrix  elements.  The  estimate  for  Lx 
determines  the  effective  Nyquist  interval  and  thereby  also  determines 
the  amount  of  domain  extension  required  to  reformulate  the  data 
into  Nyquist  intervals  with  an  epoch  of  averaging  at  the  origin. 

The  present  algorithm  was  designed  specifically  for  an  averaging 
epoch  at  the  origin;  for  the  case  where  two  averaging  intervals  meet 
at  the  origin,  singularities  were  found  to  occur  and  the  algorithm 
breaks  down.  Furthermore,  as  indicated  by  Figure  5.8,  the  domain 
extension  may  be  sufficiently  radical  so  as  to  cause  the  matrix 
algorithm  to  yield  a  poorer  estimate  for  B(x)  than  does  the  standard 
Fourier  inversion  method.  Using  an  average  value  based  on  a  Gaussian 
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taper  in  the  domain  extension  regions  may  be  of  some  benefit  here, 
and  in  any  extent  should  give  a  better  estimate  for  B(x)  than  does 
the  arbitrary  assignment  of  a  zero  average  value  in  the  domain 
extension  regions,  which  was  the  technique  used  in  the  present 
study. 

Another  problem  concerns  the  use  of  the  matrix  algorithm 
for  cases  where  the  number  of  Nyquist  intervals  in  the  reformulated 
problem  exceeds  about  seven  Intervals.  It  was  found  that,  in  such 
cases,  some  of  the  matrix  elements  become  very  large  and  the  order 
of  magnitude  differences  between  some  elements  was  as  high  as  10^. 

For  cases  where  the  matrix  elements  have  such  a  wide  range  of 
values,  the  matrix  inversion  process  often  breaks  down  (at  least 
for  the  two  different  inversion  methods  which  were  tried  in  this 
study).  Additionally,  accuracy  problems  were  encountered, 
especially  for  cases  where  the  number  of  Nyquist  intervals  exceeded 
about  seven  intervals.  Each  matrix  element  is  calculated  by  an 
iterative  process,  and  accuracy  may  be  lost  in  the  calculation  of 
individual  elements  as  well  as  in  the  matrix  inversion  process.  These 
inaccuracies  generally  result  in  poorly  determined  values  for  the 
Fourier  series  coefficients  for  B(x).  For  example.  If  V(u)  is  a 
real  function,  then  we  expect  the  Fourier  series  coefficients  for 
B(x)  to  be  real  and  symmetric  (l.e.,  Cp=C_p);  however,  the  symmetry 
between  corresponding  coefficients  is  often  lost  due  to  the  accuracy 
problems,  and  B(x)  undesirably  becomes  a  complex  function.  Modifying 
the  Fortran  program  to  operate  in  double-precision  and  perhaps  using 
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better  matrix  inversion  techniques  may  help  alleviate  the  accuracy 
probl ems . 

The  2-D  matrix  algorithm  will  undoubtably  be  even  more 
sensitive  to  accuracy  problems  than  the  1-D  case,  expecially  for 

2 

those  problems  related  to  the  matrix  inversion.  If  there  are  N 

Nyquist  cells  in  the  reformulated  problem  (N  intervals  along  both 

the  u  and  v  axes),  then  the  matrix  which  must  be  inverted  will 
4 

have  N  elements.  Accurate  inversion  of  a  large  matrix  then  becomes 
a  further  problem  to  be  considered. 

For  successful  application  of  the  matrix  algorithm  to 
"real-world"  data,  the  above  cited  problems  would  have  to  be  adequately 
alleviated.  At  this  time  it  Is  uncertain  whether  or  not  the 
algorithm.  In  Its  present  form,  will  be  capable  of  application  to  a 
general  2-D  case  because  of  its  present  numerical  problems.  Thus, 
although  the  results  illustrated  in  Figures  5.1  through  5.8  are 
encouraging,  one  must  keep  In  mind  that  application  of  the  algorithm 
to  more  complicated  cases  and  to  2-D  functions  may  become  intractable. 
The  results  presented  here  should  be  taken  as  encouraging  evidence 
for  the  feasibility  of  an  inversion  technique  for  an  averaged 
visibility  function  which  yields  a  better  estimate  to  the  true  source 
brightness  function  than  does  the  standard  Fourier  inversion 
method. 
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5.4  The  Problem  of  Curvilinear  Tracks 

Direct  application  of  the  matrix  algorithm  to  averaging 
along  a  curvilinear  track  may  be  possible.  Time  has  not  permitted 
a  detailed  analysis  or  model  study  for  this  technique,  and  we 
therefore  confine  ourselves  to  a  brief  speculative  discussion 
(admittedly  lacking  details)  of  a  possible  approach  to  this  problem. 

If  such  an  algorithm  could  be  developed,  it  would  offer  the  advantage 
of  dealing  directly  with  the  experimental  data  obtained  along 
tracks  in  the  u-v  plane  without  the  preliminary  requirement  of 
"gridding"  the  data  Into  rectangular  cells. 

We  have  seen  that  the  response  of  the  interferometer  is 
a  sequence  of  average  values  of  the  visibility  function  along  a 
curvilinear  track  in  the  u-v  plane  which  is  determined  by  the 
relative  motion  of  the  baseline  vector  with  respect  to  the  source 
vector.  For  the  case  of  an  interferometer  with  two  elements  on  the 
earth's  surface  the  track  is  elliptical,  whereas  for  a  system 
Incorporating  elements  on  satellites  the  curvilinear  track  is  not 
a  simple  closed  curve,  although  the  "open-ended"  trajectory  in  the 
u-v  plane  can  be  calculated  as  a  function  of  time.  For  simplicity, 
we  will  consider  the  case  of  an  elliptical  track;  the  generalization 
to  an  arbitrary  track  replaces  the  elliptical  curve  by  the  appropriate 
relation  v=v(u). 

Let  £  denote  a  curvilinear  track  In  the  u-v  plane,  which 
for  simplicity  we  choose  to  be  an  ellipse: 


The  differential  line  element  for  this  track  is 


a?--** 

and  the  arclength  Sj  over  the  averaging  segment  is 

si  *Je.ds  =  \  juV~ 


Denote  the  average  value  of  V(u,v)  along  the  jtfl  averaging  segment 


by  [V]s.,  so  that 


DOSj - jTasU^  e  Bo?,yl] 


where  we  have  assumed  a  region-limited  brightness  function  and 
replaced  V(u,v)  by  the  Fourier  transform  of  Its  conjugate  function. 
As  above,  we  assume  that  B(x,y)  has  an  exact  Fourier  expansion  given 
by 


B(x,y)=2L  2  ^vf.% 
-f—  S*’" 


and  that  an  approximation  to  B(x,y)  is  given  by 


S(*l!),=£f|Uc™e 
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Substituting  B(x,y)  for  B(x,.y)  in  the  expression  for  [V]s  gives 

,  i-i  Z  c, .  f  ds 

n,Jei  7r(-J?rK-u)  ’r(-fc|-«- 


“O 


The  line  integral  can  be  performed  by  expressing  v  in  terms  of  u 
along  the  arc  £  j  and  expressing  ds  in  terms  of  du;  we  get 

,  is  -  rui  +  iui. 

^V—sruT- 

ri£k-“)  H-hz-'™) 

Contracting  the  indices  p,q  to  a  single  index  [(p,q)-(p-T)P+q] 
then  let  us  formulate  a  matrix  equation: 


/t»*\  /  \  /  <?•  \ 

(  )=(  K  c‘t.»y 


where  the  matrix  elements  are  given  by 


ruj+iu^'  /ft.a+  u.*(wV»f) 

H]/ri)-  j*t  dv  a-*~ut  T[-fc-u] 


The  integrals  may  be  performed  numerically.  As  discussed  in  the 

previous  cases,  we  assume  that  a-priori  estimates  for  Lv  and  L  are 

a  y 

available,  and  that  an  iterative  refinement  of  the  source  extent 
may  be  possible  if  such  a  scheme  converges. 


A  few  comments  are  in  order.  We  have  attempted  to 

outline  a  one-dimensional  treatment  for  a  2-D  Fourier  transform 

because  actual  data  is  available  along  1-D  curvilinear  tracks.  The 

variable  u  was  chosen  arbitrarily  here  as  the  1-D  independent 

variable;  actually,  the  variable  chosen  to  be  independent  should 

probably  have  the  smaller  of  the  two  Nyquist  interval  lengths 

(i.e.,  choose  u  if  Tu  NI<Ty  NJ).  It  will  probably  be  necessary  to 

reformulate  the  problem  in  some  way  as  discussed  above  for  the 

1-D  and  2-D  cases  with  rectilinear  averaging.  Since  only  one  variable 

is  left  independent,  it  may  be  necessary  to  use  TNI  (the  smaller  of 

the  two  Nyquist  interval  lengths)  instead  of  ^  and  1n  the 

_  x  y 

expression  for  B(x,y),  i.e. 


"B  =  I  L 


e 


-iZirTm 


Since  the  smaller  of  the  two  Nyquist  intervals  Is  used,  aliasing  effects 
should  not  occur  along  the  x  or  y  direction.  The  averaging  scheme 
could  then  be  reformulated  in  terms  of  the  smaller  Nyquist  Interval. 


To  check  the  feasibility  and  develop  the  details  of  a 
curvilinear  track  approach  will  require  considerable  further  analysis 
and  model  studies. 
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CHAPTER  6 

SUMMARY  AND  CONCLUSIONS 

We  have  reviewed  the  origins  of  the  basic  Fourier 
transform  relation  between  the  brightness  distribution  function  and 
the  visibility  functions.  For  VLBI  systems  In  which  some  or  all  of 
the  array  elements  are  located  on  satellites,  the  interferometric 
response  is  a  set  of  average  values  of  the  visibility  function 
obtained  along  a  curvilinear  trajectory  in  the  u-v  plane.  For  a 
pair  of  interferometer  elements  on  the  earth,  the  trajectory  through 
the  u-v  plane  is  elliptical;  however,  if  one  or  both  elements  of 
the  pair  are  on  satellites,  the  trajectory  is  in  general  an  open- 
ended  curve.  An  advantage  to  the  use  of  elements  on  satellites 
is  that  the  curvilinear  track  in  the  u-v  plane  does  not  close  back 
upon  itself  after  completion  of  an  orbit  (due  to  the  non- conmens urate 
orbital  periods  for  the  two  elements),  and  the  track's  coverage  of 
the  u-v  plane  therefore  increases  with  time.  In  this  way, 
observational  data  may  be  obtained  over  a  large  portion  of  the  u-v 
plane  inside  the  truncation  limits,  which  are  determined  by  the 
maximum  baseline  projection  normal  to  the  direction  to  the  source. 

Using  interferometer  elements  on  satellites  provides  the 
opportunity  for  longer  baselines  and  hence  higher  resolution  of 
source  structure  than  Is  possible  with  elements  on  the  earth.  How¬ 
ever,  due  to  a  number  of  factors  which  we  have  discussed,  averaging 
of  the  visibility  function  along  the  track  in  the  u-v  plane  becomes 
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a  significant  problem  for  systems  which  incorporate  elements  on 
satellites,  and  the  brightness  function  which  is  derived  from  the 
observational  data  may  be  seriously  distorted  due  to  the  influence 
of  averaging  effects. 

The  major  goal  of  this  thesis  has  been  to  elucidate  and 
understand  the  effects  on  the  brightness  function  due  to  truncating 
and  averaging  the  visibility  function.  We  first  gained  a  qualitative 
appreciation  for  the  nature  of  these  effects  by  studying  a  model 
case.  We  then  quantitatively  analyzed  these  effects.  Periodicity 
considerations  were  discussed,  and  we  concluded  that  aliasing 
effect  should  not  be  severe  unless  the  data  are  gridded  onto  a 
lattice  with  constant  (or  nearly  constant)  cell  sizes  which  were 
larger  than  the  Nyquist  cell  size.  We  then  investigated  the 
separate  effects  of  averaging  and  truncation,  as  well  as  their 
combined  effect. 

The  effect  of  symmetrically  truncating  the  visibility 
function  is  to  smooth  the  brightness  function  and  to  introduce 
"ringing",  and  this  effect  is  concisely  expressed  by  the  convolution 
of  the  true  brightness  function  with  sine  factors  whose  widths  are 
determined  by  the  truncation  limits.  The  effect  of  averaging  the 
visibility  function  with  constant  averaging  intervals  was  concisely 
expressed  by  multiplying  the  true  brightness  function  by  a  sine 
factor  whose  width  is  determined  by  the  size  of  the  averaging 
intervals,  and  periodically  extending  this  function  with  a  period 
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given  by  the  inverse  of  the  averaging  Interval  size.  Averaging  and 
truncating  the  visibility  function  produces  a  brightness  function 
which  bears  a  combination  of  the  individual  influences  cited  above. 

We  then  examined  these  results  to  check  for  consistency  in  limiting 
cases,  and  also  for  the  case  of  a  linear  visibility  function  for 
which  averaging  and  sampling  became  identical. 

It  was  found  that  averaging  and  truncating  the  visibility 
function  produced  shifts  in  the  location  of  peaks  in  the  brightness 
distribution  function.  This  effect  was  analyzed  and  we  produced  a 
result  capable  of  predicting  the  peak  shift  to  first  order.  This 
result  can  be  applied  to  actual  experimental  data;  however,  future 
work  in  this  regard  should  consider  expansions  to  higher  order 
in  peak  shift  and  then  use  the  computer  to  solve  the  resulting 
polynomial  for  the  peak  shift. 

Since  the  basic  theory  of  truncation  and  averaging  effects 
was  developed  in  one-dimension  for  fixed  averaging  Interval  size, 
we  then  generalized  our  results  to  include  irregular  averaging 
interval  sizes,  and  we  generalized  to  two-dimensions. 

The  ancillary  goal  of  this  thesis  was  to  try  to  apply 
our  knowledge  and  understanding  of  truncation  and  averaging  effects 
to  develop  an  algorithm  which,  by  accounting  for  the  effects  of 
truncation  and  averaging,  might  possibly  produce  a  better  estimate 
for  the  true  source  brightness  function  than  does  the  standard 
Fourier  Inversion  method.  We  developed  an  algorithm  in  1-D,  and 
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subsequently  generalized  it  to  2-0.  We  have  called  this  method  the 
matrix-algorithm  because  our  particular  application  of  the  theorem 
relating  an  average  value  of  the  visibility  function  to  the  bright¬ 
ness  function  resulted  in  a  matrix  equation  which,  when  solved, 
yielded  a  set  of  Fourier  series  coefficients  for  a  function  which 
approximated  the  true  brightness  function.  The  matrix  algorithm 
was  applied  to  a  few  test  cases,  and  the  results  were  compared 
with  those  of  the  standard  Fourier  inversion  method.  Unfortunately, 
the  matrix  algorithm  did  not  produce  the  best  estimate  for  the  true 
brightness  function  on  a  consistent  basis.  It  appeared  that  the 
domain  extension,  which  is  sometimes  necessary  in  order  to  be  able 
to  use  the  matrix  algorithm,  was  the  cause  of  the  poor  results 
obtained  in  some  cases.  In  those  cases,  the  visibility  function 
had  non-negligible  values  in  the  domain  extension  region,  and 
our  crude  assumption  which  assigned  a  zero  value  for  the  average 
visibility  in  the  extension  region  apparently  led  to  the  poor 
estimate  for  the  brightness  function.  Future  investigations 
should  try  alternative  techniques  for  assigning  a  value  in  the 
extension  region;  for  example,  a  value  based  on  a  Gaussian  taper 
may  be  useful . 

The  matrix  algorithm  was  found  to  have  some  numerical 
difficulties  when  attempting  to  apply  It  to  cases  where  the  number 
of  Nyquist  intervals  in  the  reformulated  problem  exceeded  about 
seven  intervals.  A  refined  computer  program  may  correct  this  diffi- 
ciency,  but  if  it  cannot  be  corrected  then  the  matrix  algorithm  will 
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not  be  feasible  for  "real -world"  applications,  and  another  approach 
should  be  considered. 

The  possibility  of  applying  the  matrix  algorithm 
directly  to  averaging  along  a  curvilinear  track  was  briefly  considered. 
A  basic  approach  was  sketched-out,  but  a  definitive  statement  as  to 
the  feasibility  and  details  of  such  a  technique  will  require 
considerably  more  study  than  has  been  possible  here.  If  such  a 
technique  could  be  developed,  based  either  on  the  matrix  algorithm 
or  some  other  algorithm  which  accounts  for  truncation  and  averaging 
effects,  it  may  possibly  provide  a  much  better  estimate  for  the 
brightness  function  than  does  the  method  of  grldding  the  data  and 
applying  the  standard  Fourier  inversion  method. 
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APPENDIX  A 
SYMBOL  GLOSSARY 


The  following  alphabetized  list  provides  a  summary  of  the 
major  symbols  used  throughout  this  thesis.  Symbols  which  are 
introduced  in  the  text  and  used  only  for  a  short  topic  are  not 
included  here. 


Baseline  vector  for  a  two-element  interferometer. 


B(x,y)  Source  brightness  function. 

B(x.y)  Brightness  function  which  is  the  Fourier  transform 

of  a  truncated  visibility  function. 

B  (x,y)  Brightness  function  which  is  the  standard  Fourier 

0  transform  of  a  truncated  and  averaged  visibility 

function. 

B(x,y)  Brightness  function  obtained  by  matrix  algorithn 

from  truncated  and  averaged  visibility  function. 

b  ,b  Extent  of  the  visibility  function  Inside  the 

truncation  limits  in  the  u  and  v  directions, 
respectively. 

c  ,c  Midpoint  coordinates  for  the  non-truncated  portion 

of  the  visibility  function. 

6(x,y)  Two  dimensional  Dirac  delta  function,  (x,y)=  (x)  (y) . 


Ax,Ay  Peak  shift  values  along  the  x  and  y  directions, 

respectively;  i.e.,  Ax  indicates  the  distance 
along  the  x  axis  from  a  peak  in  the  true  source 
brightness  function  B(x,y)  at  x0  to  the  corresoond 
ing  peak  in  B(x,y)  at  (x0+ax). 

$  Fourier  transform  operator. 

-£-l  Inverse  Fourier  transform  operator 

I  Matrix  used  In  the  matrix  algorithm  whose  elements 

p*  are  definite  integrals  over  the  source  region. 

1  /-F  . 

1,  as  subscript:  Imaginary  part. 


L  L  Distance  from  the  origin  along  the  x  and  y  axes, 

respectively,  to  the  region  boundary  of  a  region- 
limited  brightness  function.  If  B(x,y)  is  not 
truly  a  region-limited  function,  then  an  effective 
region  boundary  may  be  defined;  here,  we  use  the 
criterion  that  B(x,y)^.001  for  all  x,y  outside 
the  effective  region  boundary. 

NI, as  subscript:  Nyquist  interval,  l.e.,  Tu,ni  means  that  the 

averaging  Interval  along  ’  the  u  axis  Is  equal 
to  the  Nyquist  interval. 


r,  as  subscript:  Real  part. 


s 


0 


sine  x 


t 


o 


Tu,n*Tv,m’> 

T  T  *T 
1  U  ’  V  ’ 


U,v 


Vvm 


Unit  vector  in  direction  of  a  point  source  or 
approximately  in  the  direction  of  the  centoid  of 
an  extended  source. 

Unit  vector  in  the  direction  in  which  the  Inter¬ 
ferometer  is  pointed. 

sine  x  =  — 1 V-  ,  after  Bracewell  (1978). 

Integration  time  period. 

Averaging  interval  size;  subscripts  u  and  v  designate 
respective  axis  along  which  interval  lies;  subscripts 
n  and  m  index  a  particular  interval  along  the 
respective  axis;  no  n  or  m  index  Indicates  that  the 
averaging  intervals  all  have  the  same  length  along 
the  relevant  axis;  no  subscripts  at  all  Indicates 
constant  averaging  Interval  for  one  dimension  only. 

Orthogonal  components  of  the  interferometric 
resolution. 

Midpoint  coordinates  for  the  (n,m)^  averaging  cell 
(or  the  (n,m)^  averaging  arc  for  an  elliptical 
track  in  the  u-v  plane) . 


um*„,vm<„  Lower  truncation  limits  along  u  and  v  axes, 

m1n  min  respectively. 


V(u,v) 


V(u,v) 


Visibility  function;  exact  Fourier  transform  of 
B(x,y) . 


Truncated  visibility  function. 

The  average  value  of  the  visibility  function  in  the 
(n,m)th  cell. 
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(u,v) 


Ttu.v) 


w 

x,y 


xo*y0 


III 

TTu.v) 


The  average  value  of  V(u,v)  in  the  (n,m)^  cell 
expressed  as  a  function  of  u  and  v  by  using  a 
delta  function  at  the  cell  midpoint  (standard 
representation) . 

The  function  which  represents  the  truncated  and 
averaged  visibility  function;  Fourier  transform  of 
B.(x,y) ;  T(u,v}=  z  0  (u,v). 

6  n,m  n»m 

Width  parameter  for  the  model  source. 

Angular  position  coordinates  for  a  celestial  object; 
viewed  as  rectangular  coordinate  system  on  the  plane 
of  the  sky;  angular-rectangular  conversion  possible 
if  distance  to  source  is  known. 

Position  coordinates  for  a  peak  in  B(x,y). 

Convolution  product,  if  used  as  a  superscript, 
denotes  complex  conjugation. 

Indicates  a  Fourier  transform  pair:  V(u)3B(x). 

Sampling  function:  III(y-)  =  T  e  6(u-nT),  after 
BraceweTl  (1978).  1 

2-D  Box  function,  tt(u,v)=h(u)  t r(v),  where 

•i  |U|4 

Tl(u)  =  \  1  U  |  =  Jr 

0  | u | >1  ,  after  Bracewell, 
(1978).  c 
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SUM!  =  S  JF.V  I  ♦  (  02  *V  IHAGI  Jf  l  D»  VM  ID)  ) 


.17:, 


v j  i  L'  =  v  ^  !  u*oci  r a 
C  tT  (  iJP 

V  l  D  -  >/ Vf’l 
I  D  =  U 't  I  -* ♦  '  L  ”  • 

1170  C'lf  I  VJC 

'.-'jjn  './  ( r  t  j~_  .  «t  r  j  ) 

A  V  »  V  l  *  S  li  •"  V I  / 1  T  T  ^ *  r  T  V ) 

VVV U sVVV^IP 

is:  r  !  ic*s 

.as  r> i.  . v i  =  : 

’  strvi* 

1  a  7  (  M  C  U  =  i  »  v,  U 

cr  217  inco.v=i  ,vv 

5UMV<  =  SU:V-U(  C?- V^r  AL(  C  :i  C  .VM1PI  1 

3  s  j  :/ 1  =  su«v  i  ♦(  >  v  ■■/ i'M(  u  .  i  j,  vmid  n 

V  I  0  =  7  M  I  (J  ♦  i'  L  L  r  A 
217  linn  tU: 

V'‘ID=VVvI0 
Jv- 1  D=  I  0  >!j  C  L  T  A 

117  o  •iTlMJL 

•W .7  5  =  S'2.W''  /(  T  T  l  j  •  T  1 7  ) 

AViV[  =  SIOVI  /(  TTl- r  IV) 

u  li  .ci  ,j  a  .  t  ;i  \u".‘u 


;  unci.  li.  's  i  Ti*  jw-</ 

188  V1lAV:*(  LINC1  ,LI  NCI  . )  »<\V.,\/K 
V  IMAV3IL  TICI.L  I\w  1  1  >=AVCV  l 
P^INF  312t  JNtM  J,7M1rrV,AVGV,<,AVGVl 
8  12  F()H^4T(2X»Fi:.  t,5KfFK.o,SX,F10.6,5X,F11.6,lJX,El&.8,l>XtE16.8) 

W-UIE17,b;3)  LH  »  ITUtV.i,  TTV#AVjV-1»AVUVI 
SCI  F08MAM& (FlZ.S  )) 


X=  XMI N 

CD  iU  I  X  =  1 ,  N  X 

y=yx  n 

DJ  21  ■*  I  Y  =L  ?  MY 

F  ACCDS  =  CL)SU  2.’VI  )  M  IbV^X  >  ♦  (  VM*Y)  )  > 

F  ACS  l  'Is  S  I  4  (  (  2  .  Jl)‘((L,  X  I  ♦  (  V  M  ^  Y  I  )  ) 

hj(  i  x,  i  y)  =  (?  .«rTi:-rrv>-(  (  w  cv^faccos  )  ♦  i  avgv i  *facs  u  ) )  )*B3(  ix,  i  y) 

Y  =  Y  *  .  2  •»  ,x 

c  •  .i :  'ij  -:  77  Ps  — **  1’ 

Y  -  J  l  , 


1 1  K  i  "  S 

B7  5V 
t'r<  H 
t-ELu'S 


~VVs 

>)  rr,  /  /.  1  *■" 

/\  i(^C- 

C<  ,V  7X1.  * 
"B  .  ’TYc'i—'' 


Hi  „  .»;..) 

VVN{  LI  ICi~  )  =V*; 

LI  421  =LI'IC1  .  U 
i f i  iv:i .  .M  -.  dsu  ro  i c 7 
'  1FII4C2.  .Ml.DTO  r  j  It  7 
VP=V7*  TTV 
yj  T.)  IP*> 

187  /'•'  =;v  Mu  v/2.  ) 

idt>  IMD-V'/F  l  DM  TV 

.'  =  7010 
U X  I  Li  -  li U '*  I  f 1 

ns  c  1 4i  ;  l  ie 

i  i 

'J  i  .  1  lf  I  l  i.  L  I  •  I 

VVH  I  L’ =  VM  I D 
JM|D=UIMID*TTU 
UUMIl>=  J  N  I  0 


LEVEL 


1  I  h 
115 


146 


127 

127 

126 

125 


“At  W 


187 


CAfE  =  £’3313 


07/49/  .  1 


JH  41^1  MC  1)=UM 
LI  -1C  .  =L  I  jn  1  ♦  1 
TT'J/  .  J 
.  [  JO  .  =1 
V  >1  s  J  »j  I  >| 

c  1 4  r i nu- 
03 4T I MUC 

.  r  r*  j  t 

.  U  I  _  (  7 , 1  • 


-6  ) 


\  T  (  3  ‘>4  2  22  2  ?#£,:*  ***  *  ■>?  2  2  2  ******  ***'•*222  2  2  ) 

pint  126 

X  =XM1  i 

Di!  125  l  X  =  1  »  NX 
Y-Y«!  I 'I 


&  l* 


'1C  > 
WwVC 


i  ‘  />  0  > 

S  r  '■ 


1‘)  ‘ 


.  it-- 

~.~H  C 


ft 


•> 


i~  /;*-  *-*■ 
-  *-  .n 


PI  [  \  r  127fX,Y,.-Li(U,lY) 

F  j  <  5  1  T  (  5H  A  (  X=  «  F  j  .  2  «  3t  >  Y=  »  F  5  •  2  «  2H  )  =  t  E  1  6  •  5  ) 

HUT  Et  7,  127  )  X,Y,  ?C(  IX,  iY  ) 

Funri3i?i6.!'ii 

Y=V*. 76 
u.!  H  I  -j  JE 


\  -  X  ♦  mi'i 

in  'u:_ 

P  <  I  Y T  12  6 

p  :<  i  s  r  123 

P*fU  126 
P3.IMT  128 


■fi  t-  CTiC'.P  J> 


y{  PM  M  THfr  \'A  T^ei'.S 


-ptA(g  pa-uo.-’ im- 

-^cX.yJ  vJhiC'H  A&£~  /t*  QHfipTCrl  3 


?5T 

WE 


T7Y/f  -7/£T-:»/.s 

Pit  'IT  12c 
P  4  I M  T  128 
P.7I'ir  123 
P2I-ir  128 

EX"  HATi;t\  ?F  P.-AXO 

JS.  rn:  KA^UFL  LIEXAXY  MuTlNE  VA06A  - 

P-'.I  *:T  17s 


A.  t 


xaIaI  —  ~  l • 

XX<  2)  = 

STEP=.: 5 
^  AL^=  •  1 

MAXFU*J»K2: 

I PX I  X  T  =  -  l 

CAL.  /A,:6M'4»XXfFtJfSl.P»A00»MAXFUN»[P(<I  JT»W) 

XCM  =  XX(1  ) 

YO'l  1=  XXt  2) 

t'»:  iy  \x  =  t-  i.  Mr 

PtU'-ii  12.:  ,xcf  i »  ycp l » tM  ip .«x 

12  5  '''!rM>4  P-  .  CM*  T  1  AT  (  ,1  1 :  .5  ,3  H,  ,t  16  .  , 

’  ...  1  .^i[j».  .  —  «  i  •  ;  ) 

XX  t  l  Ml 

xxt  21  =.: 

CALL  VACbACM,  XXfF,r.,SfcP,  ACC  ,  MA  XF  UN » I  P'<I  1 T  r  W ) 


_  /  £  L  *■'  A 1  <  188  3A  rE  =  8.313  09/49/ 

HZ'' 2  =  XX(  I) 
y  x  (  e ) 
r  .  v  \  - 1  —  l  •  )  M 

P  ^  I  •  ’  12.1  »  X  .  *  2  t  Y  L  2*i).  2  "X 

1  2 ' 1  F.>-M\T(  25H  PEAK  IN  CUf:v»„.-:£  MT  2  AT ( , E 16. S , IH, , £16 .8, 

1  c  i. .  1  >  AiD  pI-A.s  '*  a  *  I  TD./i  =  »  '  16  .  3  ) 


i 

:  vrCTLu  r\or  peak  i  n.  peak  2 

^MT  12S 

,  p  n  i  -I  r  1 5  „ 

15'.'  F  \  ‘ *  T I  Toll  VECTOR  f  RV,‘-  L:  FT  PEAK.  TO  P.lOHT  PEAK  > 
i>!.>T=S2RT(  (  XC-  ,’-CTll«D  (  YCf'2-YCMl )  **2  ) 

HIT  \=  A T  AN (  (  Y  2  v  t. -  Y  L ^  1 )  /  t  X .  X 2-  XC M  1 )  ) 

PRINT  15.1  *PIST  , TP  ETA 

15  1  :  ,J ;  I  1  T  (  4?m  VECT^.k  A.A9..I  Tub-  FOR  SEPARATION  OF  PEAK;,  = ,  El  6 . 8,  1 OX, 

f.  IhH/ECT.’,  A\|Ol_  nRT  X  A.<  I  S- » F  7.  2 ) 

P *  I '  r  12c 
p  •:  i  i  r  l  2  s 


r 

c 


c 

c 

;  NALF  rflPTH  PET  f.  RM  I  A  AT  I  JNS 
3Ck3= . 

HINC=.  1 

R  1  H  Al.  F  =  .6  •-r’Pl'-  Ax 
L  2M ALF  =  . 5~  RB22 AX 


I 


r 


PRINT  1 28 
PRINT  123 

YH1T  VD=S>;’U  (  (  AL  OU  ( 4.  •  )  )/  (  A  AAA'  P  I  )  ) 
YH2T<U  =  SORT(  «AI.iK,(2.  )  )/  ( <\-*,AA*  P  I  )  ) 
XX(  1 1 =  <C “  1 
XX  (2  )*YHlTRl> 

l  8 ' 5  CALL  L \LuAR  (  N,XX  ,l,m-\LF  ,  jl 
ACK  IM  (  -1.  »*3M-\LF  l-P.  1FALF 

a  1  ,  <  .  -  '  j  ( r  c  x  ;  ) 

If  (  A  V  v,  X  1  .  L  T  .  .  ’  I  T,  L‘J  , 

IF<  4ZK  I.2.T.  >>  VO  T  j  l  8,  1 
IF  (  PCX  A)  19'"  ? ,  12.2,  lBl't 
18  .'2  XX  1 2  I  -  <X  (2  T  -H I '  !C 


UNCLASSIFIED 

3  -  3  , 


ARMY  MILITARY  PERSONNEL  CENTER  ALEXANDRIA  VA  F/S  17/9 

FOURIER  TRANSFORMATION  THEORY  FOR  AVERAGED  FUNCTIONS*  WITH  APPL— ETC(U) 
FEB  81  R  J  BONOMETTI 


i 


189  DATE 


b:  313 


CV4?/ 


LEVEL 


B  :.r*  L"  -  -  1 . 

o  • ;  r  i  : « ;  s 

i  8  '4  i ;  ;"  =  » i 

<.<(;> -  >:xt  -  i: 
r,wKa=-i . 

r.  L--.5 

l  H*  1  I  F  (  3:-;3  )  1“  6,1^5.13.1' 

la^  x < { i  )=xx(2  i+m  iu 
;ckh=i . 

31’  T3  15.5 
1 8 3  hi  ic=Hive./r. 

XX  (2  J  =  XX  (2  )  +lil  ,C 
3CK3=l. 

•J  I  j  18.5 

1825  YiNII'=xx  (2  )  -rCVl 
PUNT  ly:?,YlVIt 

18  ''9  F:HX-U(24H  C'JHPuiElT  1  HALF  WI  0  TH  =  ,  E  1 6. 8  > 
XX(1)=XCH2 
xx<  2)  =  yi  z  mu 

cUX3  = 

HI  JC=  .  1 

1  3  i  5  CALL  ^  A  L  C  9  8  (  ;*XX»i-HALF»  O) 
iZi\2-i  (  -1.  )  -i;  t  ALf  I  -8  2H.4lF 
A8C<2*AW3(  8CK2) 

IFlAHCKZ.Lr  T:j  i.5  2C 

IF< 30K2.0T.3] 30  Tu  1811 
I F  C 30K3 ) 1812, 1812,  1814 

1.812  xxm*xx(2)-nnc 

0CK9=-1. 

3D  TD  1915 
1814  HlHC=HtNC/2. 

XX(2»=XX ( 2 ) “H INC 
3CK3=-1. 

03  T  J  1815 

1811  IK  t:CK«)  1316,  19i.6,  18  18 
1816  XXI  2>=XX(  2)+Hl‘j: 

PCK B  =  1  • 
or  TO  1815 
1319  hr,;=HiNC/2. 

<;<  { :  i  -  xxi z )  *r i  i. 


182:  Ytv,iu=xx(2)-Y09t 

P'11 ‘IT  1821  ,  Y 2  v  I  0 

1821  FuHMAT(24H  CDVPUMlMT  2  HALF  Vi  I  DTH=,E1 6 . 8  ) 

pam  128 
P2.INT  128 
P9I.NT  128 
?1HI  123 


>  i _ / 


-c7-t 


1‘  *,  1  4 


iyu 


J.UC  a  3j313 


OR/tN/i 


c  co'-'^uis  ii  with  exact  i:ixtyi 

■-■  ' :  r  n 

147  I  -.M  ,  r(  =•■  lj!33  3«**«****a*3J3  33) 

P  l  i  I  12  3 

print  i?e 

l  [  F  (•  = 

l-  r  r=  - 
x  1 1  ncr=.2 
yli  n:r  =  .  2 

S  •  U  .\=.  .4 
1  P'  71  ,  1=1,4  1 

XX!  1)  =  t  {-1,  )  *XL1  M  (  I-i)  *  XL  I  NCR) 

Oil  711  J  =  1 ,  41 

X  X  (  2  )  =  (  !  -  L  •  )  *YL  I  ♦  (  (J~l)  -YL  I  NCR  ) 
w  4  L  L  w  h  H  (  N,XX,C3AR,  .t ) 

X  =  XX!  1  ) 

Y =  XX!Z ) 

C'NOTF  HOAX  IS  NE  GAT  1 7c  OB  AH  i~  ALLY  ON  RFTURNFROM  CALC88,S3  *ADD*  TO 
C  1 ( X, Y  ) ,  N  3  T  lUi’  T  <  AC  I 

ri,IFF  =  !U>IFc*{  ( ‘M  A,  Y  )  ♦  t  M  )  -  SM  AREA  ) 

3  1  ;2  l  F  =  ■>  ♦  (  l  (  r  (  X,Y  I  ♦  :  4R>*«2  )  *S  P  All  ~  A  ) 

WR  IT!  7,71*+)  X,  Y,  1  Jir-F  ,.*;>••  IF 
714  FORMAT  (  2  (  F  4 . 1  )  »  t  .o  .  £  )  ) 

7L1  CONTINUE 

71  CONTINUE 

PRINT  712, P  D I F  f 

712  FORMAT (41 U  THE  TuTaL  Dt/IATION  OF  BHAT  FROM  3EXACT=, E16.8 » 

3[»N3RM  =  6DIFF/64. 

PRINT  72j»  BDNORM 

72  3  FORMAT ( 3 IX , 1 1 HnLKM AL IZFD=, E16 .8  ) 

PRI  NT  713,1’  SUC I  F 

713  FORMAT  (  49H  ThL  TjTAL  SQUARED  DEVIATION  OF  BOAT  FROM  BEX  ACT  =  ,  cl 6  .8  I 
KSN0RM  =  8SQniF/6  4  . 

PRINT  72 1 , H SNORK 

721  FORMAT!  3  3  X , 1  IHWMAlI  ZED=  ,£  16. 8» 

PRINT  .28 
P  R I N  T  128 


w  ' 

C 


C  C  OF  PARI  S  IN  WITH  BTWID(X,Y) 

WU  Till  7,7147) 

7147  FORMAT!  35H44444*<‘*#**'i«’***44444********:*',‘44444) 
PRINT  128 
PRINT  1 26 
5  < I N  f  123 

-  i  - 

3 Tjwjf  = ; 

XUNCRa.25 
YL I  NCR  =  .c5 
SHARE A* I .23 )**2 


i-EVfc'l  21  MAI  i  191  '  >  a  I E  =  80313  37/4  1/ 

u  'l  7  71  3  1=1,11 
XXlil  =  (I-l)*XL|-.C'< 

Cv  7  7 11  J*!  ,11 
x.<  tr  )  =  (  j-i  )  yl  r 

C  ILL  CALCBi  (N,XX,biiAR,G» 

X=  XXI  1) 

y  =  xx(2  ) 

0  !OTF  r.r,M  IS  rjcr.ArWf  bUAiR  *LLY  Cl\|  RETUfUFrjY  CALC33»S3  *A03*  T3 
;  3  (  X  ,  Y )  ,  J  J  T  SUBTRACT 

2TDIFF  =  bUMFF*  (  (811:101  I  ,  J  )  ♦  ?  C  A  R  )  «  S  M  ARE  \  ) 

i  btsudf^otsodi  ♦  m»-i .*  in i,  j  >  ) **2> *smarea  ) 

VHITE(7,7714)  X,Yt,-.TC:Ff  ,  RTSUCF 
7714  FORMAT  I  2(F4.  I)»2(rl6.fc)  ) 

77ii  chtinje 
771'*j  CO  NT  l  jUH 

PRINT  713.BTDIFF 

715  FORMAT  UGH  THE  IoIAL  tu/urn)j  OF  BHAT  FROM  8TWlD  =  ,!ii6.8) 

UTUNRM  =  BTOlFF/(2.3:**2) 

PRINT  722 ,  R  TO  NR  FI 

72?  c  j  R  M  A  f  ( 27X,  1 1H  < J<* AL  l  Z  _u  =  » -  16.8) 

PRINT  7  1  o  »  .»TS  . !'( 

716  F  ip,  “l  -A  T  f  4  6H  T-t.  T..UL  S«U\  •  =  L  DEVIATION  GF  ilH  AT  FROG  0T  ri  I  D= ,  E  15  .8  ) 

l’is'k^bi  s:df/(«  . j  -Cl 

PRINT  723*  iiTS'l  <K 

723  FORMAT!  37X,11M  JORf  AL I  Z  E  C=  ,  FI  6 . 3  > 

PRINT  128 
PRINT  128 


C  CHCCn  FDR  SYMMETR Y  ACjUT  Y  AXIS, IE, ARE  SOURCES  MIRROR  IMAGES  OF  EACH 
C  OTHER 

I  $  YM  1  = * 

u j  i  ix«i«3i 
xxum  ix-i  )/i : . 

00  1  .1  I Y=1 ,31 
X  X  (  ?  I  -  (  I  Y—  1  >  /  1  . 

.  --  ,  -  f  •'»  ,  f 

:  <>-r.  ri  ;*•*  ■■  i  •  • .  i  ,r : •/ Cc  :r  tm'::;i  respect tv:  >mz 

C  AUT  0  ilL  Y  me  I  R  DU  F  ;_R_  JCt:  I  -  l  F’PQR  T  AN  T  HERt 
UI1I*XXILI’|-:.) 

CALL  C  \LCHLH\,XX,„hAT2,G) 

'  A30IFF  =  A?SI8HAT1-3HAT2) 

;  IF  (  A'30 !  FF  .L  T  •  .  *■  1  1 )  GO  TJ  1 2  2  8 

ISYMl*ISYMltl 

10  8  IF! XXI 2)  »  1  1,1,.  1,1  ,1 
1013  XX(2I  =  XX( 2  )  *<- 1.  ) 

GO  TO  1,11 
13  1  CONTINUE 
1003  C  HI  I  k'F 

print  i : : f> 

1006  F  IR.  .in.  H  jY:  .  .  •  .<  Y  1  .!»•  OCX,  REFLECT  lU.i  SYMMETRY  I  *.J  Y  AXISI 
IF!  I  GY v  l.Lc.  1  )  G  J  I  J  J4 
PRINT  1C35.ISYH1 

1035  FORMAT! 52H  ASYf'f'LT RY  WAG  FOUNO  BETWEEN  THE  2  SOJRCL  COMPONENTS, 


LE/LL  i.1 


^  A  I  i 


192 


OAIE 


0.313 


i 


•1  7 n  ISf “l=,  14,/////  ) 

ii'  rn  i.lt 

13.4  n  1 1  4 r  1  .  .7 

i::7  f  ;aHAT(3.H  'JO  SIi.^lFICA^r  ASYMMETRY  WAS  KJJ^D  BErwCi'J 
C  l73300P.Cr  CO- I'uv  JTS, /////) 

1319  OJ'.MY=l. 


MINT  129 
(Mi  <7  12  9 

3  CHCCA  MEL. Cl  I. 14  CYV^RY  fj<  CASH  SOURCE  TKUl  THL  C  JMM  3N 
ISYM2*; 

on  1)20  I Y= l , 31 
XXI  2  >  =  I  IY-l)/lj. 
uli  1-21  IX  =  1 1  31 
xxui  =  i  ix-i  )/i 

1031  CALL  CALC03(M,XX, 3HA rifO> 

XX  l  2  1=1-1.  )  "XHUI 

CALL  C  ALC9B  1  N»XX  »orl  AT2 , 0) 

AE0IFF  =  AUS(  NHATl-fOAT*-) 

IFIAB1MFF.LT.  .C-DOO  Ij  l.:28 
I SYM2= I SYV2+1 

IC2  8  IF( XX(  1) ) K  3. ,1  .21 ,  l  :21 
1030  XXI 1 1  =  XX  C 1)*(-  1.  ) 

gj  T'j  i :  3 1 
1021  CoJTI'JUt 
IjPc  C.TJTDUF 

O  '. :  :  r  :  .  ■> 

*  • '  • ,  i  '  *  r  (  '  •  v  •  v  ,  ‘  r  i  ' T  l  '  v  1  •  “  *■  *  y  r  :  •  »  j 

t  it- 1  r  -.v  '.i . . :  .  1 1 >  i  .  . 

PUMT  U<l5,  ISY'9^ 

1025  FLKi’AT (2  7H  ASY «*?'*_ T  1Y  WAS  Ft'UND,  ISYM2*,  14, ///// ) 

GU  TO  1.2^ 

10?4  P  ’  14 T  1 C  2  7 

10>7  F  oTMAT  ( 35H  MO  SIGNIFICANT  AS  Y  MMLT  RY  WAS  FUJMD. /////> 
1029  ou’:my=i. 

P  <  I  '4  T  -23 
PHirjr  128 


111  CG'ITI  1JE 


TH: 


AXIS 


t.F'/ZL  .  1  CAL  CFG  194  L'A  Tt  =  ti  -  3 1 3  09/49/1 

;t'T  I , ;  CALCFj(  'if  <»F  f  /) 

CT’R  N  VR'AVG  (  <  ?,/:>»  *Jl'\r,{2  5,25)  tTTUU(25)  ,  TTVV1Z5  )  ,UUM(20  ) , 

E  77  U  >  )  t  ‘IJVAX,  ,7  .X 

:i  ioic  i  ot2)  •  a  c .  ) 

M=  2 

-m  =  3.  1 4 1 o 

;  f  is  - wit  x,y) 

F  =  . 

o  ( : )  =  • 

1  0(2)=. 

Dll  1  ■  I=1|MU!'-'X 

L>  *  1  ?  i  J=1 

;  k*  <  I  *1  r  .  2  tUU'l(  I  )  tTTUU  t  1  )  ,V/V,(  J),TTV\M  J  ),  VRE  AVG(  It  J  ),VIYAV0(  I,  J  » 

^0!?'.  F  jlv.ATI  2X.F  1 . .  G,  j:<,F  1 . .  o,  ‘>X,F  1  3.  6 , 5X,F  10. 6,1  OX,  €16  .H  »5X  ,  Eib  .3  ) 

f  ac* c. . »  jj  ( i )  r  f ,  / (  j  ) 

A  i  j  =  (  (  UL'4(  l  )*X(  1  ) )  ♦(77f.(  J  )-~X  (2  )  )  )*TPI 

F  =  F-(  F \C*(  {  V*?AVG<  I  ,  JIMEOSt  ARG)  >  >*  ( V I  FA/GI  I  ,  J )  *  1 3 1'J  (  ARC )  >  )  )  ) 

G(i  l  =  C(  1  )♦(  ( T3  I 1  F  A0I:UU  1 1  I  ))  *(  < VRlAVGI  I,J)“tSI'UAFG))>- 
L  (  71  >;/.7G(  I  ,  J )  -  (  C  L  (  Alu)  n  n 
'i(  :»  =  i  t  ?)  ♦•  (  (  T-M  -  F*.s«  Wf  (  J)  )  *  U  VRE  AVG!  I»J)*(SI'M  ARG)  )  >- 
£  (7  (■'  Wll  I.  J  )  (  u  .  .  (  A  3  j  )  III) 

i •; i  c.  (ri  ;j: 

10  Z'.'UIJU? 


O' 


LL'/rL  I’l 


C  \LO 


t:\rr  -•  a-iij 


07/4  it  1 


-bW. 
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sl''Ujutimc  calo.-:* f- 

i  ■  x  >'i  V.-v.  ;i  .)»vr  .15,£5)f  FTUUU5)  ,TrVVU5)  ,UUMI25>  » 

a  /  /"(  2  5  )  ,  ]'.•>  -V,  .V  '.X 
:i/.-  ISI0N  S  €  2 1  v  Kft  r  » 
i=: 

pi=o. 

T  P  I  =  2  •  “PI 
F  -  , 

,  di  = 

:o  l  i  =  i  ,  tbv'v< 

Uj  l.l 

:  P  U  ‘iT  i  .?  ,Ul\<  I  )  t  T  T  U'J  (  l  ) » 7  7  V  C  J  )  »  T  T  7  7  I J  »,  7RFA7G<  I,  J  1»VIM4V0(  I,  J  » 

ao’  FD<XAT(  2X.F1  .oo'UM  .  .  f.  o<»F10.6,5X,FU.6.l0Xt6;e».8t5X,E15.3  ) 

FA>t  .  -TlU'il  n  ■  |  WV(  J  > 

(  (  0U  l(  I  )  *x  C  )  >♦(  77  r  (  J  )*xt  2  )  )  »*TPI 
C-F-(F  AC*  t  t  VXf  '-7^1  i  «  J)  *  COS<  ARGH  )  -  ( V  I KA7G  ( I  .  J 1  *  (S  1^1 1  ARG )  ))>  ) 
l  C?4l!'IU' 

io:  c.  ,ri  ij." 

x  TUX  -1 
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APPENDIX  C 


Proaram  B 


n  n 


LEVEL  21  MAIN  197  OATE  =  80325  13/27/46 

DIMENSION  L!C7) ,M( ;7) 

DIMENSION  Z ( 07»  C  7) 

DIME  IS  I  ON  ZZ!C7,C7),ZZZ(.  ',07) 

DIMENSION  VR1  (3C  )  »Vll  (3U 
DIMENSION  CR  (  30  ) ,  L I !  3C  ) 

DIMENSION  TDATAI3: >• T  C 3c ) 

DIMENSION  A<30) 

DIMENSION  VR!25),/it 25),VRR!25),VBI !25) 

DIMENSION  BRI100 >, C  II 10 0) 

DIMENSION  TUSE  T I  3 ! ) 

,  VKEALf U>=( 2./W)*!C0S! 2.*PI*U) )*EXPI It -1 . > *PI /W>*« U**2 )> 

V IMAGI U )=0*U 


W=10. 

N=7 

PI=3.1416 
POWCK=IC.**70 
READ  1, XL  »UM IN, DEL  TAU 
FORMAT  13 (E16. 8 ) ) 

PRINT  8,XL,UMIN,CtLTAU 

FORMAT  C4H  L  X=,  E  16.  6, 5X,  5H UMI N*  , E 16.  8. 5X ,7HDE LT AU  =  , E16.8  ) 

READ  2. !TDATA(J),J=1,N) 

FORMAT! E 16. 8 ) 

PRINT  9,  I  J»TDATA(  J  )«  J=l,N) 

FORMAT  C  3H  T  1 ,  12, 2H )  =  ,E 16. 3 ) 

PRINT  6 

FORMAT!//////////) 

PRINT  85 

FORMAT! 2jX,14HAVE*AGE  VALUES,/) 

THE  *Z  MATRIX*  IS  THE  INTEGRATION  *1  MATRIX* 

UONE=UMIN  - - - - - - - 

UTWO*UONE*DELTAU  \ 

00  3  J=l,N  1 

T  t  J) »T  DAT  A  (  J  )  ^ 

TUSET1  J)*T(  J)  W£.jU\v*c 

SUM VR=C  MU-JT5  Cf 

SUKV 1*0  rr 

MN*T!  J)/DELTAU*.l  V<.*a\  /N 

CO  4  I NCRU=1 ,MN 

TC  I  I  H’T  H  M  t  t  '  *  •  H  1  ' 


l  IA2 
1111 
1114 
1113 


u one cocci 

IF!UTW0)1113, 1114,1113 
UTWJ=. jCGCEl 

SUMVR=SUMVRU ! ( VREAL ! UTWD ) +VREAL ! UONE ) ) /2.)*DELTAU> 

SUMV  I*SUMV  I  ♦(  1  !  VIMAGIUTWO  ) *VI MAG! UONE )  )  /2.)*DELTAU) 

UO.NE=UONE*-OELTAU 

UTW3=UTW0+DELTAU 

CONI INUE 

VR1!J)=SUMVR/T(J) 

VI  1!  JMSUMVI/T!  J) 

PRINT  180,J,VRl!J),J,VI  l!J) 

FORMAT !4H  VR!, 12, iH )*, E 16 .8, 15X,3HVI! , 1 2, 2H ) * ,E 16. 8 ) 

CONTINUE  —  - - 

PUN,  j 

TNI-1./ !2.*XL)  e’vflt.yjfrrg  LyQu.sr  in 

RCMRaf  !-i.*UMIN)/TNI)*.5 ) 

NCHX*RCHK  \ 


LEVEL 


VAIN 
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CHK«*CHK-NCHK  ,r  ’  -’rTC>J  •-> 

rr<CHK.LE..rct:  )gu  to  esi  -  1 

RANGE  ex  TEN  SIP':  -  t> <'  xrEf-S-.-i  'l  lj  A 

NIT  or*  I  I  l-2.*tHI  NJ/TNI  )♦!.  U.995 
NMAX*NI TOT 
PRINT  853»NtTCJT 
FORMAT ( 7H  NITO T=, 14) 

D=. 5* I  (NIT0T*Tm1  )  -•- (  2 •  *UM IN  )  J 

UVINEW=UMIN-D 

PRINT  891, UM INEW 

FORMAT I8H  UM I NEW  =  , El 6 .8 , / //// ) 

UMI N=UMI NEW 

DATA  SET  MODIFICATION  BY  EXTENSION 

NPLUS1=N*1 

CO  890  1 1  =  2 » NPLliSI 

I  lMINt=  1 1-  1 

VRI I L )=VRl ( IlMINl) 

VI(I1I=VI1 { IlMINl ) 

T I  1 1  MTDATA  {  I  lf'lNl) 

TUSETI II >-T  C II ) 

CONTINUE 
VRI 1)=0 
Vlli  )*C 

NPLUS2=NPLUS1+1 
VRI NPLUS2) =0 
V  I INPLUS2 )  =  0 
Til) =D 

TUSETI ll=Tl 1) 

T INPLUS2 )=0 
TUSETI NPLUS2 )=D 

S0^T3  850 _ _ _  ^  ^  /'k'r  ' 

NIT0T= I  I  1-2. >*UMIN)/TNI I+.995 

NMAX=NI TOT 

PRINT  853,  N I  TOT 

NPLJS2=N 

DO  852  I L=1 ,N 

VRI  IL  )  =  VR  II  IL) 

VII  IL  )=VI1I  IL ) 

C'.  NT l NOE 
NC  1=  1 
NS  I  =  1 

TSUM=9  s-vhh-  oP  £€Fc£'ys\,  Sr 

DO  600  NTNI*l,NITCT - — - '  -n,  pfT  A  t'fcr 

PRINT  951, TBIT,  TSUM  ....  ^  ~ 

F0RMATI6H  TBIT  =  ,E16.8,8H  TSUM*,E16.8)  p  *  , _ rX/ 

UU~TNI  c)F  NlGM'ST  lNic.i.v. 

IFI TINCT).GE.UU)GC  TO  601 
AINCT)=TINCT)/T II 
UU*UU“T I NCT ) 

PRINT  954 ,UU,A( NCT ) ,T I NCT ) , NCT 

FORMAT  I 4H  UU=,E16.8,9H  A  I  NC  T) » ,  E  16.  8, 8H  Tl NCT) =,E16.8,4HNCT*,I4> 
IFINCT.GE.NPLUS2I30  TO  605 
jcr='r.r»  i 
GJ  TO  602 

Al NCT) *UU/TNI 

PRINT  954,UU,A(NCT),TINCT),NCT 
NST0P=NCT 


£* TEA'S/ >■/■>/  L,/  •*> 


£tffV  ^"i> 


c,p  AeFtZ-^t XrT:^ 

-  i77Kr  £  Fer  *,•*.  v<“> 

..  hKP  2  At/Ofc,*.  I * 

(!)  p  NyQo.sT  //^TBCVA1-  -> 


LEVEL  21 


MAIN 


199  CATE  =  80325 
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i 

i  *r— 

inT 


V8RI NTN I ) =0 
VnifSTMIl— 0 
DO  670  N2 I N I T=1 ,NITCT 
Z  I N  T  N I  »  N2  I  N I  U  =  0 
CONTI NUE 

PRINT  952,NST,NSTGP 

FORMAT { L8H  ENTERING  7GI  LU0Pt6H  NST=,I4,BH  NST0P=,I4) 

00  701  NJ=N$T »NSTGP  .  .  ,= 

VBRtNTNI  >=VR(NJ  )  ‘A(  N  J  )  ♦  VUR(  NTN  I)  —  VF3.i.  &  v',‘-v*  /,  *  ,',^r 

V8I I NTNI )  =V  1 1 NJ  ) ’SAI NJ  )  *VB  l  (NTN  l )  . 

N 2=  1  -  •  -  •:  L.-.:  .  V'3X  TJf£  '"*'«*•  1  rh'T 

NN=  1 1-N  ITOT  )/2  ■  2^  A-t/irP-oL 

NNCHK=NST-N J 

UNI=UMIN*TSUM*(  TUS6  T(  N  J I  /  2.  I 
PIT  =  PI*TUSET(NJ) 

A  1=1 12.*PI  )*(  (NN/l  2.*XL»»-UNin 
ALFA=P IT  +  A1 
9ET  A=P I T-A 1 

al=xl*alfa 

BL=XL«‘8ETA 

S I  GM=  1  .  *  f 

SM=AL«-BL  K  _ — 

NOOD=  3  c »V“« 


/(  'V«A 


POM  12=0  * 

IFIPOH12.LE.POWCK)GO  TO  50 

P0WAL=1. 

POHBL= 1. 

DO  799  l,NCNT  =  l,NUDD 
POWAL=POWAL* ( AL/ I NCNT ) 
POWBL=POWBL*(BL/I\CNTJ 
CONTINUE 

P0WER=tP0WAL*P0W8L)  /MUDD 
GO  TO  25 

POWER=l ( ( AL)**NOOO)+I (BL) **NOCD H /NOCD 

P0W12=P0HER 

SIGNCK=l. 

IF|N0D0.LT.54)G0  TO  22 

l  F  (  PQV.’ER )  2  3  »  2  ?  »  2^ 

S  i  -J 

PL«E*»POwER*(- 1. ) 

POWRLN=ALOG( POWER) 

DO  651  LNFACT=2  « NODO 
XLNFAC=LNFACT 

PGWRLN=POWRLN-( ALGG(XLNFAC) I 
CONTINUE 

P3WLR=EXP(P0WRLN)*SIGNCK 
GO  TO  25 

00  650  IFACT=2  »NODO 
POWER = POWER /IF ACT 
CONTINUE 


-r&jZfl) 


«,  c^} _ rl.wl 

g  &&«'■***  c 

'JZ  ,o*»< 

,  Le*s 


,c  Cf^£ 
.  c>o  I 


SIGN=S IGMM-1. ) 
SMINCR=SlGN*POWER 
SM*SM*  SMI  NCR 
AIMCR* ABSI SM  INCR I 


f 


LEVEL 


7 


60 
7  O'J 
702 

80 


601 


C 

600 


6 10 
C 

777 

800 

70 


90 


91 


61  1 


>01 


506 


21 


ma  r  i 


200  DATE  *  80325  13/27/46 


IFI  AINCR.LE..001  )GO  TO  7 
N00D=N00D*2 

io  rj  5 

Z  (  NT  'll  ,N2  >  *  ISM/P  IT  )  *A  INJ )  *Z  INTN  I  ,N2  ) 
ZZINTNI , N2 ) =Z ( N  TNI ,N2) 

IFI N2.EQ.NITOTIGO  TO  60 
NZ=N2*l 

nn=nn* 1 


GO  TO  21 

tsum*tsum*tusehnji 

CONTINUE 

PRINT  80,(NTNI,JN2,Z(NTNl,JN2),  JN2=l,NITOTI 
FORMAT  C  3H  1 1, 1  2, 1H ,, | 2, 2H I  * ,E 16. 8) 

RRCHK*T INSTOP )-UU 
IFIRRCHK.GT..CC1IG0  TO  600 
NCT=NCT*l 
NSMNSTOP  +  l 
GO  TO  600 

TSUK*  TSUM-TUSE  T I .NJ  > 

T (NST0P)=T(NST0P  l-UU 
NST=NSTOP 

NOTE  THAT  NCT  REGAINS  THE  SAME  HERE  S  MCE  THE  CURRENT 
T  OLD  INTERVAL  STRADDLES  THE  NEW  TNI  INTERVAL 
CONTINUE 


PRINT  6 

PRINT  610, { J.VRRI J) , J,VBI I Jl ,  J*1,NKAX> 

FORMAT (5H  VBR ( , 1 2, 2H I  =  , E16. 8, 1CX, 4HVB II , I  2, 2H )  =  ,£16. 8 ) 
PRINT  6 

LINEAR  SYSTEM  SOLUTION  FOLLOWS 
PRINT  800 

FORMAT  1 24H  LINEAR  SOLUTION  FOLLOWS,//! 

CALL  MINV(Z,NMAX,D,L,M)  - 

LLL  =  l 


CALL  GMPRD ( Z , VBR  ,Ck , NMAX, NMAX  »  LLL ) 
PRINT  90,(J,CR(J),  J=I,NMAX) 
F0RMATI4H  CR (, I  2 , 2H )=, E 16 . 5 ) 

CALL  GMPRD ( Z»VBI,CI»NMAX,NMAX»LLL  ) 
PRINT  91,(J,CI(J>,  J*l,NMAX> 

FORMAT  I 4H  C  I  ( ,  1 2 , ZH I * , E 16 . 8 ! 

P  6’  I  N  T  6 


c  /-tlc/v*.  >{'7tr2i 
th '.'Z  A'be  /tU\rfl)K 
lr 


r  cs. 

Sc’i.^'T'S  Pi 
>/»/ 10 

Cr 

c  -  ;  /£•*■'  ‘L  ^  S 


PjfcAL 

<n*- 
t  & 

Ci 


’  ’  l  Oil ,  (  (  Jl , ,  L  (  Jl ,  j  ' ) ,  j i  =  1 1  ..".ax  ) ,  j2=i,  < ) 

FC’.'MT(3H  Z(,  12  I  Z,  :••)=,£  16.0 

CALL  GMPRDIZ, ZZ , ZZZ , NMAX, NMAX, NMAX > 

PRINT  6 

PRINT  611,1 (J1,J2,ZZZ(J1,J2I ,  Ji*l,NMAX!,  J2*I#wMAX) 
PRINT  6 
PRINT  6 


jLunV,  C-Mi’Zfc 

t\  L-c  i- 

TH?  S C.ICK- 
TifiC  sva  - 

At.  .-fiijC? 

PriC 

rctl  Th*- 


NX=8l  _ 

PRINT  501, NX 

FORMAT ( 23H  NUMBER  OF  X  INTERVALS*, 14! 
DO  502  NXX*l , NX 
BR I N  XX ) *0 


cf 


BKNXX  )*0 

C r N T  r  IUE 


N  <X*  1 

XM-4.M  .  1  * ( NXX-  1 )  »  I 
00  503  NNX*1,NMAX 
N INT*I I l-NMAX)/2)MNNX-l) 


LEVEL  21  MAIN  201  DATE  =  80325  13/27/46 

XF4Ct=C0S( <PI*NINT*X)/XL> 

XFAC2=S  INC  {  P  I  *’i  I.'iT  *X  )  /XL  ) 

ftRI  N<  XI = ( CR ( N\ < ) *XF  AC 1 1 -I C 1 1  NNX  ) *XF AC2  ) *BR I NXX I 
ft  l  ( MXX ) * ( C I ( NNX ) *XFAC l ) ♦( CR ( NNX  I  *XF  AC2 ) *BI ( NXXl 

503  CONTINUE 

aMAS=SQRT<  (BR<NXX)**2)M6!<NXXI**2n 
PRINT  504*  X»BR( NXX) ,X#BI(NXX>  »  X* BMAG 

504  FORMAT 14H  BR  (  ,  Fb  .3*  2H  J  =  , E  16. 8, 10X, 3HB I ( • F6. 3, 2H I* »E 16. 8»1 OX, 

L  5HBMAGI ,F6.3 ,?H)=,E16.8) 

IFINXX.GE.NXJGO  TO  505 
•  NXX  =  NXX ♦! 

GO  TO  506 

505  PRINT  6 
STOP 
END 
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